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The absorption spectra of cyclohexene, cyclopentene, 1-3 cyclohexadiene and benzene have 
been examined in vapor phase in the region A = 2300-1600A. The electronic excitation responsible 
for the first group of bands in this region is practically independent of the cyclic structure of the 
molecule since the wave number of the first band agrees closely with that obtained for olefins 
having the same arrangement of hydrogen-atoms around the double bond. Characteristic 
vibrational frequencies and different electronic levels have been identified. 





ARLIER papers in this series have reported 

absorption spectra measurements between 
\=2300A and 1600A for sixteen ethylenic hydro- 
carbons! and the present paper extends this 
investigation to include three cyclic hydro- 
carbons—cyclopentene, cyclohexene and 1-3- 
cyclohexadiene. The spectrum of benzene has 
been re-examined under conditions which give a 
better definition of the vibrational bands than 
was obtained in the earlier work and these 
results are used to compare with the measure- 
ments for the other cyclic hydrocarbons. 


EXPERIMENTAL 


Cyclohexene was prepared by dehydration of 
cyclohexanol? according to the method outlined 


1E. P. Carr and M. K. Walker, J. Chem. Phys. 4, 751 
(1936); E. P. Carr and H. Stiicklen, Zeits. f. physik. 
Chemie B25, 57 (1934); J. Chem. Phys. 4, 760 (1936). 

? The preparation and purification of all the compounds 
used for spectrographic examination are under the direction 
of Professor Mary L. Sherrill of this laboratory. The 
purification of the cyclic hydrocarbons has been made by 
fractionations through a 35 cm Fenske column. The purity 
of the material is tested not only by the usual physical 
constants—boiling point, index of refraction and density 


by Kistiakowsky and co-workers’ and had a 
boiling point of 82.70°+.05 at 760 mm; an index 
of refraction, n?°D 1.4465; a density, d,?° 0.8101. 
A part of the cyclohexene was converted into 
1—3-cyclohexadiene as described by these authors. 
The diene used for the spectrographic work was a 
highly purified product obtained directly from 
the dibromide and its physical constants, b.p. 
80.13°+0.05°, n®°D 1.4742, d2° 0.8413, were in 
satisfactory agreement with the cyclohexadiene 
of Kistiakowsky prepared by way of the cyclo- 
hexenyl ether. 

Cyclopentene was prepared from cyclopenta- 
none which in turn had been prepared from 
adipic acid according to the method used by 
Kistiakowsky and co-workers.* Cyclopentanone 


but also by absorption spectrum measurements of the 
liquid hydrocarbon and its hexane solution with the quartz 

spectrograph. Cyclopentene was prepared by Elizabeth S. 
Matlack of this laboratory; cyclohexene and cyclohexa- 

diene by Kate Wilkins. This latter work was done in 
partial fulfillment of the requirements for the M.A. degree 

at Mount Holyoke College. 

( 038) B. Kistiakowsky et al., J. Am. Chem. Soc. 58, 140 
1936). 

a B. Kistiakowsky et al., J. Am. Chem, Soc. 59, 832 
1937). 
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TABLE I. purified by distillation. The resulting cyclo- ba 
pentene had a boiling point of 44.02°+.02°, gr 
eet ecnscraas nD 1.4420, dz? 0.7716. th 
= te -” Benzene, Kahlbaum special for molecular sti 
“— “— pooped weight determinations, was fractionated and the de 
48010 49245 430 middle fraction boiling at 80.0-80.1° was used ba 
260 510 730 for spectrographic examination. pl 
51100 52240 53680 The apparatus and experimental procedure for {rc 
— 4 5 aoee the absorption spectra measurements have been be 
465 860 630 described in detail.° For the earlier study of the 
ot want 55020 benzene a 12 cm absorption tube was used; all pre 
870 320 of the measurements reported in this paper were 
52100 560 made with a 60 cm tube. The exposure time 
55435 57820 59100 varied from 4 min. to 10 min. and all of the 
990 58220 270 ‘ar 
56990 670 610 Taste III. a 
57350 970 60160 a tiol 
10 
. BENZENE clas 
1/r 1/r 1/r ear 
TaBLe II. 48800 50320 52480 pri 
49035 570 670 165 
1-3 CYCLOHEXADIENE 205 850 940 as 
355 51030 53220 ah 
1/r ifr ifr 670 300 475 sho 
80-479 8560 000 780 690 700 ; 
43070-48210 ' $80 100 ’ 970 870 54040 mn 
48235-48365 790 200 50140 52030 430 the 
— - _ 55330 56740 57550 =e 
9400 9870 0410 420 830 635 qua 
480 O70 90525 772 57020 735 a 
580 50100 635 860 095 840 
680 200 765 56130 260 955 a st 
760 305 215 350 58120 h rd 
505 220 yd 
53650 590 1 
820 aaa 
51185 52605 950 55420 56220 57115 — 
280 700 54355 of te 
380 810 440 in tl 
500 900 540 TABLE IV. 7 
615 53000 650 carb 
710 100 730 aid 
820 200 820 CYCLOPENTENE 
945 310 940 exce 
52050 430 55045 oe oe ee all 
500 540 47603 48348 49168 T: 
720 490 285 ab 
57875 58125 58415 815 645 400 of tl 
960 220 510 935 772 520 ¥ 
58035 315 48062 915 653 with 
204 49040 788 vere 
900 
subd 
was reduced, using a copper-chromium oxide 50090 50870 51470 in st 
raat pier Be Iti , | had 320 51200 645 \e 
catalyst,® and the resulting cyclopentanol had a 405 60 825 As 
boiling point of 140-141° at 759 mm. The cyclo- 700 310 a band 
pentanol was dehydrated by means of fused their 
potassium acid sulphate, dried over lime and — = 55610 the § 
tive 
5 This reduction was kindly carried out for us at the . i 
University of Wisconsin through the courtesy of Dr. ————— ure 
Homer Adkins. 6 E. P. Carr and H. Stiicklen, reference 1. Was 1 
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bands which are given in the tables were photo- 
graphed while the vapor was kept flowing 
through the absorption cell. When the vapor is 
stationary or flowing very slowly photochemical 
decomposition takes place and new absorption 
bands appear. With each compound a series of 
plates were taken after irradiation of the material 
from 10-45 min. in order to distinguish clearly 
between the bands of the original material and 
those due to the presence of disintegration 
products. 


RESULTS 


The spectra of the cyclic unsaturated hydro- 
carbons show the same general regions of absorp- 
tion as were evident with the ojiefins and may be 
classified into regions A—D as was done in the 
earlier work. The present study is concerned 
primarily with the absorption between 2300 and 
1650A which is the region previously discussed 
as region C, the long wave portion and B, the 
short wave. This subdivision is quite arbitrary 
and was made to simplify the discussion. Each of 
the cyclic hydrocarbons has been examined also 
in pure liquid and hexane solution phase with the 
quartz spectrograph, and these curves will be 
published in a separate paper in connection with 
a study of photochemical decomposition of these 
hydrocarbons. 

The results of measurements with the fluorite 
vacuum spectrograph are presented in the form 
of tables and graphs. The most striking difference 
in the spectra of the aliphatic and cyclic hydro- 
carbons in this region is the presence of sharp, 
narrow bands in the latter, whereas, with the 
exception of cis- and transbutene-2 and isobutene 
all of the olefins show only diffuse bands. 
Tables I-IV give the wave numbers of the center 
of the absorption bands, which can be measured 
with an error of less than 10 cm~ for intense 
narrow bands. The values for each compound are 
subdivided into groups depending upon similarity 
in structure of the bands. 

As in the previous work, the intensities of the 
bands as given in the curves are estimated from 
their appearance on the photographic plate since 
the Schumann plates do not admit of quantita- 
tive photometric measurement. The vapor pres- 
sure of the hydrocarbon in the absorption tube 
was varied from 40 mm to less than 0.1 mm. 


The wave number at which absorption begins 
in the Schumann region was found for the 
olefins to be dependent upon the number of 
hydrogen atoms bound to the carbon atoms of 
the double bond. All of the cyclic compounds 
which have been studied, including benzene and 
diphenyl, have the grouping C—CH=CH—C 
and the analogous group for the olefins, R—CH 
=CH—R, was designated as IIa. For this 
group the beginning of absorption for the olefins 
was between 48,200 and 49,200 cm; for benzene 
this is at 48,800 cm and for the other cyclic 
hydrocarbons it is between 47,400 and 47,900 
cm~!. While there is a slight shift toward the 
visible of the first absorption band of all of the 
cyclic compounds except benzene in comparison 
with the olefins the close parallelism of the 
aliphatic and cyclic compounds is of interest. 
Evidently the energy of this electronic excitation 
which is common to all of the hydrocarbons 
containing the C=C linkage is not only in- 
dependent of the nature of the alkyl substituent, 
as was shown previously, but it would seem also 
to be practically independent of whether the 
molecule has a cyclic or an open chain structure. 

The variation in breadth of the narrow bands 
which are characteristic of these cyclic com- 
pounds is no doubt of significance but it seems 
premature to attempt any theoretical explana- 
tion. In the first group of bands of cyclohexene 
and cisbutene-2 the breadth of the bands is of 
the order of 100 cm~, while in transbutene-2, 
cyclopentene and cyclohexadiene it is from 40—50 
cm~!. In benzene these bands have broadened to 
such an extent that they overlap and give the 
effect of shallow diffuse bands and only in a 
narrow range of pressures are the single maxima 
sharp enough to be measured. In the short wave 
region benzene has a doublet band system in 
which the average width of a band is 40 cm™ 
and cyclohexene shows a comparable band 
system at somewhat longer wave-lengths. 


DIscUSSION OF RESULTS 
The first absorption bands in the Schumann 
region for ethylenic derivatives were ascribed to 
an electronic excitation corresponding to the 
transition 'A,—'B,'’ which had been predicted 


7K. P. Carr and H. Stiicklen, J. Chem. Phys. 4, 760 
(1936). 
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Fic. 1. 1,2, Cyclohexene; 3,4, cyclohexadiene; 5,6, benzene; 7, cyclopentene. 


for the ethylene molecule.’ Since the excitation 
involves an (x) orbital, localized in the region of 
the C=C bond, it would be reasonable to expect 
the vibrational bands which are characteristic of 
this linkage. Because of the breadth and diffuse- 
ness of the bands in the spectra of the higher 
olefins it was not possible to identify any of the 
vibrational bands with certainty but the appear- 
ance of fine structure in the butenes and in the 
cyclic compounds enables the measurement of 
band separations with an error of not more than 
10 cm—. Vibrational frequencies of 1510 cm™ 
for the excited molecule of cisbutene-2, 1570 cm 
for transbutene-2, and 1485 cm for cyclohexene 
have already been pointed out® together with 
their relation to the Raman frequencies of 1658 
and 1674 for cis- and transpentene-2 and 1646 
cm for cyclohexene. The spectrum of cyclo- 
pentene in the Schumann region shows only a 
long series of equally spaced narrow bands of 
gradually increasing intensity with increasing 
wave number, whereas in the other cyclic com- 
pounds there are groups of bands with sharply 
marked maxima (Fig. 1). In cyclohexene and 


8 Mulliken, Phys. Rev. 41, 751 (1932); J. Chem. Phys. 3, 


517 (1935). 
9E, P. Carr and H. Stiicklen, J. Am. Chem. Soc. (in 


publication). 


cyclohexadiene the separation between the most 
intense band in the first and second groups 
(Fig. 2) is of the order of that which is charac- 
teristic for this vibration in ethylenic compounds 
and seems, without question, to be related to the 
vibration along the axis of the C=C bond. 
In the spectrum of benzene (Fig. 2) the corre- 
sponding separation is 935 cm~! and this fre- 
quency must be ascribed to the C=C vibration 
which constitutes the so-called ‘“‘breathing”’ 
frequency of the benzene molecule and corre- 
sponds to the 992 cm~ in the normal molecule. 

In cyclohexene and cyclohexadiene there is a 
strong band, the separation of which from the 
strongest or 0—O band of the first group is of 











TABLE V. 
CoMPOUND 1/X RAMAN* 1/X ULTRAVIOLET 
Cyclohexene 1646 cm“ 1485 cm“ 
2912 2670 
3024 
1-3-Cyclohexadiene 1574 1360 +507 
3042 2920 
Benzene 992 935 
3060 2835 








* K. W. Kohlrausch and R. Seka, Ber. A68, 528 (1935). 

+ There are two bands of almost equal intensity in the center of the 
second group. In measuring the separation of the maxima of different 
groups an error of +50 cm™ is therefore possible. 
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the order of 2700-2900 cm. This band is 
present in the benzene spectrum also where it is 
one of the intense bands in the fourth group and 
has been found in cisbutene-2 in the third group. 
A second member of the C=C vibrational series 
is clearly separated from it in all of the com- 
pounds except cyclohexadiene. This band could 
be ascribed to a C—H stretching vibration, the 
frequency of which for hydrocarbons in their 
normal state is known to be of the order of 
3000 cm~!. Considering intensities, this assign- 
ment would be allowed. The band is equal or 
stronger than the neighboring higher member of 
the C=C vibrational series, while this again is 
weaker than the preceding member of its series. 
A comparison of these vibration frequencies for 
the excited molecule with the corresponding ones 
for the normal molecule from the Raman spectra 
is given in Table V. 

The absorption spectrum of cyclopentene 
(Fig. 3) is unique among the twenty-two hydro- 
carbons that have been examined in this labora- 
tory. The beginning of absorption at 47,600 cm— 
is in accord with that for the other cyclic hydro- 
carbons and the ethylenic derivatives of group 
Ila but the spectrum in this region shows only a 
series of about twenty equally spaced narrow 
bands, the average separation of which is 130 
cm. The intensity of these bands increases 
gradually but there is no evidence in this elec- 
tronic level of the discrete bands which have 
been ascribed to vibrations of the C=C bond 
in other molecules. It would seem as if an elec- 
tronic excitation which is analogous to that 
found for the other unsaturated hydrocarbons is, 
in the case of cyclopentene, unaccompanied by 
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the characteristic C=C vibration. It should be 
noted, however, that quite clearly marked band 
separations of 1520 and 2470 cm~ are shown in 
the higher electronic level which begins at about 
53,000 cm— although it has not been possible to 
identify corresponding vibrational bands in this 
electronic trinsition for other hydrocarbons 
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except isobutene. A comparative table of Raman 
spectrum measurements for cyclopentene by 
various investigators is given by Reitz.!° There 
is a strong Raman line at 1610 cm™ and a group 
of lines at 2846, 2902, 2947 and 3060 cm~, all 
of high intensity. 

The doublet system in cyclopentene which 
appears at 50,000 cm~ is undoubtedly due to 
the presence of cyclopentadiene. This band 
system appears very faintly when the ordinary 
exposures of from 4 to 6 minutes are made but 
comes out very strongly on longer irradiation. 
The curve for cyclopentadiene given by Scheibe 
and Grieneisen" is reproduced below that for 
cyclopentene (Fig. 3) and shows quite clearly 
that this system in our measurements of cyclo- 
pentene must be due to the diene which has been 
formed by photochemical decomposition. Further 
discussion of the cyclopentene spectrum must be 
postponed until the spectrum of cyclopentane 
and cyclopentadiene can be examined. 

From the general character of the bands in 
cyclohexene and their intensity relationships 
there would seem to be from 4-5 different 
electronic levels in the region between 47,000 
and 61,000 cm-. The first excitation beginning 
at 47,750 cm—! seems to be completely analogous 
to that in the ethylenic derivatives which has 
been ascribed to a 'A,—>'B, transition. The 
curves for cis- and transbutene-2 are given in 


10 A. W. Reitz, Zeits. f. physik. Chemie B33, 182 (1936). 
11 Scheibe and Grieneisen, Zeits. f. physik. Chemie B25, 
52 (1934). 
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Fig. 4 in order to show the striking parallelism 
between the spectra of cisbutene-2 and cyclo- 
hexene in the region between 47,000 cm and 
51,000 cm-!. These results have served as a 
confirmatory test for the assignment of molecular 
configurations to the isomers of butene-2.° The 
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only other electronic excitation in this spectrum 
which shows a marked resemblance to anything 
which has been found in the spectra of other 
unsaturated hydrocarbons is the doublet system 
which begins about 52,000 cm~ and extends to 
55,000 cm, the last four groups not being 
resolved. This system is quite similar to that 
found in benzene between 55,500 and 58,500 
cm~'; in each case the short wave component of 
the doublet is much stronger than the other; 
the separation between the two bands in the 
doublet is of the order of 100 cm in benzene 
and 130 cm in cyclohexene. In both spectra the 
separation between successive doublets seems to 
remain fairly constant although the separation is 
somewhat larger in cyclohexene (350 cm) 
than in benzene (250 cm~). 

In cyclohexadiene a series of more than fifty 
narrow bands, with an average separation of 
100 cm extends with but slight interruption 
from 48,500 to 55,000 cm. Underlying this 
series of bands are the diffuse bands whose 
separations correspond to vibrational frequencies 
of the C=C and C—H bonds. The narrow band 
spectrum reappears again in another electronic 
level in the region of 58,000 cm but the separa- 
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tion here seems to be slightly less than in the 
longer wave region. It is not possible at this 
time to attempt any theoretical discussion of 
these narrow bands which seem to be so charac- 
teristic a part of the spectra of cyclopentene and 
1—3-cyclohexadiene. 

In the earlier discussion of the benzene 
spectrum” very little could be said about the 
region between 48,500 and 54,000 cm™ since 
the experimental conditions were not such that 
the narrow bands in this region were resolved 
and with the exception of the first two bands 
where some splitting was evident only the 
maxima of groups of diffuse bands were given. 
The results of the present investigation of the 
benzene spectrum show that the strongest bands 
of each of the first six groups, as indicated in 
Fig. 2, fit into a vibrational series of the form: 
1/A= 49035 +9607 — 25v’, the average deviation 
of +3 cm being within the limits of experi- 
mental error. The lower intensity of the first, or 
0—0O band, would be explained by the Franck- 
Condon principle ; after the second the intensity 
decreases as must be expected for a vibrational 
series. In the case of cyclohexene and cyclo- 
hexadiene it is of interest to note that in each 
spectrum the first two bands have the same 
intensity. 

The general characteristics of the spectra of 
benzene, cyclohexene, cyclohexadiene in the 
region beginning about 48,000 cm~ and their 
analogy to the spectra of the ethylenic deriva- 
tives point to the conclusion that the electronic 
transition involved is similar in all of the 
molecules. In the ethylenic hydrocarbons this 
transition was tentatively assigned to the transi- 
tion '4,—'B, and the much weaker absorption 
in the quartz region to a 'A,—>*B, transition. 
The olefins show a step-out indicating a broad 
absorption band of very low intensity in the 
quartz region so that no study of fine structure 
is possible but the analogy to the benzene 
spectrum is shown by the fact that the absorption 
of the group II hydrocarbons in liquid and 
hexane solution is in precisely the same region as 
that of benzene solutions."* The sharp bands of 


2E. P. Carr and H. Stiicklen, Zeits. f. physik. Chemie 
B25, 57 (1934). 

13 E, P. Carr and G. F. Walter, J. Chem. Phys. 4, 760 
(1936); International Critical Tables, Vol. 5, p. 363. 
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the latter spectrum are well known and the vapor 
phase spectrum has been exhaustively investi- 
gated. If, therefore, the assumption that closely 
corresponding electronic transitions are involved 
in the excitation of all the molecules having the 
—C—CH=CH—C-— grouping is valid the ab- 
sorption spectrum of benzene in the long wave 
region between 37,000 and 44,000 cm! would be 
ascribed to a singlet-triplet transition. Some 
justification for this assumption is given by the 
recent interpretation of the fluorescence spectrum 
of benzene by Ingold and Wilson" in which a 
complete analysis of the spectrum was possible 
on the basis of transitions from three excited 
electron states to the zero level of the ground 
state. 


4C, K. Ingold and C. L. Wilson, J. Chem. Soc. 941 
(1936). 
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The vapor pressure of the benzene was too 
high and the absorption tube too short to resolve 
the benzene spectrum sufficiently to permit 
detailed analysis of the bands but the general 
characteristics of the spectrum between 48,800 
and 54,400 cm~ and in the quartz region are 
sufficiently similar to indicate closely related 
electronic transitions. The probability of the 
singlet—triplet transition for ethylene was esti- 
mated by Mulliken to be of the order of 10-4 of 
that for the singlet—singlet transition and the 
ethylenic derivatives showed close agreement 
with that prediction. If the assumption of 
corresponding states for benzene is correct, com- 
parison of the intensities in the two regions 
would show that the probability of the transition 
to the triplet state in benzene is about 10~ of 
that for the singlet—singlet transition. 
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Infrared Absorption Studies 


IV. The Infrared Spectrum of Water in an Inert Solvent 


L. B. Borst, A. M. BusweELL AND W. H. RopEBuUsH 
University of Illinois, Urbana, Illinois 
(Received November 22, 1937) 


The infrared absorption spectrum of HO, HDO and D.O in dilute solution in carbon tetra- 
chloride and carbon bisulfide has been investigated with a grating of high resolvirg power, 
in the region 2.6—3.8u. Several bands have been identified as v3; bands and the hitherto unob- 
served v; band for H,O is believed to have been found. The frequencies observed are compared 
with the same frequencies in the vapor and a displacement of frequency which is proportional 
to the frequency itself is observed in solution. A fine structure is observed in the case of H,O ~ 
which may be ascribed to rotation unless some other explanation is found. The spectrum in 
carbon bisulfide appears to be identical with that in carbon tetrachloride. 


HE infrared vibraton rotation spectrum of 

water vapor was investigated by Plyler and 
Sleator! and that of heavy water vapor has been 
recently studied by Barker and Sleator.? They 
have located the principal vibration frequencies 
for water and the two deuterium substituted 
molecules with the exception of »; which has 
not been observed in the symmetrical molecule 
and v3 for HDO. »;, however has been observed 


as ae Plyler and W. W. Sleator, Phys. Rev. 37, 1493 
on Barker and W. W. Sleator, J. Chem. Phys. 3, 660 
5). 





in the Raman spectrum for H.O.* Studies of 
hydrogen chloride absorption in solution have 
been made by Plyler and Williams,‘ and West 
and Edwards.® Recently Kinsey and Ellis® have 
observed the infrared absorption spectrum of 
water in carbon bisulfide solution with a prism 


spectrograph. 


3D. Bender, Phys. Rev. 47, 252 (1935). 

4E. K. Plyler and D. Williams, Phys. Rev. 49, 215 
(1936). 

5W. West and R. T. Edwards, J. Chem. Phys. 5, 14 
(1937). 

6E. L. Kinsey and J. W. Ellis, Phys. Rev. 51, 1074 
(1937). 













































It seemed desirable to the authors to investi- 
gate the spectrum of water in an inert solvent 
with a spectrograph of high resolving power in 
order to determine what differences were pro- 
duced in the spectrum by the presence of the 
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Fic. 1. Vapor absorption bands for mixtures of H:O and 
D.O. Prism spectrometer. 
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solvent molecules. A grating spectrograph which 
was purchased from the University of Michigan 
was available for this purpose. This instrument 
has been described.? The grating was of the 
echelette type, 3600 lines to the inch, ruled on 
pure tin. It was constructed for us in the physics 
shop of the University of Michigan under the 
direction of Professor Barker. As a light source 
we used an ordinary street lighting bulb with a 
coiled tungsten filament of 6.6 amperes rating. 
The thin glass walls transmit very satisfactorily 
up to about 4y and this type of lamp is much 
superior to the Nernst glower. The utmost 
precautions were taken in purifying the solvents 
used in preparing the solutions. As a check on 
the calibration of the grating, the absorption 
spectrum of vapor was recorded. The results 
agreed within about 3 cm~ with those of Plyler 
and Sleator. 

In measuring the absorption, a cell 7.5 cm in 
length with two compartments, one for pure 


"TW. W. Sleator, Astrophys. J. 48, 125 (1918). 
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solvent and one for solution, was used. When the 
absorption of the vapor was measured, the cells 
were heated by circulating water from a ther- 
mostat. In representing the data, wave numbers 
are plotted as abscissa and log J)/J is plotted as 
ordinate, where J is the galvanometer deflection 
through the pure solvent and J, the deflection 
when the solution is in the path. 

The vapor absorption spectra of HO, HDO, 
and D.,O were mapped for comparison and 
identification purposes with a prism spectrometer 
of low resolving power (Fig. 1). 

When a survey was made of the solution with 
the same instrument there was observed a band 
at 2.74 due to H.O and three bands between 
3.5 and 3.8u which were attributed to HDO 
and D.O. These are shown in Fig. 2. The carbon 
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Fic. 2. Absorption bands of carbon tetrachloride saturated 
with a mixture of H.O and D.O. Prism spectrometer. 


tetrachloride was saturated with water contain- 
ing varying percentages of DO. It is evident 
that there is no very close correspondence 
between the absorption in the vapor and in 
solution. The main absorption peaks are un- 
doubtedly due to the antisymmetric vibration ;. 
It is suggested that the additional small peaks in 
solution, which have no counterpart in the 
vapor, might be due to the symmetric vibration 
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v; which has not heretofore been observed in the 
infrared. 

In Fig. 3 the absorption curves are given for 
H.O in carbon tetrachloride solution as obtained 
with the grating spectrometer. The upper curve 
is obtained from a 12.5 cm layer saturated with 
water. (The solubility of water in carbon tetra- 
chloride is stated to be 0.0083 mole per liter.) 
The lower curve is obtained by diluting the 
solution to one-fourth of the original concentra- 
tion. The fine structure obtained here is believed 
to be real as it was checked repeatedly. In order 
to eliminate the possibility of water vapor in the 
system, the curves shown in Fig. 4 were obtained. 
The upper curve is the absorption curve for 
water vapor and the lower curve is for 5 cm 


| | | 
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Fic. 3. H,O in carbon tetrachloride. Grating spectrometer. 





3760 3700 3600 em“! 





Fic. 4. Comparison of vapor and solution of H,O. 
Grating spectrometer. 


saturated carbon tetrachloride solution. There is 
no correlation observable. 

Figure 5 is the absorption curve in the 2.7u 
region for D,O. The negative absorption simply 
indicates the presence of traces of water in the 
supposedly dry pure solvent. In spite of all pre- 
cautions some water was present. In the cell to 
which D,O was added this water was effectively 
converted to HDO. The peak is therefore to be 
attributed to v3 for HDO, the minima to the H,O 
in the solvent. There is evidence of fine structure. 

In Fig. 6 are given the absorption curves for 
three concentrations of D.O in the region 3.4 
to 3.9u. The upper curve (dotted line) was 
obtained from a sample saturated with 99.6 per- 
cent D,O. The middle curve was for a solution 
saturated with 85 percent DO and the lower 
curve was obtained upon diluting this solution 
one-half. The breaks in the lines are due to 
readjustment of the monochromator. These 
curves were obtained using the first order because 
of low intensity. The resolving power is therefore 
only one-half that used in obtaining the curves 
in Figs. 3, 4, and 5. 

The identification of certain of the peaks seems 
obvious enough by comparison with the work of 
Sleator and colleagues. The frequency is shifted 
in solution to lower values, the shift being 
apparently roughly proportional to the fre- 
quency. Table I gives the principal absorption 
peaks. The number in the first column serves to 
locate the peak on the accompanying curves. 
The second column gives the frequency observed 
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Fic. 5. Solution of D,O in carbon tetrachloride. 
Grating spectrometer. 





































by us in carbon tetrachloride solution. The 
third column gives the frequency as observed in 
the vapor either by infrared absorption or in the 
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Fic. 6. Solution of D.O in carbon tetrachloride. 
Grating spectrometer. 


Raman spectrum. The fourth column gives the 
difference between solution and vapor. This 
difference is seen to be roughly proportional to 
the frequency and this can be justified on 
theoretical grounds. [y;] is enclosed in brackets 
since this frequency has not been observed before 
in the infrared and the assignment is therefore 
necessarily tentative since there is no very 
obvious reason why it should be more active in 
solution than in the vapor state. [v3] has not 
been observed for HDO either in the Raman or 
infrared spectra and the calculated value for the 
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TABLE I. Frequencies (cm™) of the principal 
absorption peaks. 




















SOLUTION VAPOR DIFFERENCE 
1 2640 [v,] D,O 2666 26 
2 2689 [»,] HDO 2720 33 
3 2748 v3 DO 2784 36 
4 3611 [»,] H,O 3654 43 
5 3659 [»;] HDO [3750] [91] 
6 3702 vs; HO 3756 44 








vapor [3750] may be in error by 50 wave 
numbers. 

There seems to be no doubt as to the reality 
of the fine structure apparent in Fig. 3 and in 
the absence of other explanations it must be 
assumed to be due to some sort of submerged 
rotational structure. This fine structure was not 
observed in the deuterium substituted com- 
pounds for the reason that, assuming it to be a 
rotational structure, a higher resolving power 
would be required, whereas these spectra were 
observed in the first order with only half the 
resolving power used for ordinary water. 

The absorption spectrum of water in solution 
in carbon bisulfide was observed but the curve is 
so nearly identical with those published here for 
the carbon tetrachloride solutions that there 
seems to be no reason for reproducing it. 

The authors wish to acknowledge their in- 
debtedness to Dr. J. G. Kirkwood and Dr. 
F. T. Wall for helpful discussions in regard 
to the interpretation of this work. They wish 
also to acknowledge their indebtedness to the 
Rockefeller Foundation for financial support. 
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A Spectral Investigation of Glucose Glass 


WitL1 M. Conn 
Berkeley, California 


(Received November 12, 1937) 


The spectral transmission of glucose glasses has been investigated over the wave-length 
region from 3600 to 8000A. Glasses of different thermal history were tested, ranging from white 
to dark brown in color. The glasses show different transmission values, and shifts of the trans- 
mission maximum and the absorption limit at small wave-lengths with increasing caramelization. 





I. INTRODUCTION 


HE thermal properties of glucose glass have 
been investigated by Parks andco-workers.! 

The color of these glasses varies, depending on 
the thermal history, from white to dark brown 
indicating the formation of decomposition prod- 
ucts during caramelization. (It is generally 
assumed that small amounts of colloidal carbon 
are finally liberated.) An investigation seemed of 


1G. S. Parks and H. M. Huffman, J. Phys. Chem. 31, 
1842-55 (1927); G. S. Parks, H. M. Huftman and F. R. 
Cattoir, J. Phys. Chem. 32, 1366-79 (1928); F. R. Cattoir 
and G. S. Parks, J. Phys. Chem. 33, 879-882 (1929); 
G. S. Parks and W. A. Gilkey, J. Phys. Chem. 33, 1428-37 
(1929); G. S. Parks and S. B. Thomas, J. Am. Chem. Soc. 
56, 1423 (1934); G. S. Parks, L. E. Barton, M. E. Spaght 
and J. W. Richardson, Physics 5, 193-199 (1934); G. S. 
Parks, L. J. Snyder and F. R. Cattoir, J. Chem. Phys. 2, 
595-598 (1934); J. C. Lyman and G. S. Parks, J. Chem. 
Phys. 4, 218-219 (1936); G. S. Parks and J. D. Reagh, 
J. Chem. Phys. 5, 364-367 (1937). 





interest to determine the effect of caramelization 
on the spectral transmission of glucose glass. 


II. METHOD 


Spectra were obtained with my Hilger grating 
spectrograph, dispersion 49.4A/mm. The range 
from 3600 to 8000A was investigated using 
Wratten and Wainwright hypersensitive pan- 
chromatic plates and Eastman Kodak I-N plates. 
The first-order spectrum was used, the over- 
lapping ultraviolet of the second order being 
eliminated by a No. 12 Wratten gelatin filter. 
A straight filament 8.5-volt lamp with rheostat 
control was used as a light source. A condenser 
lens projected the light on the slit. The slit 
width was 0.1 mm. A very fine neutral wedge 
from Hilger was placed in front of the slit. 
(A separate investigation showed that this wedge 


GG2 


Fic. 1, Plate No. BL 401. Hyperpan plate. Wedge spectra for lamp alone, and glucose glasses GG 1 and GG 2. Time of 
exposure 2 minutes. 
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Lamp 


Fic. 2. Plate No. BL 410. I-N plate. Wedge spectra for lamp alone, and 
glucose glass GG 1. Time of exposure 10 minutes. 


is neutral over the entire wave-length region 
investigated.) The density of the wedge was 
such that the light end transmits 10,000 times 
the light transmitted by the dark end.? A wave- 
length and density scale was printed on the plate, 
together with the spectrum. (The wave-length 
scale needs a slight correction which was deter- 
mined from line spectra.) 

The glucose glass samples were contained in 
metal cells of 23.5 mm inside diameter with 
silica glass windows, the length of cells GG i and 
2 was 24.0 and that of all other cells was 25.0 mm. 
The samples were placed in front of the wedge. 
Several spectra were taken from each sample. 
The difference in density of the lamp and glass 
spectra on a plate was small enough to permit 
the photometric reduction. 

The photometric reduction was carried out as 
follows: A wedge spectrum of the lamp alone 
was printed on each plate, together with one or 
two glucose glass spectra. All spectra on a plate 
were taken with the same time of exposure and 
the same voltage of the lamp. Fig. 1 shows the 
panchromatic plate BL 401 with wedge spectra 
of the lamp alone and glasses GG 1 and 2, and 
Fig. 2 shows the I-N plate BL 410 with wedge 
spectra of the lamp and glass GG 1. 

The density of the glass and lamp spectra was 
measured every 20 my. A simple computation 
gave then the light transmission values. Con- 
sidering the values for the lamp alone as unity 
the transmission of the glass could be computed 
in percent. The steep decrease in sensitivity of 
the panchromatic plates makes values around 
660 mu less accurate than the remaining values. 
(All values of less weight are given in parentheses 


2C. E. K. Mees and S. H. Wratten, Brit. J. Phot. 54, 
384-385 (1907). 


in Table II.) The probable error of the 
measurements amounts to +2 percent. 


III. SAMPLES 


The samples were prepared by Dr. Parks and 
co-workers at different times, samples GG 1 and 2 
were, approximately, four years old. These 
samples had been prepared from “C. P. Pfan- 
stiehl’’ anhydrous glucose with a subsequent 
purification by the method of Hudson and Dale’ 
with glacial acetic acid. The remaining samples 
were made directly from Pfanstiehl glucose with- 
out this special purification. The samples were, 
with the exception of GG 4 and 6, remarkably 
free of bubbles. No striae appeared in any glass. 
GG 4 and 7 showed the beginning of crystalliza- 
tion starting from the edge of one, front end. 
GG 1 was practically colorless, the remaining 
samples extended from light yellow to dark 
brown in color, see Table I. This table presents 
also colorimetric values obtained by means of a 
Lovibond tintometer, series No. 52.4 Glucose 
glasses GG 3 and 5, GG 4 and 6, and GG 7 and 8, 


TABLE I. Color and Lovibond values for glucose glasses. 











LovIBOND 
Grass No. CoLor VALUE 
GG1 White 3.0 
GG 2 Light yellow 11.0 
GG 7 Yellow 12.0 
GG 8 Yellow 10.7 
GG 9 Yellow 12.0 
GG 4 Yellow-brown 15.3 
GG 6 Yellow-brown 14.7 
GG 3 Dark brown 37.3 
GG5 Dark brown 36.7 








3C. S. Hudson and J. K. Dale, J. Am. Chem. Soc. 39, 


324-328 (1917). 
4K. S. Gibson and F. K. Harris, Sci. Papers, Nat. Bur. 
Stand. No. 547, 1927. 
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respectively, show only small differences in color, 
each of these pairs having been prepared at 
about the same time. 


IV. RESULTS 


Table II contains the results of the individual 
measurements of the glucose glasses (reduced to 
a thickness of 25 mm). Fig. 3 shows the trans- 
mission curves including averages for glasses 
GG 3 and 5, GG 4 and 6, and GG 7 and 8. 
The difference in transmission values of GG 4 
and 6 is due to the larger number of bubbles and 
the beginning of crystallization in GG 4. The 
following remarks refer to the results obtained. 

(1) A careful study under varying conditions 
showed the absence of narrow absorption bands 
in the entire wave-length region investigated. 

(2) The spectral energy distribution shows flat 
maxima for all glasses. The maxima occur in the 
range from 580 to 640 my. The maxima of 
transmission are shifted toward greater wave- 
lengths with increasing caramelization. 

(3) The absorption limits at small wave- 
lengths are shifted toward the red with in- 
creasing caramelization. Table III gives the 
absorption limits for small wave-lengths. 

(4) None of the glasses shows a limit of trans- 
mission toward great wave-lengths within the 
range investigated (up to 800 mu). The trans- 
mission curves become, at great wave-lengths, 
approximately parallel to the abscissa for several 
glasses; a small increase in transmission is indi- 


cated for GG 8. 
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21G. 3. Spectral transmission of glucose glasses. 


(5) The transmission values of GG 1 are lower 
than those of the remaining glasses. This is 
probably due to the age of GG 1 as against the 
recent glass batches GG 3 to 9. It seems likely 
that crystallization has started already in GG 1 
and 2, although it is not visible. Another cause 





TABLE II. Transmission of glucose glasses in percent as a 
function of the wave-length; thickness of layer 25 mm. 
































WAVE- 
LENGTH | GG 1| GG 2| GG 7| GG 8| GG 9| GG 4| GG 6| GG 3|GG 5 
IN Mz To % % % % % % % % 
440 8.4 
460 15.2 2.1 
480 22.9} 11.8 9.6 9.2 | 11.6 
500 29.1 | 29.3 | 20.4] 17.0 | 26.2 
520 33.6 | 53.3 | 29.6 | 25.4 | 44.4 | 13.2 | 10.2 
540 36.9 | 66.9 | 37.8 | 33.6 | 59.2 | 15.0 | 19.4 
560 39.0 | 69.0 | 44.8 | 41.8 | 70.6 | 17.2 | 29.8 3.0 | (3.2) 
580 40.1 | 70.6 | 49.6 | 49.6 | 78.6 | 21.6 | 37.8 7.6 6.6 
600 38.6 | 71.3 | 55.2 | 56.2 | 84.6 | 26.8 | 45.8 | 13.6 | 14.4 
620 34.6 | 62.9 | 59.2 | 60.0 | 86.4 | 27.2 | 50.8 | 22.2 | 21.8 
640 28.2 | 48.7 | 58.4 | 57.2 | 80.2 | 25.6 | 49.2 | 28.0 | 25.4 
660 (20.4) | (37.1)| (44.8) | (46.6) | (65.4) | (24.2) | (40.2) | (28.0) | (24.4) 
680 14.5 | 28.2 34.6 | 49.6 | 23.4 | 36.8 | 27.0 | 19.4 
700 11.0 | 20.7 27.8 | 37.2 | 18.6 | 33.8 | 25.6 | 16.6 
715 10.4 | 15.2 27.0 | 31.2 | 13.3 | 31.2 | 24.8 | 16.9 
735 13.2 28.3 | 24.6 | 11.6 | 26.9 | 24.7 | 15.7 



































TABLE III. Absorption limits in my at short wave-lengths for 
glucose glasses. 








GG1,/GG2,/GG7,|GG8,|GG9,|GG4,|GG6,|GG3,|GGS5, 


Giass No. mye | mp | mp] mp | mp] me] me] mp] mp 
Absorption 
Limit 397 | 426 | 447 | 441 | 441 | 479 | 481 | 539 | 539 








of this difference may lie in the different purifi- 
cation process of the two groups of samples. A 
recrystallization would tend to decrease the 
transmission values over the entire spectrum. 

The results indicate that the spectral trans- 
mission of glucose glass depends on the thermal 
history of the glass, in a similar manner as do 
the thermal properties. A spectral investigation 
presents, if carried out under standardized con- 
ditions, a quick and sensitive method for a 
classification of samples, and a checking of the 
thermal history of the samples. 


V. SUMMARY 


The spectral transmission of glucose glasses of 
different thermal history has been investigated 
between 3600 and 8000A. All glasses show flat 
transmission maxima which are shifted toward 
greater wave-lengths with increasing carameliza- 
tion. The absorption limits at small wave-lengths 
are shifted toward the red with increasing 
caramelization. 

In conclusion I should like to express my 
thanks to Dr. George S. Parks of Stanford 
University for preparing the samples and for 
helpful discussions I have had with him. 
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The Heat Capacity of Diamond from 70 to 300°K 


KENNETH S. PITZER 
Department of Chemistry, University of California, Berkeley, California 
(Received November 20, 1937) 


The specific heat of diamond has been measured over the temperature range 70 to 300°K, 
and is compared with the theoretical formula of Debye. Small but definite deviations are 
observed, which in both direction and order of magnitude are those to be expected when the 
approximations of the Debye formula are taken into consideration. The entropy of diamond 
at 298.1°K is 0.585+0.005 cal. per degree per mole. 





HE temperature variation of the heat capac- 

ity of crystals received a satisfactory gen- 
eral explanation shortly after the development of 
quantum statistics. The papers of Einstein,’ 
Born and v. Karman,? and Debye’ are particu- 
larly noteworthy, and it is the formula developed 
by the latter that has received most attention. 
Debye assumed that the frequencies of the 
normal modes of vibration of the crystal are 
distributed just as the vibration frequencies of a 
continuum, except that the spectrum is cut off 
at a certain maximum frequency (vp) so chosen 
that there are exactly the proper total number 
of frequencies (3N if N be the total number of 
atoms). 

TABLE I. Molal heat capacity of diamond. 








| 











5 Cp 6D 
iv | (cal./deg.) (Cp/T*) X 10% | (Debye) 
70.16 | 0.0224 650 1930 
75.37 * (0296 692 1890 
81.59 0362 667 1900 
88.65 0452 650 1930 
96.68 0592 655 1920 
105.10 0785 | 677 1900 
113.04 0991 686 1890 
125.28 138 | 703 1874 
134.29 | 174 | 719 1863 
144.10 | 218 | 730 1851 
153.71 | 267 | 736 1845 
162.76 | 318 | 740 1840 
173.33 | 386 743 1833 
181.96 | 444 | 739 1832 
191.44 | 518 740 1830 
200.94 595 | 735 1818 
211.84 685 | 722 1820 
231.06 858 | 694 1817 
241.09 918 653 1845 
252.37 1.033 642 1840 
264.31 1.137 | 616 1852 
276.61 1.257 595 1855 
287.96 | 1.355 | 568 1865 








1 Einstein, Ann. d. Physik 22, 180, 800 (1907); 34, 170, 
590 (1911). 

? Born and v. Karman, Physik. Zeits. 13, 297 (1912). 

* Debye, Ann. d. Physik 39, 789 (1912). 
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As seems likely in view of its derivation, the 
Debye formula only approximates the experi- 
mental data for most substances. Blackman‘ has 
discussed the nature of the deviations to be 
expected, and has shown that the specific heat 
curve for the vibrations of a given crystal will 
always approach the Debye formula at very low 
temperatures, but that it will be a rare and 
accidental occurrence if the Debye formula holds 
precisely over the entire temperature range. 

It was the purpose of this investigation to 
consider the heat capacity curve for a particular 
substance in the light of these discussions. The 
choice of diamond as the subject was favored by 
the high symmetry of the crystal, the presence of 
only one kind of atom, the fact that the inter- 
atomic forces are reasonably well known and of 
short range, and the experimental consideration 
that, in a comparative sense, very low tem- 
peratures are easily attained. 


EXPERIMENTAL 


About twenty grams of diamonds of a grade 
known as “‘bortz’’ were obtained for these experi- 
ments. The stones, while imperfect, were reason- 
ably clear and were free from foreign matter. It 
is believed that their heat capacity was the same 
as that of perfect diamonds within the experi- 
mental errors of this investigation. 

The calorimeter employed in this investigation 
was of similar character to that described by 
Latimer and Greensfelder® but was in size just 
large enough to contain the diamonds. The 
cryostat and other equipment were the same as 
has been used previously. The measurements 

* Blackman, Proc. Roy. Soc. London A148, 365, 384 
(1935); A149, 117, 126 (1935). 


’ Latimer and Greensfelder, J. Am. Chem. Soc. 50, 2202 
(1928). 
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were of similar accuracy except at the lowest 
temperatures where the heat capacity of the 
diamonds was only about one-tenth of the heat 
capacity of the empty calorimeter. The error in 
the results may be as large as 5 percent at the 
lowest temperatures but rapidly diminishes to 
about 0.2 percent at 150°K. At the highest tem- 
peratures the errors become somewhat larger 
due to the rapid heat exchange by radiation. 
The results are listed in Table I and are presented 
graphically in Fig. 1, together with those of 
other investigators. ® 

As can be seen from Fig. 1, the results of 
various investigations are in reasonably good 
agreement with the possible exception of the 
data of Nernst below 100°K. That these values 
may be in error seems not improbable when one 
considers that Nernst’s results at even lower 
temperatures show a smaller heat capacity for 
the filled calorimeter than for the empty. 

The entropy of diamond was calculated by a 
graphical integration of C, vs. log T, together 
with an extrapolation to the absolute zero based 
on the 7* law. The results are summarized in 


Table II. 
DISCUSSION 


Columns three and four of Table I give the 
values of (C,/T*) and @p (Debye) respectively 
for the various temperatures. The heat capacity 
of diamond is found to be proportional to the 
third power of the temperature up to about 
100°K, after which it rises faster than 7* for 
about 75°. This increase of the heat capacity, 
faster even than 7%, is of course unexplainable 
by the Debye theory but is quite possible in 


TABLE I]. The molal entropy of diamond. 








0 - 70.8°K 


T? extrapolation 0.0077 cal./deg. 
70.8-298.1° 


Graphical from data 0.5776 
Entropy at 298.1°K =0.585 +.005 








view of Blackman’s considerations of lattices. 
Similarly the values of 6p are constant at about 
1910 up to 100°K, then decrease gradually to 
about 1820 at 200°K, after which they appear 





° Nernst, Ann. d. Physik 36, 395 (1911); Weber, Phil. 
Mag. 49, 161 (1875); Robertson, Fox and Martin, Proc. 
Roy. Soc. London A157, 579 (1936). 
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Fic. 1. The molal heat capacity of diamond in cal. per 
degree per mole. Solid circles, this research; crosses, Nernst; 
triangles, Weber; open circles, Robertson, Fox and Martin. 


again to rise slightly. This variation of about 
5 percent is, of course, not very great but is 
entirely outside the experimental error. 

Some rough calculations made by the writer 
indicate that this variation of @p is entirely 
reasonable. These calculations followed the 
general methods of Blackman‘ but were specif- 
ically set up for the diamond lattice. The only 
forces considered were central forces between 
adjacent atoms and forces tending to keep the 
bond angles those of a tetrahedron. The ratio 
between the two force constants was estimated 
from the normal vibration frequencies of propane. 
The results of this calculation deviate from the 
experimental values by about the same amount 
as the Debye theory, but the deviations are in 
the opposite direction. The true distribution of 
vibration frequencies is thus somewhere between 
that given by the rough lattice calculations and 
that of the Debye “continuum.” It is now ap- 
parent that there is nothing anomalous about the 
specific heat of diamond, and that the assumption’ 


7Simon, Sitzber. preuss. Akad. Wiss., Physik. Math. 
Klasse 33, 477 (1926). 





of “thermal quantum jumps’’ is quite unnecessary. 

In conclusion it may be said that there is 
agreement between experimental and theoretical 
specific heats of crystals, within the uncertainty 
in the latter. This uncertainty is, nevertheless, 
quite large, and it is hoped that in the near future 
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the experimental data presented above can be 
compared with a more precise theory. 

The writer wishes to thank Professor W. M. 
Latimer for his assistance in obtaining the 
sample of diamonds and for his interest and 
advice throughout the investigation. 
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University of Chicago, Chicago, Illinois and 
Cornell University, Ithaca, New York 
(Received November 12, 1937) 


A theory of order and disorder in solid solutions, based upon a direct evaluation of the crys- 
talline partition function, is described. As in the theory of Bethe, the fluctuation in energy within 
states of a fixed degree of order, is taken into account. The method promises to be of use in the 


study of other cooperative phenomena. 


I 


T has been found that the atoms of certain 
solid solutions, although distributed at random 

among the lattice sites of the crystal at high 
temperatures, tend to arrange themselves into an 
ordered structure or superlattice at low tempera- 
tures. Moreover, there exists a critical temper- 
ature, analogous to the Curie temperature in 
ferromagnetism, below which the ordering process 
first begins to set in. Experimental evidence for 
the existence of order in such alloys as Cu-Au and 
Cu-Zn is provided by x-ray analysis of their 
crystal structure, and by the observation of 
discontinuities in certain of their physical con- 
stants or in their derivatives, which occur at the 
critical temperature. 

Bragg and Williams! have treated the problem 
of order and disorder from the _ theoretical 
standpoint. Their theory accounts for the exist- 
ence of a critical temperature and provides an 
approximate description of the influence of 
order on the thermodynamic functions of binary 
solid solutions. In the development of the theory 
they assume that the energy of the crystal is 
determined by its degree of long range order. 


* The outlines of the present theory were sketched by the 
author at the Symposium on the Structure of Metallic 
Phases, held by the department of physics of Cornell 
University, July 1, 1937. 

1 Bragg and Williams, Proc. Roy. Soc. London A145, 699 
(1934); 151, 540 (1935); 152, 231 (1935). 





It was pointed out by Bethe? that a given degree 
of order includes many crystalline configurations 
of different energy. Bethe developed a more 
elaborate theory in which the fluctuations in 
energy were taken into account. His theory, later 
extensively applied by Peierls,’ yields results 
which are more nearly in accord with experiment 
than does the theory of Bragg and Williams. 

We shall describe here an alternative theory of 
order and disorder, essentially equivalent to 
Bethe’s in its degree of approximation, but 
somewhat less unwieldy in form. It is based 
upon a straightforward evaluation of the parti- 
tion function of the crystal by the methods of 
statistical mechanics. The method of calculation 
bears a certain resemblance to that of the 
Heisenberg theory of ferromagnetism. However, 
the convergence difficulties, encountered by the 
ferromagnetism theory at low temperatures, do 
not arise, at least for solid solutions of rational 
composition. We believe that the theory may 
prove useful not only in the study of order and 
disorder in solids, but also in the treatment of 
cooperative phenomena in general. 


II 


We consider a binary solid solution containing 
equal numbers, N, of atoms of the components 


2 Bethe, Proc. Roy. Soc. London A150, 552 (1935). 
3 Peierls, Proc. Roy. Soc. London A154, 207 (1936). 
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ORDER AND DISORDER 


We suppose that the crystal is built up of two 
interpenetrating simple lattices, a and b, as for 
example, a body-centered cubic structure. Each 
lattice point of a has z nearest neighbor sites on 
lattice b, and each point of b has an equal number 
of nearest neighbor sites on a. In enumerating 
the microscopic states of the crystal we must 
consider the possible distributions of the N 
atoms of type 1 and the N atoms of type 2 
among the sites of the two lattices, a and b. We 
shall refer to each such distribution as a con- 
figuration of the crystal. We shall make the 
assumption that the partition function associated 
with lattice vibrations is the same for all con- 
figurations of the crystal. With this simplification 
the partition function may be expressed as the 
product of two factors, one arising from configu- 
ration and the other from lattice vibrations. 
Correspondingly, the thermodynamic functions 
split up into a configurational term and a lattice 
vibrational term. The latter may be calculated 
in the usual manner by the methods of Debye 
and of Born and Karman.‘ We shall be interested 
here only in the configurational contributions. 

The various configurations of the crystal may 
be partially characterized by a parameter, s, 
called the long range order, which is defined as 
follows. 


s=(N, —N,)/N, (1) 


where Ni, is the number of atoms of type 1 
situated on lattice a, and N, the number on 
lattice b. We remark that N,™ is equal to N2\, 
the number of atoms of type 2 on lattice 6, and 
that NV, is equal to N2™, the number of type 
2 on a. Further, since VN, +N, = N, we may 
write 
Ni, =N, = N(A+s)/2, 


N, = N2 = N(1—s)/2. (2) 


We notice that when s has the values, +1, the 
crystal is in a completely ordered state with all 
atoms of type 1 on one of the lattices and all 
atoms of type 2 on the other lattice. On the other 
hand, when s is zero, each type of atom is 
randomly distributed between the two lattices. 
In the theory of Bragg and Williams the 
energy of the crystal was assumed to be de- 
termined by the long range order, s. As was 





*See R. H. Fowler, Statistical Mechanics, Chap. IV 
(Cambridge University Press, 1936). 
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pointed out by Bethe, this assumption is equiva- 
lent to ignoring the influence of configurational 
fluctuations on the value of the partition func- 
tion. If the interaction of pairs of atoms which 
are not nearest neighbors can be neglected, the 
energy may be expressed as a linear function of a 
parameter, p, which is a measure of the local 
order of the crystal. 


E=E,)—2NakT(i-—p), 
a=2V/2kT, 

Vo=(Virt+ Vo2)/2— Vie, 

where V1, Vee, and Vj. are the mutual energies of 


the indicated types of nearest neighbor pair, 
and ? is given by, 


p=(2/N2) vata, 


a=1 


(3) 


(4) 


where »y, is the number of atoms of type 1 
occupying lattice point a, either zero or unity, 
and n, is the number of sites of lattice 6, neigh- 
boring site a, which are occupied by atoms of 
type 1, a number ranging from zero to z. The 
parameter, ~, may be alternatively expressed in 
the form, 


p=(1/N2)2™ va, (5) 
a,b 

where », is the number of atoms of type 1 occu- 
pying site b, and the sum extends over all pairs of 
a and 6 sites which are neighbors. The constant, 
Ep, is the energy of the two pure crystals, 1 and 2, 
before mixing, provided there is no change in 
the lattice parameters on mixing. We remark that 
Bethe’s local order parameter, o, is equal to 
1—2p. It is immaterial whether we work with 
the variable o or p. The parameter, a, has the 
intuitive advantage of vanishing when the 
average environment of any lattice point is 
randomly distributed between atoms of the two 
types. Specification of a pair of values, (s, p) 
does not of course completely characterize the 
configuration of the crystal. We shall designate 
by w(s, p) the number of different configurations 
corresponding to the pair of values of the .order 
parameters (s, p). 

By statistical mechanics, the free energy, F, of 
the crystal is related to the partition function, 
f(8), in the following manner. 


F=—(1/8) log f(6), (6) 















































where 8 is equal to 1/kT. We shall find it con- 
venient to calculate the configurational partition 
function for a fixed value of the long range order, 
s, and then to minimize the free energy with 
respect to s. The configurational partition func- 
tion for a fixed s may be written as follows. 


1—s 


f(8, s)= dws, or, (7) 


where the sum extends from 0 to 1—s, the maxi- 
mum value of p. The total number of configu- 
rations included by a fixed s for all values of p is 
readily seen to be the product of the number of 
ways the N sites of lattice a may be distributed 
among NV, atoms of type 1 and N2™ atoms of 
type 2, by the number of ways the N sites of 
lattice b may be distributed among NV, atoms 1 
and NV,“ atoms 2. We distribute the sites among 
the atoms rather than the atoms among the 
sites in order to avoid distinguishing between 
configurations which differ only by the inter- 
change of atoms of the same type. We may 
therefore write, 


1l—s N! N! 
eee ee 
p=0 No IN{®! N, INA! 


It is convenient to define a distribution function, 
¢(p), which is normalized to unity. 





Delp) =1; 


9 
No IN, 1N, INQ! | ( ) 
o(p) = — hat. 





Making use of Eqs. (2) and (9), we may write 
(7) in the following form. 





. N! 2 
— 0 


1l—s 


X Le *Nery(p). (10) 


p=0 


The sum in (10) may be expressed in terms of the 
semi-invariants of Thiele,» which are con- 
veniently used to characterize a distribution 
function by its moments. The semi-invariants 


5 See A. Fisher, Mathematical Theory of Probability, Vol. 
I, second edition (Macmillan, 1922). 
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are defined as follows, 


exp > (An/n px"|- > (M,/n Nx, 
E n=0 


‘“ (11) 
1—s 


My=Zp"e(6). 


By logarithmic differentiation the }, may be 
shown to satisfy the following set of linear 
equations, in which the moments M, of the 
distribution function appear as coefficients. 


nsn-—i1 
» ¥ 


m=1\m—1 


\raMen= Mas m=1,2,---. (12) 


Solution of these equations for the first three \, 
yields, 
Ay => M,, 
Ae= M,—M;?, (13) 
A3 = M3- 3M,M24+2M;'. 


By expansion in series and comparison with 
Eq. (11), it is readily verified that the sum in 
(10) may be expressed in the form, 


1— 


= eter g(p) exp ) = (-2va)*/n Da] (14) 


Pp 


Using Eqs. (6), (10), and (14), and approximating 
the factorials of large numbers by Stirling’s 
formula, we obtain the following expression for 
the configurational free energy of the crystal, 


— F(s)/kT=N[2 log 2—(1+s) log (1+s) 


—(1-s) log (1-s)]J+E((—2Na)"/n!)dn (15) 
+2Na. 


The equilibrium value of the long range order 
corresponds to a minimum in the free energy, 
F(s), and is to be determined by the condition ; 


dF(s)/ds=0. (16) 


Before proceeding to the calculation of the 
moments of the distribution function, we remark 
the retention of only the term involving ): 
amounts to assuming ¢(p) to have the form of a 
delta function, 6(p—.,). This leads to the Bragg 
and Williams approximation. Retention of the 
first two terms, involving \,; and de is equivalent 
to assuming that (p—)1)/A2! has a continuous 
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ORDER AND DISORDER 


range of values between — © and + ~, and that 
¢(p) is Gaussian in form. 


o(p) = (1/(ade)#)e~@-A*/2 2, (17) 


Since we shall find by calculation that 1/d2' is 
of the order, N?, the error involved in assuming 
the infinite interval for (b— ,)/A2! is of negligible 
order. The Gaussian form for ¢(p)proves itself to 
be a rather good approximation, and, as we shall 
show, leads to results in substantial agreement 
with those obtained by the more complicated 
method of Bethe. 

For the calculation of the moments of the 
distribution function we make use of the follow- 
ing relations. 


N 
Lme=2N,. 


a=1 


N N 
La=NiM; VnH=Ni; (18) 
a=1 b=1 

A calculation, the details of which are given in 
the appendix, leads to the result, 


Ay = (1—s?)/2, 
he= (1—s?)?/4Nz, 
A3=5?(1 —s*)?/2N22?. 


(19) 


In the expressions for the thermodynamic func- 
tions of the crystal we shall retain only the terms 
in \; and dA». The reader may ascertain with a 
little calculation that the term in \3 makes only 
very small contributions. We obtain for the free 


energy, 


— F(s)/NkT =2 log 2—(1+s) log (1+s) 
—(1-—s) log (1—s)+a(1+s?) 


+o%(1—s%)?/2z++++. (20) 


The equilibrium value of the long range order, s, 
corresponding to a minimum in the free energy 
at constant temperature, (0F/ds).=0, is found 
to satisfy the following trancendental equation, 


s=tanh ys 
y=a—o?(1—s*)/z. 


(21) 


The Bragg and Williams approximation is ob- 
tained if the terms in a/z are neglected in Eqs. 
(20) and (21). We note first of all that Eq. (21) 
has solutions other than s=0, only if @ is 
positive, which we shall see corresponds to a 
negative energy of mixing. If a is negative there 
will be a critical solution temperature and 
separation into two phases. However, treatment 
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of this case requires the discussion of solutions 
of arbitrary composition. When a is positive, 
Eq. (21) has a solution other than s=0, only 
below a certain critical temperature, 7., corre- 

s 


Cc, /k 
1.0 4.0F 











kT/V. —— 
2.37 KT. 2.37 





Fic. 1. Fic. 2. 


sponding to a value of a,, satisfying the following 
equation. 
a?—Za-+z=0, 
a,=(2/2)[1—(1—4/z)!}, 
T.=2V0/2ka.. 


(22) 


For the simple cubic structure, z=6, we find 
T. equal to 2.366 Vo/k, a value which differs 
only slightly from that of Bethe, 2.372 Vo/k. In 
Fig. 1, s is plotted as a function of kT/Vo 
for z=6. The average or thermodynamic 
energy, E=0(F/T)/d(1/T) and the entropy, 
S=(E-—F)/T, may be written in the form, 


E/NkT= —a(1+s?) —a*(1—s*)*/z, 
S/Nk=2 log 2—(1+5) log (1+s) 


—(1-—s) log (1—s)—a?(1—s?)?/2z. (23) 


The configurational heat capacity per atom, 
equal to (1/2N)(dE/dT), has the form, 


i. 21 1—(2a/z)(1-—s?) |? 1 : 
: =f avin ili + 1-si}] (24) 


cosh? ys—y—2a’s?/z 22 





& 


The configurational heat capacity is plotted as a 
function of RT / Vo for z=6 in Fig. 2. It undergoes 
a jump at the critical temperature, c,/k having 
the value, 4.23 just below 7., and the value 0.134 
just above. The peak of c,/k is thus higher than 
Bethe’s value, 2.14. The inclusion of the term in 
hs would lower our value somewhat. The 
configurational heat capacity does not vanish 
above the critical temperature, because a certain 
average degree of local order persists even after 
the long range order vanishes, and its gradual 
destruction with increasing temperature con- 
tributes to the heat capacity. Since s is a con- 
tinuous function of the temperature, E is also a 
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continuous function of the temperature. Thus 
there is no heat of transition, and we have to do 
with a phase change of the second kind at the 
critical temperature. In contrast with a phase 
change of the first kind, there is only one equi- 
librium phase present at the transition tem- 
perature. 

It is possible to apply the present theory to 
face-centered crystals with an atomic ratio, 1 : 3, 
a case considered by Bragg and Williams and by 
Peierls, with satisfactory results. The theory may 
be readily generalized for a crystal constructed 
of two interpenetrating lattices with any compo- 
sition instead of the special 1 : 1 atomic ratio 
considered here. However, when a is positive, 
certain difficulties are encountered at very low 
temperatures. They arise from the failure of the 
series in \, to converge rapidly. If the Gaussian 
approximation, involving only \; and ¢ is used, 
the transcendental equation for s fails to have a 
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nonvanishing solution near the absolute zero, 
except for the atomic ratio, 1 : 1, treated in the 
present article. This is analogous to the difficulty 
encountered in the Heisenberg theory of ferro- 
magnetism, according to which the spontaneous 
magnetization would vanish at low temperatures. 
In our case the difficulty may be removed for a 
region of composition in the neighborhood of the 
atomic ratio, 1 : 1, by including the term in )s. 
Outside this region, terms in the X, of still higher 
order are required. The extension of the theory 
to solid solutions of arbitrary compositions will 
be discussed in detail at a later time. Its limiting 
form for vanishing long range order promises to 
be of use in the theory of liquid solutions. 

In conclusion, the writer wishes to express his 
thanks to Professor H. Bethe of Cornell Uni- 
versity and to Professor F. Bitter of the Massa- 
chusetts Institute of Technology for many 
helpful discussions. 


APPENDIX 


Calculation of the moments 


From Eq. (4) we observe that the moment, M,,, of the 
distribution function, g(p), may be written in the form, 


N 
Mn= (2/Nz)™(( > vats)™)av, (25) 


a=l 


where ay denotes an a priori average over all crystalline 
configurations. Since the a priori distributions on the two 
different lattices a and b are independent, an average of the 
form, ((fi(vi-++vw)fo(mi-+-mNn))a, may be written as 
(fi(vi- ++ v~))av(fo(mi+++2n))ay, inasmuch as the vg refer to 
atoms on lattice a and the m, to atoms on lattice 6. To 
avoid confusion, we remark that this is true since we are 
concerned only with a@ priori averages. It would not be 
true of statistical-mechanical averages, where the Boltz- 
mann factors, of course, introduce correlation between the 
distributions on the separate lattices. We shall employ the 
following sum relations in the subsequent calculations, 


N N 
dra" =rN; Dn. =2wN, 


a=l a=l 
r=N,/N=(1+s)/2, w=N,/N=(1—s)/2. (26) 


These relations follow from Eqs. (2) and (18), and from 
the fact that va«”=va when m is not zero, since the va are 
either zero or unity. 

Since the average values of (va™)a, and (ma) are of 
course the same for all lattice sites, we obtain by averaging 
both sides of Eqs. (26), 


(v™)w=r, (n)w=sw, (27) 


where (v™)a, and (m)ay refer to any lattice site. Similarly, 
(vava')ay has the same value, (vv’)av, for any pair of sites (not 
true of (#a7a’)~y). Thus we have from Eq. (26), 


N 
PN? = ((>ova)*)w=Nv)wt+N(N—1)(vv’)w, (28) 
a=l 
. _r(rN—1) _ _rw 3 ; 
ponatBD an Mo( 1) 


We shall also require the value of (n™)w, the average of 
nq” for any site a. The number of ways in which a set of z 
sites of lattice b may be assigned to m atoms of type 1 
and the remaining N—z sites to N° —n atoms of type 1 and 


N 
N2® —z+n atoms of type2 isevidently equal to ( 
N,O—n 


z 
( ) , while the total number of configurations of lattice b is 
n 


N,o 
P.(n), of configurations in which the z neighbors of a site 
of lattice a include » atoms of type 1, 


N \" N 2 
P.(n) = Ni, Lain n 
5) -»+0( 4). 
-(%), w+0(<) (29) 


The average value (n”),, may then be expressed as follows, 


N 
( ) Thus we may write for the normalized probability, 


(n™) y= > n™P(n) 
n=0 


=[(a"/dx™)(we*+1)* ]rao. (30) 
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The relations stated above suffice for the calculation of Mi 
and M2. For the higher moments others of a similar nature 
would be needed. 
We may calculate M; at once, using Eqs. (25) and (27), 
and obtain the result, 
2 2(v) av) av 
mM=— a) Ay (a) ay = ———- = 
1 wz (v dav a) Av " 
For the moment, M2, we may write, by virtue of the 
independence of averages on the separate lattices, 


(1—s*)/2. (31) 


4 N 
M.= pee (vv" ay) (n? ay (v0 YwX(2 Ma)” wy}. (32) 


From Eq. (30) we find that m is equal to 2’r?+2rw. This 
result together with Eqs. (27), (28), and (32) yields, 
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M.= M+ (4/NZ)rw*+O0(1/N?) (33) 
= M;?+(1—s?)?/4Nz+O(1/N?). 


The third moment, M;, may be obtained in a similar 
fashion, although the calculation is somewhat more 
tedious. When substituted in Eqs. (13), the moments so 
obtained yield Eqs. (19). 

The moments may be calculated in another way if p is 
expressed in the form of Eq. (5). A method analogous to 
that employed by Van Vleck® in his treatment of the 
Heisenberg theory of ferromagnetism may then be em- 
ployed. This method, though somewhat simpler to state, 
involves about the same amount of computational labor 
as the one described above. 


6 J. H. Van Vleck, Electric and Magnetic Susceptibilities, Vol. 12, 
§ 78 (Clarendon Press, Oxford, 1932). 
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The Electron Diffraction Investigation of the Molecular Structures of Ketene and 
Thiophosphoryl! Chloride’ 


J. Y. Beacn? anp D. P. STEVENSON 
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(Received November 23, 1937) 


The structures of the gas molecules of ketene and thiophosphoryl chloride have been de- 
termined by the diffraction of 57 kilovolt electrons. The interpretation of the photographs was 
made by the visual method in conjunction with the radial distribution method. The results for 
ketene are: carbon-carbon distance 1.35+0.02A., and carbon-oxygen distance 1.17+0.02A. 
The carbon-hydrogen distance was assumed to be 1.07A. The results for thiophosphory] chloride 
are: phosphorus-chlorine distance 2.01+0.02A, phosphorus-sulphur distance 1.94+0.03A, 
and the Cl— P—Cl angle 107+3°. The results are discussed in terms of covalent radii and the 
resonance of molecules among several electronic structures. 


HE apparatus used in this investigation is 

closely similar to the one described by 
Brockway’ and used by Pauling, Brockway et al. 
for determining the structures of gas molecules. 
The electrons were accelerated by a potential of 
roughly 57 kilovolts. The wave-length of the 
electrons was 0.0515A. This wave-length was 
found accurately by taking transmission pictures 
of gold foil, the structure of which is known. 
In interpreting the photographs we have used 
the visual method and the radial distribution 
method. Both of these procedures have been 





‘Presented before the American Chemical Society, 
Rochester, September, 1937. 

* National Research Fellow in chemistry. ’ 

* For a general review of the experimental procedure and 
a summary of electron diffraction results see L. O. Brock- 
way, Rev. Mod. Phys. 8, 231 (1936). 





described by Pauling and Brockway.‘ Bauer?® 
has developed an analytic method, based on 
visual measurement, for interpreting photo- 
graphs. The method consists of differentiating 
the theoretical scattering curve and setting the 
derivative equal to zero at each of the measured 
maxima and minima. The interatomic distances 
are then obtained by solving a set of simultaneous 
equations. 

The formula used in the visual method for 
calculating the intensity of the scattered elec- 


4The visual method of interpretation, originally de- 
veloped by R. Wierl, Ann. d. Physik 8, 531 (1931) and 
13, 453 (1932), has been tested and discussed by L. Pauling 
and L. O. Brockway, J. Chem. Phys. 2, 867 (1934). The 
radial distribution method was developed, tested and 
applied to a large number of compounds by L. Pauling 
and L. O. Brockway, J. Am. Chem. Soc. 57, 2684 (1935). 

5S. H. Bauer, J. Chem. Phys. 4, 406 (1936). 
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Fic. 1. Radial distribution curves for ketene and 
thiophosphory] chloride. 


trons as a function of the scattering angle is 
I= ,2;Z,Z;(sin $15;/SPij), (1) 


where Z; and Z; are the atomic numbers of the 
ith and jth atoms, 7;; is the distance between 
the ith and jth atoms, and s is 


(4m sin 0/2)/X. 


Both summations are taken over every atom in 
the molecule. \ is the wave-length of the electrons 
and @ is the angle through which the electrons 
are scattered. To evaluate Eq. (1) a model for 
the molecule is assumed, the 7;;’s are calculated 
from the geometry of the model and put into 
the formula. The intensity, in arbitrary units, is 
calculated and plotted in terms of s. This pro- 
cedure is repeated for a number of molecular 
models, the number depending on the com- 
plexity of the molecule and on the pre-existent 
ideas regarding its structure. The shape of the 
intensity curve, that is the relative heights and 
spacings of the maxima and minima, is de- 
pendent only on the shape assumed for the 
molecule and not on the assumed size. The size 
assumed affects only a scale factor which is 
automatically corrected when the ring diameters 
of the maxima and minima on the films are com- 
pared quantitatively with the positions of the 
maxima and minima on the theoretical curve. 
In general many of the models for the molecule 
can be eliminated from a qualitative comparison 
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of the photographs with the calculated curves. 
By qualitative comparison is meant a com- 
parison of the relative heights and spacings of 
maxima and minima between photographs and 
curves without resorting to measurement of the 
photographs. The final results are obtained from 
a quantitative comparison of the photographs 
with the qualitatively acceptible calculated 
curves. 

The radial distribution method consists of 
inverting the intensity formula. Instead of calcu- 
lating the intensity of scattering as a function of 
the structure, the density of atoms in the mole- 
cule is calculated in terms of the observed 
scattering. This distribution function is obtained 
directly from measurements on the photographs, 
no assumption of molecular model being in- 
volved. These functions for ketene and _thio- 
phosphory! chloride are shown in Fig. 1. 


KETENE 


The ketene used in this investigation was pre- 
pared by Dr. W. E. Vaughn, Dr. Elgene A. 
Smith and Mr. M. A. Dolliver of Harvard 
University. From the time the preparation was 
made until the photographs were taken it was 
kept in an evacuated vessel at the temperature 
of solid carbon dioxide and alcohol. During all 
of the time the ketene was kept cold it remained 
a clear colorless liquid. Just before it was dis- 
carded the dimer began to deposit as a colorless 
solid on the bottom of the container. The 
polymer remained colorless until allowed to 
warm up to room temperature in the presence 
of air, at which time it turned brown. The 
presence of some of the dimer at the time the last 
pictures were taken could not influence the pho- 
tographs because of the low vapor pressure of 
the dimer at low temperatures. The boiling 
point of the dimer is about 128°C. The pictures 
were taken with the sample at —75°C. 

Visual inspection of the photographs shows 
five maxima. In addition, on some of the photo- 
graphs, there is a spurious maximum close to the 
intense central image. When this occurs it is 
easy to detect as its position is well inside the 
possible position of diffraction maxima and is 
quite variable. The first maximum is very 
intense. The second appears to be a small 





ne w re | 


Res 


T 
sho 
and 
carl 
tan 
long 

T 
it is 
moc 
posi 
inte 
hyd: 
med 
in a 
dist 
the | 
tetré 

aton 
that 
cann 
dista 
dista 
scatt 
deter 
in th 
tance 





eee COC 





ELECTRON DIFFRACTION 


maximum or shelf on the outside of the first. 
The third maximum is sharp. The fourth is less 
sharp and is close to the third. The fifth is sharp 
but is very faint and is farther from the fourth 
than the fourth is from the third. The visually 
estimated intensities along with the measured 
values of s are given in the third and fourth 
columns on Table I. 


TABLE I. Ketene. 











MAX. MIN. I Sobs Scale C-CA 

1 20 5.657 5.93 1.442 
2 (shelf) 7 8.193 

3 9.578 8.96 1.291 

3 10 10.53 10.43 1.367 

4 12.04 11.42 1.308 

4 4 12.91 12.30 1.315 

5 1 15.45 15.35 1.371 

Ave. 1.349 























Results: C—C 1.35+0.02A 
C—O 1.17+0.02A 
C—H 1.07A (assumed) H—C—H angle 109° 28’ 
(assumed) 








The five term radial distribution function is 
shown in Fig. 1. It possesses maxima at 1.37A 
and at 2.51A. The first peak represents the 
carbon-carbon and the short carbon-oxygen dis- 
tances unresolved. The second peak gives the 
long carbon-oxygen distance. 

To calculate the theoretical scattering curves 
it is necessary to evaluate Eq. (1) for several 
models. No attempt is made to determine the 
position of the hydrogen atoms. In calculating 
intensities it was assumed that the carbon- 
hydrogen distance is 1.07A, which is inter- 
mediate between the carbon-hydrogen distance 
in acetylene (1.06A) and the carbon-hydrogen 
distance in methane (1.09A). The angle between 
the carbon-hydrogen bonds was assumed to be 
tetrahedral. The scattering from the hydrogen 
atoms is such a small part of the total scattering 
that the assumptions regarding their position 
cannot affect the results obtained. There are two 
distances to be determined: the carbon-carbon 
distance and the carbon-oxygen distance. The 
Scattering curve, except for a scale factor, is 
determined entirely by the shape of the molecule, 
in this case the ratio of the carbon-carbon dis- 
tance to the carbon-oxygen distance. Accord- 





curve F, 1.38/1.13. 
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ingly several intensity curves are calculated for 
ratios of the carbon-carbon distance to the 
carbon-oxygen distance varying from 1.38/1.13 
to 1.38/1.36. These curves are shown in Fig. 2. 
From a qualitative comparison of the curves 
with photographs it is found that the most 
satisfactory curve is the one for which the ratio 
is 1.38/1.20. The fifth maximum, counting the 
shelf on the first maximum as a separate maxi- 
mum, on the 1.38/1.16 curve is too low, while 
the second maximum (or shelf) on the 1.38/1.24 
curve is too low. Quantitative comparison of the 
photographs with the 1.38/1.20 curve is given in 
Table I. The final values are: for the carbon- 
carbon distance, 1.35+-0.02A and for the carbon- 
oxygen distance 1.17+0.02A. 
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Fic. 2. Theoretical intensity curves for ketene. The 
ratios of the carbon-carbon distance to the carbon-oxygen 
distance are: curve A, 138/132; curve B, 138/1.28; curve 
C, 1.38/1.24; curve D, 1.38/1.20; curve E, 1.38/1.16; 






THIOPHOSPHORYL CHLORIDE 


The thiophosphoryl chloride used in this in- 
vestigation was an Eastman product purified by 
repeated fractional distillation in vacuo by Mr. 
F. B. Jennings. The photographs, taken with 


the sample at 50°C, show six rings. The first - 


three are sharp, the second being the most 
intense. The fourth maximum is less intense than 
the third and is close to the third. The fifth is 
intense and sharp. The sixth is weak and is 
difficult to measure. The visually estimated 
intensities are given in Table II. 
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Fic. 3. Theoretical intensity curves for thiophosphoryl 
chloride. Denoting the Cl— P—Cl angle by @ and the ratio 
of the phosphorus-chlorine distance to the phosphorus- 
sulfur distance by p, the curves are for the following 
models: curve A, a=100° and p=2.02/1.93; curve B, 
a@=105° and p=2.02/1.93; curve C, a=1f09° 28’ and 
p=2.02/1.93; curve D, a=115° and p=202/1.93; curve E, 
a=100° and p=2.02/1.85; curve F, a=105° and p=2.02/ 
1.85; curve G, a=109° 28’ and p=2.02/1.85; curve H, 
a@=115° and p=202/1.85; curve J, a=105° and p=2.02/ 
1.89. 
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TABLE II. Thiophosphoryl chloride. 














Max. | Min. I ww Sette P-CIA 
"3 1.691 1.59 1.899 
1 14 2.442 2.43 2.010 
2 3.388 3.17 1.890 
2 25 4.238 4.18 1.992 
3 5.449 5.33 1.976 
3 17 6.290 6.42 2.062 
4 7.326 7.29 2.010 
4 7 8.036 7.98 2.006 
5 9.132 8.98 1.986 
5 12 10.12 10.11 2.017 
6 11.83 
6 2 13.19 
Ave. 2.007 























(1st max. and 1st and 2nd min. omitted) 
Results: P—Cl 2.01+0.02A 


P—S_ 1.94+0.03A 
Cl—P—Cl angle 107°+3°. 








The five term radial distribution function is 
shown in Fig. 1. In calculating the radial dis- 
tribution function the sixth maximum was 
omitted because of the difficulty in measuring it. 
The peak at 2.02A represents the phosphorus- 
chlorine and the phosphorus-sulfur distances 
unresolved. The peak at 3.27A represents the 
chlorine-chlorine and the sulfur-chlorine dis- 
tances unresolved. The small peak at 1.38A is 
without significance. 

in calculating theoretical intensity curves the 
only models considered were those having the 
symmetry C;,. Having made this assumption— 
that the three chlorine atoms are equivalent— 
two parameters are necessary to specify the shape 
of the molecule. These parameters are taken to 
be the angle between two phosphorus-chlorine 
bonds, a, and the ratio of the phosphorus- 
chlorine distance to the phosphorus-sulfur dis- 
tance, p. Theoretical intensity curves were 
calculated for models having values of a ranging 
from 100° to 115° and p ranging from 2.02/2.10 
to 2.02/1.77. Sixteen curves were calculated. 
Nine of these are shown in Fig. 3. The top four 
show the effect of changing a from 100° to 115° 
keeping p constant at 2.02/1.93. The next four 
show the effect of varying a with p constant at 
2.02/1.85. Curve I is one of the qualitatively 
acceptable curves. From the curves shown and 
other curves it is seen that values of a of 100° or 
less and 115° or greater are unsatisfactory be- 
cause of the prominence of the fourth maximum 
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ELECTRON DIFFRACTION 


on the photographs. The 105° models are unsatis- 
factory when p is greater than 2.02/2.02 for the 
same reason. The tetrahedral models are unsatis- 
factory when p is greater than 2.02/2.10 because 
the average of the chlorine-chlorine and sulfur- 
chlorine distances is greater than the 3.27A 
required by the radial distribution curve. The 
tetrahedral and 105° curves are unsatisfactory 
when p is less than 2.02/1.77 because the average 
of the chlorine-chlorine and sulfur-chlorine dis- 
tances is less than the 3.27A required by the 
radial distribution curve. All the curves within 
these limits are qualitatively acceptable. As the 
most probably correct model we have chosen the 
one for which a is 107° and p is 2.02/1.95. This 
curve is not shown in Fig. 3 but is practically 
identical in appearance with curve I. Quantita- 
tive comparison of the photographs with the 
107°, 2.02/1.95 curve is shown in Table II. 
The final values are: for the phosphorus-chlorine 
distance 2.01+0.02A, for the phosphorus-sulfur 
distance 1.94+0.03A and for the Cl—P—Cl 
angie 107°+3°. 


DISCUSSION 


The distances found in ketene are not those 
that would be expected if ketene has the elec- 
tronic structure 


The expected distances for this structure, from 
Pauling and Brockway’s® recent revision of the 
covalent radius table of Pauling and Huggins,’ 
are 1.34A for the carbon-carbon bond and 1.24A 
for the carbon-oxygen bond. It is seen that the 
carbon-carbon distance is greater and the carbon- 
oxygen distance is less than expected. This is 
immediately accounted for by resonance of the 
molecule between the above structure and the 
structure 


H 
_ + 
C—C=0:. 


/ 
H 


°L. Pauling and L. O. Brockway, J. Am. Chem. Soc. 59, 
1223 (1937). 
(s8 ; Pauling and M. L. Huggins, Zeits. f. Krist. A87, 205 
4). 
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It is possible to estimate the extent of the reso- 
nance by using the relation between bond type 
and bond distance given by Pauling, Brockway 
and Beach’ and revised by Pauling and Brock- 
way.’ This procedure leads to the result that 
the structure is three quarters double-bond, 
double-bond and one quarter single-bond, triple- 
bond. 

The P—Cl distance in PSCI; (2.01A) is about 
the same as the P—Cl distance in PCl; (2.00A) 
obtained by Brockway and Wall® and Pauling 
and Brockway.'® The sum of the single covalent 
radii for phosphorus and chlorine is 2.09A. The 
shortening of the bond is ascribed by Brockway 
and Jenkins" to partial double bond character 
due to resonance of the molecule to structures of 


the type 
P 
@ \\ 


+Cl Cl. 
Cl 





Structures of the type 
P 


4 
+Cl Cl— 
Cl 


may also contribute, as suggested by Gregg, 
Hampson, Jenkins, Jones and Sutton.” The pos- 
sibility that the shortening is due to ionic 
character of the bonds was eliminated by 
Brockway and Jenkins by a study of several 
metal methyls. In these compounds, the atoms 
of which have greater relative electronegativities 
than phosphorus and chlorine, the distances are 
normal. Here, double bond formation cannot 
occur because of the necessity for observing the 
octet rule on carbon. It was therefore concluded 
that the shortening in PCl; (and other similar 
compounds) is due to partial double bond char- 
acter and not to ionic character. In our discussion 
of PSCl; we shall accordingly assume that ionic 
character cannot affect the length of the bonds. 

8L. Pauling, L. O. Brockway and J. Y. Beach, J. Am. 
Chem. Soc. 56, 2706 (1935). 

9L. Brockway, and F. T. Wall, J. Am. Chem. Soc. 56, 
2373 (1934). 

10L, Pauling and L. O. Brockway, J. Am. Chem. Soc. 
57, 2684 (1935). 

1L. O. Brockway and H. O. Jenkins, J. Am. Chem. Soc. 
58, 2036 (1936). 


12 A. H. Gregg, G. C. Hampson, G. I. Jenkins, P. F. L. 
Jones and L. E. Sutton, Trans. Faraday Soc. 33, 852 (1937). 
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To a first approximation the structure of 
PSC]; is 


Cl 

—- + 

$C. (1) 
Cl 


The contributing structures which will tend to 
shorten the P—Cl bond are of the following type: 


Cit 
+O 
-S—P—Cl (II) 
Cl_ 
Cit 
VA 
-S—P—Cl . (III) 
Cl 


Structure (II) will not contribute, however, 
because of the adjacent charge rule.’** Probably 
a more important (lower energy) structure, 
tending to shorten the P—S bond, is 


Cl 
/ 
S=P—Cl. (IV) 
Cl 


The structure 
Cl+ 


VA 
S=P—Cl (V) 
cK. 


may also contribute to some extent. 

In these last three structures the octet rule is 
not observed for phosphorus. However, violation 
of the octet rule in atoms in other than the first 
row of the periodic table is permissible; for 
example PF;, etc. The structures that shorten 
the P—Cl bond in PSCI; are less numerous and 
probably less important than the ones for PCls. 
It is therefore surprising that the P—Cl distance 
in PSCl; is not substantially greater than in 
PCl;. The answer to this question is found in the 
discussion of the short N —O distance in nitrate 


2 [,, Pauling and L. O. Brockway, J. Am. Chem. Soc. 
59, 13 (1937). 
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ion given by Elliott.1* The short distance is 
attributed to the resultant positive charge on 
the nitrogen atom. In the most important 
structure for PSCl; there is a resultant charge of 
plus one on the phosphorus atom which by 
analogy with the small radius for nitrogen in 
nitrate ion leads one to expect a short P—Cl 
distance in PSCl;. This effect apparently just 
offsets the lack of suitable contributing double 
P—Cl bond structures mentioned above, leaving 
the P—ClI distance in PSCI; the same as in PCls. 

If the electronic structure of PSCl; is simply 
structure (I) written above, the P—S bond being 
a single bond between oppositely charged atoms ; 
the considerations of Elliott would predict that 
the bond distance is the single bond distance. 
The increase in the single bond radius of the 
negative sulfur atom would be offset by the 
decrease in the radius of the positive phosphorus 
atom. That this prediction leads to correct 
results is seen in the case of a large number of 
tetrahedral crystals. In these crystals the most 
important electronic structure, responsible for 
the tetrahedral orientation of bonds, is one in 
which bonded atoms have equal and opposite 
formal charges. In these cases the interatomic 
distances are observed to be just the interatomic 
distances for normal valence bonds for which 
there are no resultant charges on bonded atoms. 
The observed P—S distance in PSCI; is 1.94A, 
nearly the same as the double bond distance for 
this bond (1.95A). Having ruled out any effect 
of ionic character and formal charges on the 
length of this bond, the shortening below the 
single bond distance must be due to a consider- 
able (about one-half) contribution of structure 
(III). It is difficult to see how the P—S distance 
can be less than the double bond distance. This 
effect is also observed in the case of sulfate and 
phosphate crystals. In the cases investigated the 
distance between the central atom of the anion 
and the oxygen atom is found to be roughly the 
double bond distance. 

We wish to express our appreciation to Pro- 
fessor C. P. Smyth for suggesting that the appa- 
ratus be built and for his continued interest in 
the research. 


18 Norman Elliott, J. Am. Chem. Soc. 59, 1380 (1937). 
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Titanium salts should behave more simply than most 
materials under magnetic cooling, since the ground level 
of the paramagnetic ion is here inherently doubly de- 
generate and so unaffected by crystalline fields. De Haas 
and Wiersma show that as long as the magnetic moment 
remains constant during adiabatic demagnetization, the 
temperature is a linear function of a field strength, a most 
valuable fact since it permits the establishment of a thermo- 
dynamic temperature scale. The present paper computes 
when deviation of the moment from constancy should 
commence because of the perturbing effect of second-order 


terms (the first nonvanishing order) in the dipole-dipole 
coupling energy. The calculated moment does not remain 
constant to as low field strengths as reported experimen- 
tally. The discrepancy cannot be ascribed to a counter- 
balancing effect of exchange coupling but may be due to 
the neglect of third and higher order dipole terms. Because 
of these large dipole perturbations, the status of the experi- 
mental temperature scale in titanium alum becomes rather 
uncertain. The present stage of the theory may be tested 
more accurately by using much lower initial fields than in 
the existing Leiden experiments. 





INTRODUCTION 


HE production of very low temperatures by 

the method of adiabatic demagnetization is 
a subject already classic. However, the great 
difficulty is the determination of the thermo- 
dynamic temperature. This is usually estimated 
by assuming Curie’s law. The latter, however, 
cannot be accurate, as it would require zero 
specific heat in the absence of an applied mag- 
netic field, and so allows no heat capacity to 
prevent the substance from cooling clear down to 
the absolute zero! 

A neat way out of the difficulty, though of 
somewhat limited applicability, has been pro- 
posed by de Haas and Wiersma,' who suggest and 
try the following procedure. As the field is 
diminished adiabatically from the high value 
which produced the initial decrease of entropy at 
helium temperatures, there is a certain interval 
during which the moment retains the same value 
as at the original temperature and field strength. 
It is purely a matter of thermodynamics to 
show, as they do, that as long as this is the case, 
the thermodynamic temperature is a linear 
function of field strength, so that 


T=aH+H0. (1) 


The constants a and 6b in general involve the 
intensity of magnetization M as a parameter, 
and can be determined by finding two known 





* At present, visiting professor, Princeton University. 
(1936) J. de Haas and E. C. Wiersma, Physica 3, 491 


pairs of values of H, T in the helium range 
which yield the same given value of M. Once a 
and b have been fixed, (1) can be used to yield T 
in terms of H at low temperatures where T 
cannot be measured by ordinary means. If the 
field is reduced below a certain critical point, 
however, the moment deviates from its original 
value, in general decreasing. Below this point, 
Eq. (1) ceases to apply. 

To illustrate, de Haas and Wiersma find that 
caesium titanium alum has the same moment 
when demagnetized to 105 gauss as under the 
initial conditions of 1.20°K and 24,075 gauss. 
They thus conclude that the temperature as- 
sociated with 105 gauss is 1.20 (105/24,075) or 
0.0052°K. Below 105 gauss the moment decreases 
rapidly, and at 1.05 gauss has less than 1/50 of 
its initial value. 

In selecting a material to be measured, it is 
important to find one which comes as close as 
possible to an “‘ideal’’ behavior, i.e., to zero 
specific heat when H=0. When H0 there is a 
finite, socalled magnetic specific heat even for an 
ideal substance. It is when H is reduced so far 
that the heat capacity due to other causes 
begins to be comparable with this magnetic 
specific heat that (1) starts to lose validity. 

De Haas and Wiersma’s choice of caesium 
titanium alum was an exceedingly wise one. In 
the first place it is magnetically dilute. Secondly, 
it has a spin of S=} and susceptibility measure- 
ments show that it is in a nondegenerate orbital 
state, at least below hydrogen temperatures or 
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so. On superposing both spin and orbit, one thus 
sees that the Ti+++ ion is in a doubly degenerate 
state. Now a fundamental theorem of Kramers? 
tells us that no crystalline field whatsoever is 
capable of separating the levels in a doubly 
degenerate state with an odd number of electrons. 
Hence there is no Stark splitting to create 
deviations from Curie’s law. The only apparent 
cause of deviation from an ideal behavior is 
therefore the magnetic dipole-dipole coupling, 
and possibly a little exchange interaction, be- 
tween the widely spaced paramagnetic ions. In 
the present paper we calculate when these effects 
should begin to enter. 

The.computation of the influence of the dipole- 
dipole coupling involves no unknown constants, 
as the lattice spacing is known. We find that 
when terms of the second order (the first non- 
vanishing order) in this coupling are included, the 
moment at 105 gauss should be about 18 percent 
lower than the value associated with the initial 
field 24,000 gauss. Experimentally, the moments 
are the same at these two field strengths to 
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within the experimental error, estimated by de 
Haas and Wiersma as five percent. Conceivably 
the discrepancy arises from the neglect of third 
and higher order dipole terms, but in any event 
the existence of such large dipole perturbations 
casts serious doubts on the propriety of using (1), 
without further investigation, to determine tem- 
peratures as low as 0.0052°K, the value yielded 
by (1) at 105 gauss. Quite irrespective of quanti- 
tative estimates, simple qualitative consider- 
ations show that one should expect ideal con- 
ditions probably not to persist down to 0.0052°K. 
Namely, the dipole-dipole perturbations should 
be appreciable when the temperature ceases to 
be high compared with the characteristic temper- 
ature* +=3N?/k associated with dipole-dipole 
interaction. Here N is the number of atoms/cc, 
and 6 is the Bohr magneton he/4rmc. Now 
caesium titanium alum has r=0.0038°K and so 
one would a priori expect the dipole corrections 
to be appreciable at temperatures considerably 
above 0.0052°K. 


MATHEMATICAL DETAILS 


It is necessary to have formulas for the entropy and magnetic moment of an ensemble of coupled 





atoms. Practically all the requisite calculations have been made in a previous paper by the writer,‘ 
henceforth called |.c., in which the partition function was computed in Eqs. (21), (23), (27) to the 
second approximation in the coupling parameters, with arbitrarily large external applied fields. As 
the demagnetization measurements are made on a powder, the anisotropic term proportional to 
Q, in Eq. (21) of l.c. can be omitted, since it averages to zero. The effective field H is now to be taken 
as identical with the applied field, rather than as primarily an inner Weiss field such as was used in the 
previous application of the calculations to ferromagnetism. 

Once the partition function Z has thus been determined, it is a matter of simple differentiation to 
find the entropy S and magnetic moment WM, viz., 


=-—0F/dT, M=-—0F/dH, F=—kT log Z. 
Thus it is found that 


S/kN =log 2+log cosh (g8H/2kT) —(g8H/2kT) tanh (g8H/2kT) 
+(r/ T)fs'+(7/T)(fs"+Fst+Gs), (2) 


M/gBN =B,+(7/T) fu’ +(7/T)*(fu"’ + Fut+Gm), (3) 


where B,, fs:--Gy are functions only of the argument x=g8H/2kT and are defined as follows: 


2H. A. Kramers, Proc. Amsterdam Acad. 33, 959 (1930). 

8Cf. J. H. Van Vleck, J. Chem. Phys. 5, 320 (1937); M. H. Hebb and E. M. Purcell, ibid., p. 338. ; 

4J. H. Van Vleck, Phys. Rev. 52, 1189 (1937). The following errata may be noted in this paper. In the first line ot 
Eq. (21), read 4 in place of y. In the third relation of Eq. (27), Q should read Qo and in the formula for 2’v on the first 
line of p. 1193, N should read N-. Because we are usually dealing with measuring fields large enough to produce saturation, 
it is necessary to use the calculations of this article rather than the treatment appropriate to weak final fields given in 
reference 3. Hebb and Purcell (reference 3) have already discussed calorimetric measurements on the field-free specific 
heat of caesium titanium alum. 
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ADIABATIC DEMAGNETIZATION 


B,=} tanh (g8H/2kT), fs’ =}(@—v2’)(—gBH/kT)(1—4B,")B,, 
Of 5" = (& —v2')?[8B,4— 2B, — (g8H/kT)(1—4B,")(Bi—8B;') ], 
9F s= —(2v)[4B,4+3 —4B,?+ (g8H/kT) (4B —2B,)(1—4B,’) ], 
9Gs= —QoLt —B?J[3—16B?+ }(g8H/kT)(10B, —32B,*) |. 
fu’ = 3 (—2'v)(1—-4B) By, Of" =(@—2'v)?(1—4B,”)(56,—8B;'), 9Fy=2''v?(1—4B,*)(4B,3—2B)), 
9G u = QoL} (RT /g8H)*(24B,8 — 12B,) }4+Qo(4—16B,?)[ — §Bi+3Bi3+ 3(kT/g8H) (i —$B,’) ]. 


Here g=2, Qo=144, 2s=12, 


7 =3Np?/k =0.0038°K, (4) 


B=he/4rmc, v= —2vV2J/g¢°8°N, 2'=6v2, 2'’=r*N*z=6, (5) 


where J is the usual exchange integral between 
adjacent atoms. Exchange coupling between the 
atoms is included as a possibility in order to 
examine whether this can sometimes fortuitously 
neutralize the perturbing effect of the dipole- 
dipole forces. In (4), (5) we give the values of 
various expressions appropriate to caesium 
titanium alum, i.e., a face-centered lattice with 
practically all the magnetism arising from spin. 
For increased generality in connection with 
possible future uses, we mention that (2), (3) 
still apply to other types of cubic lattices 
provided the definitions (4), (5) of the symbols 
are replaced by more general ones given in l.c.*: 4 

The G terms in (2) and (3) are the pertur- 
bations due to dipole-dipole interaction, which 
enter only in the second order (r?). The terms f’ 
and f’’+F give respectively the first and second- 
order effect of the exchange (v) as well as 
demagnetization (®) corrections. The significance 
of the demagnetization plus Lorentz factor ® is 
explained in l.c. and earlier papers,*: * and need 
not be elaborated here. Anyway, ® vanishes for a 
sphere, the shape of the test body used in the 
Leiden experiments. 

If we omit all the correction terms of the 
various types f, F, G, Eqs.- (2), (3) reduce to 


M/g8N =} tanh x, (6) 
S/kN=log 2+log cosh x—xtanhx, (7) 


where x=g8H/2kT. Then (1) clearly applies 
with b=0. We show below in fine print that 
(6), (7) are also valid when the f terms are 
retained, provided the argument is taken as 


x= g6[ H+ (b—2'v)M |/2kT. (8) 
In this case (1) is still correct, with now 


b/a=(&—2'v) M. (9) 











One way to prove (6), (7), (8) is to expand (6), (7) as 
Taylor’s series in (#—z’v), using successively more ac- 
curate values of V on the right side. Then (2), (3) and (6), 
(7) agree to terms in (@—v’z)? inclusive, provided F=G=0. 
Another way of establishing the result is to utilize the well- 
known fact that the demagnetization plus Lorentz or ® 
corrections are equivalent to use of an effective local field 
H+#4M in the formula for the moment. If the v as well as ® 
terms are included, it is clear that the effective field is 
H+(&—2'v) M, since and —2’z enter in a simple additive 
fashion in f. One immediately suspects that as far as the 
entropy is concerned, the f terms are also incorporated by 
employing a field of this character. An analytical proof that 
(6, 8) imply (7), except for an easily determined function 
of T alone, is furnished by the thermodynamic relation 


(0M/dT)y=(0S/0H)r 
inasmuch as by (6) and (8) 


(ax/8T)y/(dx/8H) r= —(8H/aT). 
= —(8H/dT) yu = —2kx/gp. 


In view of the foregoing, we can drop the f 
terms in (2), (3), as they do not affect the 
constancy of the moment, or the validity of (1). 
We henceforth interpret 7 not as the applied field 
H, but rather as the local field H,+(@—2'v)M. 
Our proof of the equivalence of the f terms to this 
effective field ceases to apply with complete 
rigor when the F and G corrections are super- 
posed, but the resulting error is unimportant for 
our purposes, as to a sufficient approximation the 
various contributions may be considered additive. 

Since S is constant under adiabatic conditions, 
Eq. (2) determines T as a function of H. The 
procedure is to solve (2) for T in terms of H, or 
more simply, for 7 in terms of x=g8H/2kT. 
The resulting pairs of values of 7, H (or T, x) are 
then substituted in (3) and in this way the course 
of the moment during the adiabatic demag- 
netization can be calculated. 

In computing the influence oi the F and G terms, it is 


often possible to make one of two simplifications appro- 
priate to different limiting conditions. 
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TABLE I. 








Ho=24,000 gauss, Tp>=1.17°K (xo=1.37) 





H 430 250 105 1.05 
(T/T) cale 5.50 3.15 1.20 1.61 
(x0/X)eale 993 972 .29 | 120 
(M/Mo)ecaie 99 97 82 001 
(M/Mo)obs 1.00 016 








H,.=12,000 gauss, To=1.25° (xo =0.64) 





H 200 115 64 5.23 
(T/t)eale 5.7 3.60 2.81 2.68 
(x0/X)cale 1.04 1.14 1.60 18.8 
(M/Mo)ecale .93 79 BS | .040 
(M/Mo)obs .050 














H,=7000 gauss, To>=1.17° (xo=0.4) 





ni 145 95 50 5.23 
(T/r)ecale 7.20 a5 4.4 4.02 
(x0/X)ecale 1.14 1.33 2.00 17.4 
(M/Mo)eale 85 BY | A7 .052 











Down to fields of about 200 gauss one can usually with 
fair approximation assume that the argument of the F and 
G terms, which are perturbations, has the same value xo 
as at the initial helium temperature, while the leading 
terms in S and M, i.e. those of structure (6), (7) can be 
developed as a Taylor’s series in x—%Xo, and only terms 
through the first order retained. On remembering that S 
remains constant, it is thus found that 


x—x0=[(Fs°+Gs°)/x0(1—-4Bi") ](r/T0)*(Ho/H)? (10) 


M—M°=gBN[}(Fs°+Gs°) xo 
+ Fu°+Gu°](1/T0)*(Ho/H)? (11) 

where Ho, T are the original helium values of H, T, and 
where F° means that F is to be evaluated with the ar.su- 
ment x=Xo, etc. 

At very low field strengths (x about 0.3 or lower), a 
Taylor's development in ascending powers of H can be 
made, v7z. 


S/RN =log 2—4x2—}(1/2T)?(Qo+2'”) 
+3(1/2T)*x*(3Q0+42"'v"), (12) 
M/gBN = }x—(1/3T)?x(3Q0+2""0*) 
— 4x34 090 9(9772/675 72) +.282/"7(1172/5472). (13) 


Eq. (12) is readily solved for (7/T)?, and the result sub- 
stituted in (13). 


EFFECT OF EXCHANGE COUPLING 


The F terms in (10), (11), or v terms in 
(12), (13) represent the part of the exchange 
corrections which cannot be incorporated by the 
artifice of a local field proportional to M. The 
expression }Fs°x~!+ Fy in (11) is easily seen to 
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be always negative. Hence exchange effects tend 
to decrease the moment in adiabatic demag- 
netization, regardless of the sign of 9, i.e., of the 
exchange integral. The constancy of the moment 
reported in the Leiden experiments thus cannot 
be attributed to a fortuitous cancelation of 
exchange and dipole-dipole corrections, a hy- 
pothesis which looked hopeful before the requisite 
mathematical calculations were made. We hence- 
forth omit the exchange terms F, as they do not 
help matters any. It is, of course, reasonable that 
the exchange integrals should be very small, 
because of the large spacing between para- 
magnetic ions. The most important effect of any 
possible exchange coupling is to give rise to a 
nonvanishing constant b in (1), since Dd is by (9) 
of the first rather than second order in the 
exchange integral. Such a } term, of course, does 
not affect the constancy of the moment, but 
might prove troublesome in calculating the 
temperature, inasmuch as 0 is negligible com- 
pared with aH at helium temperatures and so 
cannot be detected by calibration experiments 
made there, but may be appreciable at the final 
low temperatures, where aH is small. 


EFFECT OF DIpoLE-DIPOLE COUPLING 


Table I illustrates the effect of dipole-dipole 
forces as calculated from Eqs. (2), (3). The 
values of H given in the table often are not round 
numbers because x0/xo rather than H is con- 
veniently taken as the independent variable in 
numerical solution. The ratio x0/x=(THo/HT») 
is the same as the factor by which the tempera- 
ture deviates from the “‘ideal’’ value (1). 

One immediately speculates as to the cause of 
the rather glaring discrepancy between the 
observed and calculated behavior in the case of 
the demagnetization from 24,000 gauss. As far 
as we can see, the only possibility is to attribute 
it to neglect of third and higher powers of the 
dipole-dipole interaction, as Eqs. (2), (3) are 
correct only to the second power, the first 
nonvanishing order, in this coupling. The effect 
of each successive power is dimensionally of the 
order +/T compared with the preceding. Un- 
fortunately the ratio r/T is not very small, and 
so the convergence of the development of the 
partition function is notoriously bad in problems 
of the present sort. 
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The ultimate inadequacy of Eqs. (2), (3) can be seen in 
a particularly striking way by the graphical construction 
given in Fig. 1. If we plot S against 1/H?, Eq. (2) shows 
that the curves of constant x=g8H/2kT are straight 
lines of slope 7(k?x?/8?)Gs, as illustrated in Fig. 1. Since the 
entropy retains a constant value So during demagnetiza- 
tion, the magnitude of x and hence of the temperature 
appropriate to the final field strength H can be found by 
locating the point in Fig. 1 whose ordinate is So and whose 
abscissa is 1/H?, and then finding which line passes 
through this point. Fig. 1 shows that often the lines cross, 
so that many points yield more than one value of x. 
Such a situation is impossible, as it would require that the 
temperature be a two-valued function of the field strength, 
an absurdity. Also, if the curves crossed, to the right of 
the crossing we would have higher S correlated with higher 
x, and hence (0S/0T)x <0, i.e., a negative specific heat, 
something obviously impossible physically, and also im- 
possible mathematically because of the exponential 
character of the partition function. Hence the lines of 
large x which start with positive slope doubtless ultimately 
tend to curve downwards when the neglected higher powers 
of the dipole-dipole interaction are included, as we have 
indicated schematically by the dashed lines in Fig. 1, 
whose exact position is, of course, only conjectural. The 
lines corresponding to small x are reliable further to the 
right than those of large x because the former correspond 
to smaller values of +/T for given H. 

To the extreme left of Fig. 1, the effect of the neglected 
higher members is unimportant and the lines are reliable. 
Consequently Fig. 1 shows that for xo greater than about 
0.9, the ratio of temperature to field strength actually 
decreases for a while during adiabatic demagnetization. 
(This is reflected by some values of xo/x being less than 
unity in the table for x=1.37.) In the vicinity of x=.9, 
there is a line of zero initial slope, and for this particular 
field strength, Eq. (1) is valid even when second order 
terms in the dipole-dipole interaction are included. 

The table for x = 1.37 shows actually a lower temperature 
for H= 105 gauss than for H = 1.05. This absurd situation is 
because x has been determined from a rising line for 
H=105, the only choice available here (as So/kN=0.22), 
while a descending line has been employed for H=1.05. 
The values of +/T are here so high that obviously no 
quantitative significance can be attached to our calcula- 
tions with straight lines. 


Of course the mere fact that Fig. 1, or general 
considerations, show that third and higher order 
dipole terms cannot really be neglected, by no 
means makes it clear that the neglected members 
will just cancel the included second-order effects, 
so that the moment has the same value at 105 
and 24,000 gauss. Indeed, the substitution of the 
dashed curves for the solid lines in Fig. 1 raises 
the temperature, and so accentuates the falling 
off of the moment, unless there is a counter- 
balancing upward bending of the curves of 


constant x in the M—(1/H”) plane. (According 
to (3) these curves are straight lines with 
downward slope for all x.) If there is an actual 
near equality of moments at these two field 
strengths, the question immediately arises as to 
whether the moment remains constant at inter- 
vening fields, or instead passes through a 
maximum somewhere between 105 and 24,000 
gauss. It would obviously be desirable to have 
measurements at intermediate fields. Unfortu- 
nately such measurements are not available and 
do not appear very feasible with de Haas and 
Wiersma’s apparatus since at low temperatures, 
a weak solenoid capable of producing only low 
“‘measuring fields’’ is employed, while the very 
large magnet used at helium temperatures is only 
fit for high fields. So the gap in data apparently 
cannot be closed. However, even a small increase 
in the measuring field above 105 gauss would be 
illuminating. If there is an intervening maximum, 
which we rather doubt, Eq. (1) obviously cannot 
be used to determine the thermodynamic temper- 
ature, since its derivation requires that the 
moment remain constant over an interval rather 
than be merely fortuitously equal at two field 
strengths. If, however, the high temperature 


| x-2 
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value of the moment is being approached 
asymptotically, (1) might be approximately 
correct. In this event, M would needs be slightly 
lower at 105 than at 24,000 gauss, but the 
deviation might conceivably be so much lower 
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than our estimate of 18 percent as to escape 
detection. Two other causes besides the neglected 
higher order members which might reduce the 
apparent discrepancy are (a) a g factor slightly 
less than 2, inasmuch as the second-order dipole 
perturbation is proportional to g‘, and (b) experi- 
mental errors in the relative calibrations of the 
magnets used at high and low temperatures. The 
error due to the latter cause might amount to 5 
percent or so. 

It is clear that additional theoretical and 
experimental work is much to be desired, because 
of the light it will shed on the mooted question 
of how far dipole-dipole interaction can be 
represented by an inner field of the Lorentz type. 
If the Lorentz hypothesis were strictly correct, 
then Eqs. (6), (7), (8) would always be fulfilled, 
and (1) would be valid at all temperatures. The 
latter is, of course, impossible, but the Lorentz 
assumption is to a certain extent the equivalent 
of the simplest form of the Heisenberg theory of 
magnetism, in which the constant of pro- 
portionality in the Weiss molecular field is 
independent of temperature, but which works 
much better than a second-order (so-called 
“Gaussian”) perturbation calculation. All the 
same, it is impossible to have a power series 
in t/T with nonvanishing coefficients always 
summing to zero, but this situation is what 
would be required if the perturbations of various 
orders are to neutralize each other so as to 
restore the accurate validity of (1). 

It seemed worth while to publish the calcu- 
lations in their present form, inclusive only of 
second order dipole effects, not merely because a 
third order computation would be quite long, but 
primarily because an experimental test of the 
second order corrections is possible by using 
much lower initial demagnetizing fields than 
24,000 gauss. The lower this initial field Ho, the 
smaller the final value of +/T and so by de- 
creasing Hy the ambiguity due to the neglected 
members becomes asymptotically small. Table I 
shows that the order of magnitude of the critical 
field strength H, at which appreciable deviations 
of M/M, from unity first appear is not greatly 
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influenced by the size of the initial field, so that 
experiments with lower Hy remain feasible. The 
advantage is that H, can be computed with 
increasing certainty, the smaller Ho. Unfortu- 
nately the only published measurement with 
H,)=12,000 gauss is for a very weak final field, 
5.23 gauss. The discrepancy of about 25 percent 
between the observed and calculated value of M 
at this field is far from excessive for a possible 
influence of higher order members, since 7/T is 
still not small enough to make the series converge 
at all adequately,® and, after all, a deviation of 
25 percent is small compared with the total 
reduction in moment by a factor 20. It would be 
of interest if measurements were available with 
still smaller initial fields, say 2000 to 7000 gauss. 

It may be noted that in his most interesting 
experiments on paramagnetic relaxation, Gorter® 
finds that in caesium titanium alum the ab- 
sorption per titanic ion does not increase greatly 
on dilution with nonmagnetic ingredients. This 
he interprets as showing that some other. factor 
than coupling between paramagnetic ions inter- 
feres with the free orientation of the spin. 
Because of Kramers’ theorem on double de- 
generacy,? and because no adequate disturbing 
nuclear moments can be invoked, it is hard to see 
what this mysterious factor can be. If it really 
exists, it perhaps modifies the present calculations 
for extremely low temperatures, in some unknown 
way. 

The writer is much indebted to Dr. C. J. 
Gorter for valuable discussion. 


5 That our neglected higher order terms are very im- 
portant even in demagnetization from 12,000 gauss can 
be seen by integrating Hebb and Purcell’s (reference 3) 
specific heat formula (13) so as to obtain an expression 
for the entropy. This formula is more special than our 
(2) in that it applies only to very weak fields, but it does 
include third and fourth as well as second order terms in 
the dipole interaction. The inclusion of these two extra 
orders makes the value of (7'/r) appropriate to zero field 
strength imaginary, an absurdity, rather than 2.68 as 
given in the table. The situation is reminiscent of that in 
the Heisenberg theory of ferromagnetism, where Opechow- 
ski shows that an imaginary Curie point is obtained if 
extra terms are included in place of stopping with the 
usual approximations (Physica 4, 181 (1937)). 

6 Cf. C. J. Gorter and R. de L. Kronig, Physica 3, 1018 
(1936). 
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A general development of the Gibbs phase integral for a system of chemically saturated 
molecules gives equations for the thermodynamic properties of the system. The equations 
predict the usual condensation phenomenon of the gas to form a condensed phase with a surface 
tension at temperatures below a characteristic temperature 7;,. At 7, the surface tension of the 
condensed phase is zero, and this phase undergoes a discontinuous change without heat of transi- 
tion. Between 7,, and the true critical temperature 7., there exists a finite volume interval for 
which P and the Gibbs free energy F is independent of the volume. Above 7, the usual P— V 
diagram predicted above the critical point is found, with no volume for which (8P/dV)7r is zero. 





THE GENERAL METHOD 


N a preceding paper! one of the authors has 

shown that the method commonly used in 
statistical mechanics to obtain the second virial 
coefficient, if followed through more completely, 
leads to equations predicting the existence of a 
condensed phase. In this paper we propose to 
develop a more concise method of handling the 
equations, and to show that the equations also 
predict the existence of some unexpected phe- 
nomena at the critical point. The method is a 
development, without arbitrary physical or 
mathematical assumptions, of the equations for 
the thermodynamic properties in terms of several 
series of integrals. In the temperature pressure 
range for which the liquid is stable the integral 
playing the predominant role cannot be numeri- 
cally evaluated. 

The system with which we wish to work must 
be limited to have certain properties if the 
development is to be applicable. 

Three of these limitations are absolutely 
essential. They are: first, that the masses of the 
particles making up the system are large enough 
so that the classical laws of mechanics are 
applicable at the temperatures and volumes we 
wish to investigate; second, that the mutual 
potential between particle pairs falls off with the 
fourth or higher power of the distance between 
the pairs; and third, that the total potential 
energy of the system is representable as the sum 
of the potentials between pairs of particles which 
depend only on the differences of the coordinates 





*Contribution from the chemical laboratory of the 
Johns Hopkins University. 
' Joseph E. Mayer, J. Chem. Phys. 5, 67 (1937). 
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of the pairs. The last limitation is equivalent 
to restricting the system to one composed of 
chemically saturated molecules incapable of 
forming true valence bonds with each other. 
It does not mean, of course, that we neglect 
higher than two-body collisions: on the contrary, 
we specifically treat multiple interactions, but 
we regard these as being made up of the sum 
of the interactions between all conceivable pairs 
of molecules. 

Certain other unessential limitations are im- 
posed on the system, limitations which are not 
necessary for the applicability of the general 
method, and which could be eliminated one by 
one at the cost of more complicated equations 
and a more complicated terminology and nota- 
tion. These limitations are: that the system 
consist of N identical molecules; that N is 
large enough so that terms of the order of N-! or 
smaller can be neglected compared to unity; 
and that the individual molecules have three 
translational degrees of freedom only. This latter 
limitation means that internal degrees of freedom 
are absent or neglected, and that the mutual 
potential of a pair of molecules depends on their 
distance apart only, and not on mutual orienta- 
tion angles. 

The problem is attacked by means of the 
normalized phase integral, which is 1/N! 45% 
times the integral (over the 3N momenta and 3N 
coordinates of N particles), of the exponent of 
minus the Hamiltonian divided by kT. The 
integration over the momentum space can be 
performed immediately if Cartesian coordinates 
are used. There remains the integration over the 
coordinate or configuration space. The integrand 
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is the exponent of the negative of the potential 
energy divided by kT. Since this potential energy 
is the sum of the potential energies of the 
N(N-—1)/2 possible pairs, the integrand can be 
written as the product of this number of terms. 
As in the usual method of obtaining the second 
virial coefficient, a function of the distance 7;; 
between molecules 7 and j, the quantity f(7:;) +1 
(which defines the function f), is substituted for 
each exponential. The product is then expanded. 
The steps are indicated by the following equa- 
tions. 


1 
o=-— Jf ; » ferme DITA D,- » «dn 


2armkT \ *%!2 Q, 
-(——) oa) 
h? N! 


a= ff. e s ferontdry ‘ -dtn 
= we II en ?(rig kT 7,- . -dtn 
N2>i>j21 








(2) 
=f. . 2 ec ° -dtN 
=f f-- +2 fitLfifert---J 
Xdr1: ‘ -dtyn, 

with dr,=dxdydzi, (3) 

fis=f(rizs) = (EO FP —1). (4) 


A plot of f(r) for A/r® law of attraction and 
B/r® repulsion potential at a temperature corre- 
sponding to the stable existence of a liquid is 
shown in Fig. 1. Since f drops rapidly to zero as 
the value of its argument 7 increases, the value 
of any term in the sum of Eq. (2) which involves 
a certain f;;, is zero if the distance r;; between the 
two molecules i and j is large. This term con- 
tributes to the total phase integral only in that 
part of configuration space for which the two 
molecules i and j are close together. For this 
reason we shall refer to the two molecules 7 and 7 
as being “‘bound’”’ to each other in any term 
containing the function f;;. If, in any particular 
term, a molecule 7 is not bound to any other 
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molecule (i.e., its index does not occur anywhere 
as a subscript in the product) we shall refer to 
it as a free molecule, or as being in a cluster 
consisting of one molecule. If, in a given term, 
the molecule 7 is bound to j, and neither 7 nor j 
is bound to any other molecule, we shall say that 
i and j are in a cluster of two molecules. If 7 is 
bound to j and j bound to k, whether or not 7 is 
bound to & we shall say 7, j7, k form a cluster of 
three, etc. Or in general, if / molecules are 
bound directly or indirectly to each other, but 
none of them is bound to any molecule not 
occurring in the group we shall say that they 
form acluster of / molecules in that particular term. 

We now propose to collect together all terms 
characterized by having the same number, m, 
of clusters of one molecule (free molecules), mz» 
clusters of two molecules, . . ., m, clusters of / 
molecules each. This method of collecting terms 
is purely arbitrary and by no means the only 
logical method of procedure from this point; it 
appears, however, to be justified, a posteriori, by 
the simplicity with which important results can 
be deduced from the procedure. Other equally 
logical procedures have failed so far to lead to 
interpretable results. 

The physical significance of: this particular 
grouping of terms is the following. All terms 
characterized by the same set of numbers, m1, 
of clusters of size 1, have in common the property 
that they contribute to the total phase integral 
Q only in that part of the configuration space 
where mz pairs of molecules are spatially close to 
their mates, . . ., mclusters of / molecules each 
are close to the other molecules of the clusters, 
etc., whereas the contribution to the integral is 
independent of the coordinates of the my, free 
molecules. 
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STATISTICAL MECHANICS 


The quantity 5), the integral over a cluster of 
size 1, is defined by the following equation : 


nS S 


(all products consistent with a single cluster). 


> I fijdry---dry (5) 


l>i>j>1 


The normalization factor 1//! is merely for 
convenience. The integration is extended over 
the configuration volumes of all of the particles. 
The integrand is the sum of all possible products 
of the f’s which are consistent with the existence 
of a single cluster. For small values of /, and for 
volumes of macroscopic dimensions, }; is inde- 
pendent of the volume, due to the property of f 
that f(*)=0. For values of / of the order of 
the number of particles in the system, the volume 
below which }, becomes (significantly) volume 
dependent we shall later see is the volume of 
the condensed phase. It can be seen that in any 
term, the integration over the last particle, the 
I’th particle, in every cluster yields the volume, 
V, as a factor in the term, which is the reason 
we have defined }; with V in the denominator. 
If 6; is taken by definition to be unity this is 
also true of the m clusters of single molecules. 

It is then easily seen that the expression for 
Q,/N! becomes 





Q: (Vd,)™ (Nvb,)™ 

=) Il——=L I : (6) 
N! me m,! mi it m,! 

= lmi=N = lmi|=N 


with o> (the volume per molecule). (7) 


THE GENERAL EQUATIONS 


The logarithm of the sum Q,/N! of Eq. (6) 
may be approximated with sufficient accuracy by 
the logarithm of the largest term,! giving Eq. (9). 
This method is only valid if all the },’s are 
positive. Professor Born has pointed out to the 
authors that the function Q,/N!, which can be 
seen to be the coefficient of Z” in the expansion 
of the exponent }°vb,Z', can be expressed as a 
Cauchy integral and evaluated by means of the 
well-known method of steepest descents, leading 
also to Eq. (9). We believe that Professor Born 
will publish this proof shortly along with several 
other comments on the method of the present 
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article. Another method of obtaining the asymp- 
totic Eq. (9) is discussed in the appendix of 
this article. 

If a function Z(vb,, vbe, ---, vb;, «+ +) is defined 
as the lowest positive real root of the equation 
((A4) of the appendix), 


Gi(vb;) = > lvb,Z t= i. (8) 


l21 


and the functions G; are defined as 


G;(vb,) = >olivb,Z', (9) 


21 


so that G;=1, then we find that, (Eq. (A57) of 
the appendix), 


‘. QO; 
7 n "oe In Z, (10) 


plus terms of order of magnitude 1/J, as long as 
Q, is positive and real for all large values of N, 
a condition which it obviously always fulfills. 

The Helmholtz*free energy A, the free energy 
at constant volume, or the work function of 
Lewis and Randall, is derived immediately from 
the phase integral as 


A=-—kTInQ 
and from (1), 


A= rz] mn (— “yk 
2amkT N N! 


which, with (10) becomes 





A=RTT[ In (h2/2rmkT)!—Go(vb,) + In Z]. (11) 


The negative of the derivative of this with re- 
spect to the volume gives the pressure 


P=-—(0A/dV)r=kT (0/dv)[Go(vb,) — In Z] (12) 
=(RT/V)[Go(vb:1) +Goo(vbi) ], 
where 


G;»(vb1) = Dlvb,Z"(d In b/d lnv). —(13) 


l>1 


(the details of the differentiation are given in the 
appendix, preceding Eq. (A18)). The free energy 
at constant pressure, F, the Gibbs, and Lewis 
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and Randall, free energy, is given by 


F=A+PV 


4 
= RT[In (h?/2rmkT)'+1n Z+Go, |. vale 
We shall defer a more detailed discussion of 
the nature of these equations until the properties 
of the arguments }b; have been discussed more 
fully. However, one interpretation will be intro- 
duced at this point. The original Eq. (6) for Q, 
can be handled more simply than by the method 
of the appendix if all the },’s are positive. In this 
case the logarithm of the total sum of Q,/N! can 
be shown to become asymptotically equal to the 
logarithm of the largest term as N increases 
indefinitely. Now the largest term is charac- 
terized by the set of numbers m,*, the number of 
clusters of size 1, with Jm,* the number of 
molecules in clusters of size /, and /m,*/N the 
fraction of material in clusters of size 1. lm,*/N 
can be shown to be equal to /vb,Z' so that the /’th 
term in the sum G, of Eq. (8) is the (most 
probable) fraction of material in clusters of / 
molecules. The /’th term in Gp is the number of 
clusters of / molecules each, divided by the total 
number of molecules, and from the Eq. (12) for 
the pressure it can be seen that each cluster 
contributes exactly as much to the pressure as a 
single molecule, (if Go, is zero). For large values 
of v, 0 In b;,/d In v is zero and Go, must be zero. 
The volume below which Go, differs from zero 
will be called vy and later be shown to be that 
of the condensed phase. 


EVALUATION OF THE ),’s 


We shall now investigate the integrals over the 
clusters, the b,’s of Eq. (5). It can be shown! that 
normalized irreducible integrals 8, can be so 
defined that the cluster integral b; is the sum of 
coefficients times products of the 6;’s raised to 
various powers. 

The first three normalized irreducible integrals 
are 


6=(1/V) f a f Anr*f(r)dr, (15) 


Bo=(1/V2!) f f Dihtddidbesn, (16) 
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a= (1/V3)) f ff [3fafsfafat 6fafofafufe 
+ fasfsoferfarfsefs:ldtidredrsdry. (17) 


The coefficients 3 and 6 in $83 occur because 
permutation of the indices gives, respectively, 
3 and 6 different products leading to the same 
integrals. The general k’th irreducible integral is 
an integral over the configuration space of k+1 
particles. The integrand is the sum of all products 
of f between k+1 particles which connect every 
particle to every other particle by more than a 
single path, 


(4/0 _ IT figdri- + -dri41. 
a=ci/vey f f f= Bs “*48) 


Products having all particles more than singly connected. 


If n, is the power to which f; is raised in a 
certain term in 5; we have the condition that 


kn, =l1—-1. 
k 


The coefficient with which a term, defined by 
the powers 7,;, occurs in b; is very simple, namely 


1-*TIl"*/n;, I, 
k 


The proof that this is the correct coefficient 
caused the authors a great deal of difficulty. 
They are greatly indebted to Dr. Goeppert 
Mayer for the final completion of the proof. 
An indication of the general method of this 
proof is given in the second part of the appendix. 
Due to its complicated nature the proof is not 
given in detail. 
The equation for /?b; becomes 


rb.= > II (1B ,.)"*/nx!. (19) 
n k 
Sant 


Comparison with (6) shows that it is similar in 
form to the equation for Q,/N!, indeed it is 
identical in form with the Eqs. (A58) and (A36) 
of the appendix defining the function F,(M, x). 
We may use Eq. (A59) of the appendix which, 
with the arbitrary introduction of the function 
p(l, B), is also an exact equation for b:, 


[0(8) }’n(8) 


= ’ (20) 
(2rG2(8))*2*/? pil, B) 
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where p(/, 8) is a function which takes erratic 
values for small values of / but which becomes 
asymptotically equal to unity as / becomes 
large. For bo(8) the equation, (A42), 


In bo(8) =Go(8) —In p (21) 


is obtained. p is determined by the £;’s in the 
same way that Z is determined by the vb,’s 
(Eq. 8), i-e., 


Eko rot =1, (22) 
k>1 
_ Go(8) = EBipt. (23) 
k>1 


For small values of / it is more convenient to 
use Eq. (19) to determine 0, for large values of / 
the form (20) is of course preferable. The con- 
dition for Eq. (20) to be valid is that Eq. (22) 
has a real positive root of p, the lowest real 
positive root being the correct one to use. 


DEVELOPMENTS FOR THE GAS 


Returning now to Eq. (8), we see that at 
infinite volumes, since }; is unity, Z=1/v, and 
Go=1. Using these values in (11), (12) and (14) 
the equations of a perfect gas are obtained, 
remembering that Go, is zero for large volumes. 

As long as the sum G, is absolutely convergent, 
(see the discussion in the appendix), Z can be 
developed as a power series in 1/v. Use of such a 
power series in (8) enables one to equate coeffi- 
cients, and if b; is expressed by (19), one obtains 


Z=(1/v) {1-8 —[Be— (B1?/2) Jo 
—[Bs—B:B2+(B1°/6) Ju: - =} 


“=yte Bae (24) 


This may, in turn, be used in (9) to calculate 
G(vb;) as a power series in descending powers 
of v. The following simple equations are obtained : 














Go(vb vin 1 ak 25 
o(v =" ~2 ae ; (25) 

4 he \i1 1 

A=RT in ( ) ——), au] (26) 
q 2rmkT/ ev & k+1 

. r hn? 31 

F=RT in ( ) -— Au (27) 
L 2rmkT/ vk 
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It is seen that —k®,/(k+1) is the (k+1)th 
virial coefficient. 

The developments were first made by algebraic 
evaluation of the first few terms. They can how- 
ever be proved by inserting the solution of Eq. 
(24) for Z in Eq. (8), and by manipulation of 
the sum obtained on expanding the exponential 
it can be shown that the coefficients of all 
powers of 1/v except the zeroth power are equal 
to zero. Professor Max Born has written to the 
authors that there exists a still simpler proof of 
the equations. The developments are, however, 
only valid if the sum G,(vb;) converges absolutely, 
so that the last terms are negligible in their 
effect on the determination of Z. 

From (25) we may evaluate (0P/dV)¢7 as 


(9P/8V)r= —(RT/V*)[1—DkBw*], (28) 


and if we define p as the reciprocal of the largest 
volume for which this becomes zero, it becomes 
the smallest positive root of the equation 
gi(8, p) —1=0, identically the Eq. (22) by which 
p was previously defined. We see, then, that the 
quantity p has the physical significance of being 
the reciprocal of the largest volume, v,, for which 
(0P/dV)7=0, 


o4™ 1/p (29) 


From Eq. (24) for Z when v=v,=1/p, and (21) 
for bo, we see that 


Z,=bc, v=0,=p-1. (30) 


By differentiating the logarithm of Z in (24) 
with respect to the logarithm of the volume, one 
obtains 


dln Z/d In v= >-kBw-* —1, (31) 


which, in view of the definition of v, as the 
largest root of this equation, is zero at v=v,, 
and negative at all higher volumes. It follows 
that for v>v,, bbZ<1, and the sums G;,(vb,) 
converge absolutely (see the discussion of the 
second section of the appendix). In the range of 
volumes, v, <v <«, then, the equations (24) 
to (31) are valid. 

The quantity Z has the dimensions of the 
reciprocal of a volume. From the equation for F, 
(14), it can be seen that for large volumes when 
Go, is zero, Z—! is the volume per molecule of a 
perfect gas having the same free energy F as the 















real gas in question. The quantity (vZ) =y, may 
be called the concentration activity coefficient 
of the gas. It is the dimensionless quantity by 
which the concentration of the gas must be 
multiplied to obtain the concentration of an 
ideal gas having the same free energy F. From 
(24) this is seen to have the value 


Ye=0Z = em =He0*, (32) 


CONDENSATION 


For volumes lower than v, the sums G;(vb) 
are not absolutely convergent, (see the discussion 
in the second section of the appendix), and the 
equations (24) to (32) will no longer be valid. In 
this range we may use (20) for b;, which, inserted 
in (8) leads to the equation for the determination 
of Z, 





v (b0Z) ‘p(8) 
=1, (33)? 
(2m)! 121 13(G2(8))*p(Z, B) 
Since p(l, 8) becomes asymptotically unity for 
large values of /, it can be readily seen that boZ 
cannot become more than infinitesimally greater 
than unity. For values of the volume v equal to, 


or less than v,, we can, to within terms of the 
order of magnitude of 1/N, write 


Z=1/bo, 9InZ/dv=0, v <vs. (34) 


By taking the derivative of Go(vb,)/v with 
respect to v one obtains 


0 [ Go(vb:) =, dlnZ 
= |- 





+Go(o)| (35) 


v vL ov 


If now 2; is defined as the volume above which 
Go, is essentially zero, and since v, is the volume 
below which @1n Z/dv is zero, in the volume 
range vy <v <v, we find that Go/v is a constant. 
From (12) for P we see that the pressure is a 
constant in this volume range. Since Z is inde- 


2 We have written G2(8) as though it were independent 
of /. In the case (A26) of the appendix for which the sums 
G;(8) are not absolutely convergent G2(8) will depend on /. 
For this reason we have written G2(8) under the summation 
sign. 

Since we have_previously found that Zbp=1 at v=1/p 
it follows that ym 1/19/2(G2(B))*p(1, 8) = (27), an equation 
which we should be able to verify independently, but have 
been unable to do. Essentially this equation, however, 
has been carefully checked numerically for a wide variety 
of different values of the (;’s. 
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pendent of v, and Go, zero in this volume range, 
we see from (14) that the free energy, F, is also 
independent of the volume. This constancy of P 
and F for v; <v <v, may be considered to be the 
thermodynamic criterion of a change of phase. 
v; is seen to be the volume per molecule of the 
condensed phase, and v, the volume per molecule 
of the saturated vapor. 

The decrease of volume from v, to vs is accom- 
panied by an increase of the fraction of material 
in the very large clusters, lvb,Z', from zero to 
almost unity,’ since the sum becomes divergent 
in this range (see the discussion in the appendix), 
but there is no corresponding increase in the 
fraction of material in the moderately large 
clusters ; this remains essentially zero during the 
whole process. The fraction of material in the 
small clusters decreases proportionally with v in 
this range. Even at v; the density of the small 
clusters, the number per unit volume, 0,2’, is 
the same as in the saturated vapor at v,. Insofar 
as the clusters have a physical meaning this 
states that the number of free molecules in the 
liquid per unit volume is equal to the number 
per unit volume in the vapor in equilibrium with 
the liquid. This same relationship of equal den- 
sities in liquid and saturated vapor holds for all 
the small clusters. 

With the knowledge that the volume of the 
saturated vapor per molecule is v,=1/p, and that 
Eq. (25) holds down to exactly this volume we 
may write the equation for the vapor pressure 
over the condensed phase as 


P,=kT p[i— 2 (k/k+ 1)Bxp*], 
(vsKvKv,). (36) 


Since Z=b,"! for the range vy <v <vs, and 
Go»(vb:) is zero, we may write from (14), for the 
free energy F of the saturated vapor, and of the 


3 The treatment of Fo(M, x) in the appendix is not fully 
valid without further discussion if G.(x) diverges. The 
method of steepest descents is also invalid in this case. An 
extension of the method of the appendix shows that the 
equations used in the article are correct, in this case, to 
the degree of accuracy to which they are used. The method 
of using the logarithm of the largest term shows the 
difficulty most clearly. In the previous treatment the 
numbers m, are treated as continuous, although they are 
obviously integers. If this simplification is not made, it 
becomes apparent that in the range of volumes vs=0=—2, 
for which condensation takes place, only one of the large 
clusters is present for which /= N(v,—v)/v,. The equations 
as used in this article are in no way affected. 
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liquid under its own vapor pressure, 
F,=RTT[ In (h?/2rmkT)}— In bo | 
= RT[In (h?/2mmkT)!—G)(8)+ In p], 
(vsQ vv). (37) 


THE CritTicaL Pornt 


The {;’s are all positive at very low tempera- 
tures. As the temperature increases the 8;’s of 
higher k values tend to become negative. As dis- 
cussed in the appendix, above a certain tem- 
perature, 7, there exists no real root of equation 
(22) for p. Above this temperature Eq. (20) for 
b; is invalid, and the complete argument that 
predicts the condensation phenomenon breaks 
down. One is immediately tempted to identify 
T, with the critical temperature. This identi- 
fication is made complete when we observe that 
p has also been defined as the reciprocal of the 
volume for which (0P/dV)7, as given by Eq. 
(28), becomes zero. Above T,., where no real 
root for p in Eq. (22) exists, there is no positive 
volume for which (0P/dV)+ is equal to zero. 

From Eq. (28) we may calculate (d?P/dV7?)7, 
and find 


(0°P/dV?)r=(RT/V*)[2—2>- kB wo * 
—Dk*Bw*], (v2). (38) 
At v,=1/p this becomes, since G,(8) = 1, 


(0°P/AV*)r, v=o, = —(RT/V*)G2(8), (v=vs). (39) 


As is shown (Eq. (A 60)), in the appendix G2(8) 
is zero at T=T,. 


CONTINUITY OR DISCONTINUITY AT THE 
CONDENSATION POINT 


For the sake of completeness we may return 
to the general Eq. (12) for the pressure P, and, 
limiting ourselves to v2 v;, for which Go, is zero, 
calculate (@P/dV)7 and (6?P/dV?)7 as 


aP od RT RT 1 
(—) =— —Gr(0b,) = -—__—__,, 
OV/ 7 AV V V? Go(vb:) 


(v2 vs), (40) 








a?P RT G;(vb:) — 3G2(vb1) 
(—) Pn , (v>v,). (41) 
OV2/ y3 


[G2 (vb;) Ny 
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If v>v, all the sums in these expressions are 
finite, and the expressions are more readily cal- 
culated by the equivalent equations (28) and 
(38). For v<v,, however, the sums diverge (in the 
sense discussed in the second section of the 
appendix) and G2(vb;) becomes of order of mag- 
nitude N, and G;(vb;) of order N?; both (@P/dV)r 
and (d?P/dV?)r are seen to be zero. 

It is more difficult to ascertain from Eq. (42) 
whether (0?P/dV?)7 becomes “‘infinite’’ (i.e., of 
order of magnitude JV), or not, as the volume v, 
is crossed. Eq. (39), however, answers this 
question. If the sums G,(8) and G2(8) converge, 
(@?P/dV?)¢ is seen to be finite at v=v,; whereas 
if they diverge, this function becomes “‘infinite’’ 
as the volume crosses v,. 

It appears then that the question of whether 
(@P/dV)r goes continuously or discontinuously 
to zero at v=v, is determined by whether the 
sums G;(8) and G2(8) converge or diverge. If the 
sums are convergent (0P/dV) 7 goes continuously 
to zero as v decreases through v,, and we must 
assume that the vapor cannot be supersaturated 
in the sense that a higher pressure than P, could 
be observed. If, however, the sums are divergent, 
there is a sudden break in the shape of the P-V 
curve at v,, and supersaturation of the vapor 
would be expected. 


SURFACE TENSION AND THE TEMPERATURE 7, 


For large values of k we have been unable te 
find a method of calculating the 8;’s numerically. 
It is, however, fairly obvious that they must 
obey the relationship 


(lim kR- ©) In Bg = In Bo, (42) 


in which 8, of the dimensions of a volume, 
represents the contribution as a factor to 8; of 
the integration over an average particle. It is 
also fairly obvious that in the important terms 
of the products of fi;’s of Eq. (18) a number of 
molecules proportional to k! must have fewer 
“‘bonds”’ to neighbors than the average particles. 
These are the particles which are physically on 
the surface of the mechanical frame which may 
be associated to each product of f’s. We may 
now amend (42) to 


In B.=k In Bo—yk'+1n a(k, Bo) —nInk, (43) 
B.=Boke-*a(k, Bo)/k", (44) 
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where 8» is a positive quantity with the dimen- 
sions of volume, the dimensionless y may be 
positive or negative, a(k, Bo) is presumed to be 
independent of k for large values of k, (a property 
which we shall not really use), and k-” is in- 
troduced merely because it seems rather probable 
that some such term should be present. The only 
portions of (44) which are of importance are the 
already justified Bo and y. 

Now for low values of the temperature, 7, y 
must be positive, the contribution of the “‘sur- 
face’’ molecules to 8; must be less than the con- 
tribution of the ‘‘inner’’ molecules, and we must 
expect that dy/dT <0. Let us assume y>0 for 
T<T,,, dy/dT <0, and y=0 at T=T,,. We now 
propose to show that, for T<Tm, kTy/(3620/)! 
is the surface tension, S, or surface free energy 
of the condensed phase. 

The case (44) is the same as the case II dis- 
cussed in the appendix under the title “‘Con- 
vergence of the Auxiliary Functions.”’ Below 7», 
with y positive the sums G(8) and G,(8) diverge 
in the sense used in that section. For tempera- 
tures above 7,, the sums converge.‘ We have 
already seen in the preceding section that the 
cases of divergent and convergent G2(8) cor- 
respond respectively to discontinuous and con- 
tinuous transitions of (@P/dV)7r to zero at the 
volume v, of the saturated vapor. The existence 
of a positive surface tension below 7, assures 
the instability of the small clusters; so an 
arbitrary restriction to integration over that 
region of phase space where no large numbers of 
molecules are close together (corresponding to 
the arbitrary assumption that the 5,’s of large / 
values are zero or small) would result in a P-V 
curve showing higher values of P than P, for 
v<v,. 

The proof that k7Ty/(367v7)! is the surface 
tension is as follows. In the appendix it is shown 
that p, (the same as z of Eq. (A25II)) depends on 


4 Actually, although the sums certainly converge if y is 
negative, it is not necessarily true that they diverge for all 
positive values of y, and they would not if Yak, Bo)k™™ 
were greater than unity. We should define 7,, as the 
temperature at which the sums change from divergent 
to convergent, which could be at a finite positive value 
of y. This would change nothing in the subsequent argu- 
ment except that the surface tension S would extrapolate 
to a finite value at T7,, and would become discontinuously 
zero at this temperature. This appears not to be the case 
experimentally. 
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l, the upper limit of summation, if the sums are 
divergent, and indeed, (A29II), 


p=Bo ter!" (45) 


If this value of p is inserted in Eq. (20) for b,, 
with (21) for bo, we find a multiplicative factor 
e-v'! in b;. Now consider the chemical reaction 
in which clusters of size / dissociate to form 
clusters of size ]/n. The concentration of clusters 
of size / is given by c,=},Z', since vb,Z' gives the 
number of such clusters divided by N. The con- 
centration of clusters of size 1/n is byj,Z'!". In- 
serting these values in the reaction constant 
K=Cijn"/ci, and using the relationship AF°= 
—kT\n K we find, upon neglecting lower order 
terms that 


AF° = —kTyl(1—n!), (46) 


where AF® is the difference in free energy (at 
hypothetical unit concentration) between 1 
clusters of size //n and one cluster of size /. This 
difference in free energy must be the product 
of the surface tension S multiplied by the increase 
of surface in the reaction. The area of the surface 
of a (spherical) drop of / molecules will be 
(367) *(v,l)!, since vs is the volume per molecule 
of the condensed phase. The area increase is then 
seen to be (367v//*)'(n'—1) and one obtains 
with (46) the equation for the surface tension in 
terms of vy, 


S=kTy/(3600/)!. (47) 


In this section we have inserted the quantity 
y into the equation for 6, as a reasonable as- 
sumption and have shown that its reality is in 
agreement with the experimental observation of 
the existence of a surface tension. For our equa- 
tions to have a physical meaning, however, it is 
important that the quantity a(k, Bo) be positive 
for all large values of k, as long as y is positive. 
In this case G2(8) must be very large and 
positive. At 7, the quantity y is zero, and 
although it may and presumably does become 
negative for higher temperatures, because of the 
necessary convergence of the sums G,(8) and 
G2(8) above T,,, the quantity p does not depend 
on the upper limit of summation /, and In 0; will 
contain no term proportional to /', Therefore, 
the surface tension S is zero for T>T,,. How- 





ev 
tel 


Ty, 


ter 
an 
dis 
the 


tan 
tur 
by 


clu: 


the 
exis 
Bel 
ten: 
disc 
whi 
is a 
curt 
but 
Bet 
ther 


Ir 
dyn: 
it w 
as f 
Alth 
any 
large 
more 


By. L 


‘ll 


w- 





STATISTICAL MECHANICS 


ever, G2(8) is positive at T=T7,,, and so the 
temperature 7, the true critical temperature, 
at which G,(8)=0, must lie higher than 7,,, 
Ta<le 

Thus 7, the temperature at which the surface 
tension of the condensed phase becomes zero, 
and where the meniscus between liquid and gas 
disappears, is not identical with, but is lower 
than the true critical temperature T., the tem- 
perature at which the P-V curve has a horizontal 
tangent for only one volume. The two tempera- 
tures have been hitherto regarded as identical 
by most theoretical workers. 

To recapitulate: we are forced to the con- 
clusion that there exist two characteristic tem- 
peratures, T, and 7,,, with T.>T7,,. Above T, 
the P—V curve is nowhere horizontal and there 
exists no difference between gas and liquid. 
Below T7;,, the condensed phase has a surface 
tension and the P—V curve shows (at least) one 
discontinuity in slope at the volume v,, below 
which it is horizontal. Between T7,, and 7, there 
is a finite temperature range in which the P-V 
curve is horizontal over a finite volume interval, 
but has no discontinuous change in_ slope. 
Between T,, and 7, and of course also above T,, 
there exists no surface tension. 


TEMPERATURE DERIVATIVES 


In order to obtain equations for the thermo- 
dynamic quantities S, E, H, and (8P/dT)y, etc., 
it will be necessary to express the functions 6, 
as functions of T and to differentiate them. 
Although this might be done formally without 
any special assumption as to the form of 6, for 
large values of k, certain relationships become 
more obvious if Eqs. (43) and (44) are used for 
6,. Using the definitions 


—d In Bo/d\ln T=Bo, (48) 
—dy/d\n T=Cy, (49) 
d\n alk, Bo)/d In T=B,, (50) 


we may write 
d\n B,/d In T= —RkBo+k!Cot+ Bx. (51) 


By the use of (51) in Eq. (26) for A with 
v2v, one obtains 
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k 
(— “) =- S=—— “—an+n| 2 Ds “TT eal 


ki 


Bi By 
—Cod, —bw*-> —r] (52) 
k k+1 k k+1 


or 
2rmkTe\ } By. 
mae Cr) ev+)>> ——v* 
h? k k+1 
kip, BiB 
yk y-* |. 
k k+1 k k+1. (53) 


Using E=A+TS=A-—T(0A/0T)y, and H=E 
+PV, with (25) for PV one obtains 
3 RB; kiB,, 
E= RI] -— BE ——y kt Co>. ——y-* 
2 k+1 k k+1 


BB: 


+z, r+ (54) 
e k+1 





ett Od 








kip, 
v*+Cy>, 
k k+1 


y—* 








5 
H=RI| 
2 


BiB: 
+E] (55) 
k k+1 


Similarly (25), if differentiated, yields 








wi 
(— =) == “9-4 BS 
k+1 k k+1 
R*!8B,. RB By - 
Cod, 0 -5 v4], (56) 
k+1 k k+1 


For dP,,/dT we may use 
dP,/dT =(0P/8T),=»,+(0P/OV)r, v=v,(0Vs/AT) 


but since (0P/dV)r, .-», is zero, we find that 
dP,/dT is equal to (@P/dT)y when v=»,.5 


5 This apparently high-handed neglect of the ‘“‘infinite’’ 
value of (d?P/d V?)r when the sum G2(8) diverges is justified 
by the fact that G2(8), and hence (6?P/dV?)7, is actually 
finite in a mathematical sense, although very large. We 
have used the expression ‘‘discontinuous” for (@P/dV)r at 
v=v, in the physical sense since magnitudes of order N 
are pretty big. This is, of course, mathematically unjusti- 
fiable, and all our equations are strictly speaking contin- 
uous. The equation (56) for (@P/dT)y and (57) for dP,/dT 
lead to different values for these two quantities if v is 
infinitesimally greater than vs. We shall discuss the values 
of the thermodynamic properties of the super-saturated 
vapor in the neighborhood of v=2, in the following section. 
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Using the condition that G,(@)=1 at v=v, 
=p, and that k?/(k+1)=k—k/(k+1) we 
obtain 


R RBx 
=| 143) (1-5 *) 
V;, k k+1 


5/38, , RB iB: 


~O5 pt— 
k k+1 k k+1 


P, 


dT 














at (57) 


Finally, by noting that G:(8)=1 at v=v,=p™ 
so that dG,/dT=0, we obtain for d In v,/d In T, 


d\n v,/d In T= —Bo+[G2(8) 
XLCoVk**Bip' + VRBiBip*], (58) 
or, since G2(8) is “infinite” when T<Ty, 
dinv./dln T=—By, (T<Tn). (59) 


It is seen from Eq. (58) for d In v,/d In T that 
this quantity increases discontinuously as T is 
increased through the value T=T,,. dP,/dT 
is, however, seen to be continuous at 7=T>,. 
Further differentiation shows that d?P,/dT? is 
discontinuous at 7,,. Since dP,/dT is continuous 
at 7, there is no discontinuous change in the 
heat content of the condensed phase at this 
temperature; therefore no true change of phase 
occurs at Tn. 

At T=T;, since G2(8) is zero it is seen from 
(58) that dlInv,/dln T is infinity. Numerical 
calculations of the 6;’s at this temperature have 
shown that the derivative approaches infinity 
from the negative side as is to be expected. 


EQUATIONS OF THE SUPERSATURATED VAPOR 


Equation (28) for (@P/dV)7 and (56) for 
(8P/dT)y show discontinuities at v=v, if T<T» 
and the sums G;,(8) and G2(8) are divergent. The 
point v=B,o~!, Bo/v=1 lies on the circle of con- 
vergence of the sum G;,() if that sum is regarded 
as having an infinite upper limit of summation. 
The value of v, is infinitesimally smaller than Bo 
so that 8o/v, is actually larger than unity by a 
very small amount. If we set Bo/v exactly equal 
to unity (a value which it has when 7 is negligibly 
larger than v,) we obtain a value of G;(8) which 
may be used in (28) and (56) for the super- 
saturated vapor in the neighborhood of the 
volume v,. Using (44) for B;, this leads to 
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gi(8,v')= Do RB 
k 


Bo/ v=1 


=Z alk, Bo)e-7**'Ri-" =G9(8). (60) 


The symbol G;°(8) has been used to indicate the 
sum g;(8,v~') when B)/v=1, G;,(8) being the 
value when G,;=1 which for divergent sums, 
T<Tm, gives Bo/v=it+e. It is obvious that 
Go(8) =Go(8) always, but G,°(8)<1 if T<Tn. 
Since @ In g;(8, v~')/d In v= —g;.1(8, v-) we may 
develop g;(8, v') in terms of a semi-convergent 
Taylor’s series as 


gi(8, 0") = [Inv—In», » 


X (0/0 In v”)g;(B, 0") v=, 
= DL (— In vp)’G?;,(8) J/v!. (61) 


v>0 


This equation could be used for values of 
v<v,=p~'! with any of the equations (25) to (32) 
or (52) to (59) to calculate the thermodynamic 
functions of the supersaturated vapor. We shall 
not carry out the obvious operations. 

It is interesting to note that the sum g;,(@, v7) 
is itself semi-convergent for values of v smaller 
than v,, that is the first terms of g.(8, v~') con- 
verge, and it is only after a comparatively large 
value of the index of summation of k has been 
reached that the terms commence to grow larger 
again. Using (44) for 6, and differentiating 
B.(v,/v)* with respect to k, we find that the term 
of minimum value is given by the equation 


1dlna(Bo,k) n 2¥ Vs 
— ——_—_—__—--—-—+4+In—=0. (62) 
k adlnk k 3k v 


For large values of k this is approximately 
8 7 


k=— ——___, (63) 
27 [In (v,/v) ]* 


which should give approximately the size of 
droplet which forms least easily at a volume 
v<v,. After a droplet of this size has once been 
formed in a supersaturated vapor of volume 2, 
the drop should increase rapidly in size until the 
pressure attains the value P,. 
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CONCLUSION 


Except for Eq. (44) assuming a special form 
for 8, no assumptions have been introduced into 
the general deductions other than those listed 
in the third and fourth paragraphs of the intro- 
duction. The assumption in (44) that In p; 
varies linearly with k for large k values can 
scarcely be questioned. The addition of a term 
proportional to k! appears, a priori, to be 
reasonable, and is later justified by its appearance 
in the final equations as the surface free energy. 

We conclude that a liquid-saturated vapor 
system has two characteristic temperatures, T,, 
and 7., instead of one critical temperature as 
usually assumed. Above T, and below T,, the 
properties correspond to those predicted classi- 
cally above and below the critical point, but 
between T,, and 7. for a certain volume range 
we predict properties which are not given by any 
classical description of the critical phenomenon. 
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The authors will show in a subsequent paper that 
published experimental results amply confirm 
this conclusion. 

Numerical evaluation of the first few £;’s 
shows that the experimental properties of the 
critical point, especially its relative position with 
respect to the temperature for which the second 
virial coefficient is zero (the Boyle point), and the 
value of PV/RT, are in agreement with the 
assumption that only #; and #2 play a role at 
these temperatures and pressures. It is probable 
that a calculation of 81, B2 and 83, which is not 
impracticable, would suffice to cover the range 
above T,,,. As yet no theoretical calculation of Bo 
appears to be practicable. Until this is accom- 
plished, no numerical values below T,, can be 
theoretically calculated. 

It is possible, by means of the formal intro- 
duction of new quantities, to extend the equa- 
tions to volumes below v;. This appears to be 
unfruitful at the present time. 


APPEND1X 


The auxiliary functions 


We shall first investigate the properties of a set of 
functions which occur repeatedly in the article. The func- 
tions g; of the variables x;, x2, ---, x,---, g will be defined 
as 


gi(*1, si Mad Q) = Lvixyg’, (Al) 
v21 


so that gg; /Oq = gj+1. (A2) 


If the value of the lowest real root of the equation for gq, 
g:=1, be designated as z, then this z is a function of the 
arguments x,; and if this z replaces g in the functions g; 
they then become new functions.G; of the arguments x, 
alone. 


G,(x) = > vix,2”, (A3) 
v21 

Gi(x) =) vx,2"=1 defines z, (A4) 
v21 


where the symbol x is used to represent all the x,’s together. 
Now 


IG;/dx, = Ag;/dx,+ (A In z/dx,) -q(dg;/Aq) quz, (AS) 
= viz” +Gj.1(0 In 2/dx,). 


Since 0G,/dx, =0 = v2” +G,.0 In 2/dx,, (A6) 
0 In 2/dx, = —v2”/Ge, (A7) 
and from (A5) that 

0G;/dx, = v2”(vi1—G;41/Gs), 


it follows that 


(A8) 








8? In 2/dx,? = (v2”)2(2vG2—G;3) /G:3. (A9) 
Also 8(Go—In 2) /dx, =z", (A10) 
02(Go—In 2) /dx,? = — (v2")?/Ge. (A11) 


If the x,’s are all functions of one variable y, x,=x,(y) 
from (A8) 


0G;/dy= >. x,dG;/adx,(d In x,/dy) 


= > vix,2”(0 In x,/dy) : 
— (Gj41/G2) D0 »x,2"(9 In x,/dy), (A12) 


8 In 2/ay= —[ > vx,2"( In x,/dy) 1/G2. (A13) 


It is seen from (A12) that 9G,/dy is identically zero. 

We shall have use for the special case that the x,’s are 
explicitly written in the form x,=vb,, but the b,’s may also 
be functions of v. If x,=vb,, 


d In x,/d In v=1+ In 3, /d In 0, (A14) 
and we define 


Gjv= >, vivb,2”(8 In b,/a In 2), (A15) 


then from (A12) and (A13) 
0G;/d In v= Gj— Gj41/GetGj »—Gj41/G2Gi», (A16) 

8 In 2/8 In v= —(1+G;»)/G2, (A17) 
0(Go—In z)/d In v=Go+Gov. 
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In the special case that all x,’s are proportionatl to u, 
p p 


i.e., if d In x,/d In u=1 (A19) 
0G;/d In u=G;—Gj41/Go, (A20) 
0 In 2/0 In u= —1/Ge. (A21) 


These functions play a role in handling the properties of 
the function Fo(M, x) which will be defined later. 


Convergence of the auxiliary functions 


In the definitions (A3) and (A4) of G; and z, the sums 
were written as if v varied from unity to infinity. We shall, 
however, be interested in cases in which the summation 
is explicitly limited to 1 < »v < M with M an extremely large 
number, although the x,’s are defined to »= ~. Actually 
we shall only be interested in cases for which 

( “ ) In x,=v In xo, (A22) 


ve? Oo 


and in two different modifications of this dependence of 
x, on v, namely 
(A231) 


(A23I1) 


in which xo, @ and a are positive, m need not necessarily 
be an integer, and da(xo, v)/dv approaches zero for very 
large values of ». 

We shall refer to the cases (A23I) and (A23II) as cases 
I and II, respectively. 


Xp =ae(Xo, v)xX0"/v™, 


and x, =a(Xo, v)xo’a”?/v”, 


v=M 
Now if Zz a(xo, v)/v™ 1 >1, (A241) 
v=] 
v=M 
or >Y a(xo, v)a’?/y™1 >1, (A24I1) 
v=] 
respectively, the equation determining z, namely 
y=M 
Gi=1=D a(xo, v)(x08)”/v™, (A251) 
v=] 
v=M 
or Gi=1= Yoa(xo, v)a’¥(xo2)"/v", (A25I1) 
v=] 


can only be satisfied by (xoz) <1, and the sum G; as well 
as all the other sums Gz and G;, etc., will be absolutely 
convergent; so that if the upper limit of summation, M, 
is sufficiently large, their values, and that of z will be 
independent of M. This case will be referred to as the 
case of convergent sums. 


y=M 


i a(xo, v)/vy™1<1, 


v=] 


If, however, (A261) 


v=M 
or _ a(Xo, va’ /yn- <1, 


v=] 


(A26II) 


which can only be satisfied in case I if m >2, and in case II 
if a<1, we find that the solution of (A25) will be xz=1+e 
where ¢ is an extremely small quantity. We may estimate 
the value of ¢ in the following manner, The order of mag- 
nitude of the last term, the term »=M in (A25), may be 
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expected to lie between 1/M? ard unity (actually it seems 
probable that it will be of order 1/M In M). Using this 
inequality, and inserting e* for xoz in the last term of 
(A25) we obtain 

1/M?<a(xo, M)e”*£/M""1 <1, 
or 1/M? <a(xo, M)a"te’*/M"" <1, 
resulting in 


(1/M)[(n—3) In M—In a]<e 
<(1/M)[(n—1) In M—Ina@], (A28I) 
(1/M)[(n—3) In M—In aJ]—(1/M!) Ina<e 
<(1/M)[(n—1) In M—In a@]—(1/M!) Ina. (A28IT) 
In case I we see that ¢ is of the order of magnitude 1/M, 
which will be sufficient for the purposes of our investiga- 
tion. In case II, however, we find that, to within terms of 
definitely smaller order of magnitude, we may write 


e=(—In a)/Mt=y7/M}, (A29I1) 


(A271) 
(A27I1) 


where y=—In a is positive because the case of divergent 
sums can occur only if a is less than unity. 

The condition of Eq. (A26) will always lead to the 
equation z=xX 9"! to within terms of order of magnitude 
M~- in I, and M- in II. This will be referred to as the case 
of divergent sums. 

If the sums are divergent in the sense discussed in the 
previous paragraph, the values of the individual terms of 
the sum G; will not decrease rapidly with increasing values 
of v. It is obvious that the sum Go, the individual terms of 
which have 1/v times the value of the terms of Gi, must 
be smaller than G; and independent of the value of M. 
The terms of Gz and G;, which are v, and v? fold greater, 
respectively, than the corresponding terms of G; will cer- 
tainly not converge. The values of the sums will be ex- 
pected to be dependent on the value of M. Various 
methods of estimation lead to the same qualitative result 
for the order of magnitude of the sums G;, namely 


Gj=Mr,  (j>1). (A30) 


The equation (A30) is nowhere of great importance to 
the argument of the article. It is of importance, but fairly 
obvious, that Gz and G; are extremely large and dependent 
on the limit of summation M, in the case of divergent 
sums. 

The case that 


M 
> a(xo, v)/y™1=1, (A311) 


v=] 
is used once in the article. In this case the solution of 
(A25I) is obviously that xo2 is exactly equal to unity. 
In this case the summation of G2 and G; offers no dil- 
ficulties. If »=5/2 we obtain 


y=M 
G:2> 1/v4=2M}, (n=5/2), 


v=] 


(A321) 


v=M 
G2> v=Mi, (n=5/2), 


v=] 


(A331) 
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whereas if n=3: 


Gi2>d 1/r2In M, (n=3), (A341) 
v= 
M 

G2>12M, (n=3). (A351) 


v=] 


The function Fo(Mx) 


The function F(M, J, x1, +++, Xp, --+) of the variable 
M, a very large integer, J, a parameter, and of the x,’s, 
is defined as 

F(M, J, x) = >> TI (Jx,)#0/uy!, (A36) 
My v 


Zvpy=M 


where, as before, x is taken to represent all the x,’s to- 
gether. If J is set equal to M we shall call 


F(M, M, x) = Fo(M, x). (A37) 

Direct differentiation of (A36) with respect to x, leads to 

10F(M, J, 1 o’F(M, J,x 

et Mn 5 ioe, CO 
J OX, Rai Ox,” 


which, upon setting J= M, gives the equation for Fo 
- A oF OF _ “(* -*) = a 1 3” Fo 
M Fo ax, ~ Ox, M Mw Fy Ox,” 


Since Fy is a homogeneous function of degree M in the 
variables y,=x,!/” it follows from Euler’s theorem that 


1 1 OFo 


0 (in F, 
2 aay, LL vty (=) =1. (A40) 


M Fy OVy v M 





(A39) 


We now propose to show that we can write approxi- 
mately*® 


In Fo(M, x)= M In f(x)+1n h(x)+In g(M), (A41) 


and that 


In f(x) =Go(x)—In 2(x), (A42) 
In h(x) = —[In G2(x)]/2, (A43) 
In g(M) = —[ln 27M ]/2. (A44) 


The approximation of Eq. (A42) holds as long as there is a 
real positive root to Eq. (A4) for z(x), and as long as 2(x) 
is independent of M as discussed in the section under 
convergence of the auxiliary functions. 
The combination of the assumption (A41) with (A40) 
shows that 
> »x,d In f(x)/ax, =1, (A45) 


> »x,d In h(x)/ax,=0, (A46) 


since f(x) and h(x) are assumed to be independent of M. 
We shall use the symbols 


2(x) =0 In f(x)/dx, (A47) 
21(x) =02(x)/dx1, (A48) 
a(x) =0 In h(x)/dx, (A49) 





* Professor M. Born_of Edinburgh has kindly communicated to the 
authors that this theof’m may be proved by using a Cauchy integral 
and the method of steepest descents. 
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and later identify 2(x) with the 2(x) defined by (A4). 
Upon using the assumption of (A42) in Eq. (A39), then 

developing according to powers of 1/M, and inserting the 

above symbols, one obtains 

a ( In Fo ) _ 9 In Se) 1 dlnhk(x) 1 1 dF 

ax,\ M ] ~~ ax, "M dx» — M” Fy dx,’ 


oo ae ba a 
7 ny \(R!)"* oxi* M 


n 














Dknp=v 


had = [ 2 “(< 2 ) aes ae | 4 (A50) 
=2 ~<a i? — ree, 
M s ae)” 22 


Equating powers of 1/M leads to 


0 In f(x) /dx,=2’, (A51) 
0 In h(x) /dx, = v2"(o/2—2;/227) + v22"s,/22?. (A52) 


The use of (A51) in (A45) gives us Eq. (A4) which 
identifies z(x) as defined by (A47) with 2(x) defined by 
(A4). Comparison of (A51) with (A10) shows that Eq. 
(A42) holds to an additive constant which we shall later 
show to be zero. Having identified z(x) we may now use 
(A7) to show that 2,/22?= —1/2G,2 and inserting this value 
in (A52), and (A52) in (A46) one obtains 


> *x,2"(6/2+1/2G2) — > v8x,2"(1/2G2) 


=G2(0/2+1/2G2)—G;/2G:=0, (A53) 
(G3 — G2) /2G,?. (A54) 
Using (A54) and the equation 2; in (A52) finally yields 

0 In h(x) /dx, = v2"G3/2G2 — v?2"(1/2G2). (A55) 


Comparison of this equation with (A8) shows that (A43) 
holds to within a constant of integration, which we shall 
now show can be taken as zero. 

We shall consider the special case that x1.>0 and 
x,=0, v1. In this case we may evaluate Eq. (A36) with 
J=M directly by the use of the Stirling approximation 
obtaining, since only the term wi = M occurs, 


In Fo(M, x)= M(In x, +1)—(In 27M)/2, 
(x, >0; x,=0, v¥1). (A56) 


or o/z= 


In this special case Go and G2 are both unity and z=1/x, 
so that the equations (A42) and (A44) inclusive lead to the 
same value for In F)(M, x).7 This suffices to show that the 
integration constants properly present in deriving (A42) 
and (A43) for In f(x) and In h(x) from (A51) and (A55) are 
zero and the In g(M) of (A44) respectively. 


Validity of the asymptotic expressions for F 


Combining equations (A41) to (A44) inclusive we ob- 
tain the final asymptotic expression for Fo(M, x) which is 


7 This special case x, >0, x» =0 if vk, may also be readily evalu- 
ated by both using (A37) directly and (A42) to (A45). In this case 
In s = —(In kxx)/k and Go=1/k, G2 =k. Evaluation by means of (A42) 
to (A45) leads to In Fo(M, x) =( Mp —In kx~) —(In 2xkM)/2. Direct 
evaluation of (A37) with J=M using the Stirling approximation for 
the factorial gives the above expression plus In & if M/k is an integer. 
If M/k is not an integer Fo(M, x) is zero. It is obvious that this extreme 
case is not covered by Eqs. (A42) to (A45), but it is of interest to note 
that if the Fo is averaged over a range of M's the above expression is 
obtained; since only every k’th M leads to Fo other than zero, and if 
this value for M/k an integer is divided by k we obtain agreement 
with the above. 
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correct for very large values of M, 


Fy(M, x) =e™ 0) /[2(x) ]” (24 MG2(x))3, 
In Fo(M, x) = M[Go(x) —In 2(x)] 
—(In 27MG.(x)]/2. (A57) 


If the function F, (M, x) is defined as 
F\(M, x)= F(M—1, M, x)=(1/M)dF)(M, x)/dx, (A58) 
then by direct differentiation 


F,(M, x) =e” G0) /[2(x) ]¥-1(24 MG2(x))}, 
In F\(M, x) = MGo(x) —(M—1) In 2(x) 
—[In 27MG.(x)]/2. (A59) 


The function Fo(M, x) is the same as Qr/N! if the 
x,’s are the vb;’s of the article and M is N, the total number 
of particles. The function F;(M, x) is the same as /?0, if J re- 
places M and the x,’s are the irreducible integrals 6;. 

The expressions (A57) for /o(M, x) and (A59) for 
F,(M, x) are correct asymptotic expressions for these 
functions as long as they can be approximated in the gen- 
eral form of (A41). We have already discussed the case that 
one of the x,’s alone is different from zero. Unless » is unity, 
the approximation (A41) is invalid. Even then (A57) and 
(A59) give correct average values for a range of M values. 
It appears to be sufficient for the validity of (A57) and 
(A59) that two x,’s alone be different from zero, and that 
the ratio of the smaller to the greater. x, be larger than 
v/Mv’. 

One other case that (A57) and (A59) as well as the ex- 
pression (A41) are invalid is that Go(x) <0. If all the x,’s 
are positive the equation Gi—1=0 has only one positive 
root for z. If the x,’s are all functions of one variable, say 
the temperature 7, and are all positive at low values of 7, 
it may be that some of them become increasingly negative 
with increasing 7’s. Because Gi—1 is a finite polynomial 
in z with its zeroth coefficient minus unity, the single 
positive root of the equation at low T’s cannot become zero 
or infinity, but may, and does, disappear into the complex 
plane with increasing values of 7. At this temperature 
T=T., z must be a double root, and 


0(G,—1)/d In z=G2.=0. (A60) 


At T=T, the real root z disappears, G2 becomes zero, and 
the expressions (A41), (A57) and (A59) become invalid. 
This temperature 7, is shown in the article to be the critical 
temperature. 


The coefficients of I18,”” in b; 


It can be shown that there is a one to one correspondence 
between the coefficient of the term characterized by the 
n,'’s, the powers of the v’th normalized irreducible integral 
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in a term of b;, (> »m,=/—1), and the answer to the prob- 
lem discussed in the following paragraph. 

Suppose we have , identical mechanical frames of 
type 1, m2 of type 2, ---, m, of type v. The frame of type v 
has v+1 holes into which bolts can be inserted, the frames 
being so built that every one of the »+1 holes is exactly 
equivalently situated with respect to every other one of 
the holes (a geometrical impossibility in three dimensions 
for most large values of v). Suppose we have / numbered 
bolts with /= }-»n,+1. The problem is to determine how 

v 


many ways we can insert the / bolts into the /—1+ on, 
holes so that the frames are all singly connected to each 
other by common bolts. 

Each bolted arrangement can be broken down into / 
different dissociated arrangements consisting of frames 
and one free bolt, where each frame has » bolts in it and one 
free hole, and no frames are connected by a common bolt. 
We can dissociate in exactly / ways because we can choose 
from the bolted arrangement any one of the bolts as the 
free bolt, and the dissociation process is then uniquely 
determined. The number of ways that we can make a 
dissociated arrangement is 


1!/TIn, '(v!)"*. 


If D is the number of bolted arrangements we can make 
make from one dissociated arrangement then the desired 
answer is DI!/IIn,!(v!)"1. 

v 


For the small values of n= Yn, we know by counting 
that D=/""!, We now prove that if the equation D=/""! 
holds for ” it also holds for n+1. To do this we add a 
frame v to a system of ” frames and bolts and show that the 
new D becomes (/+ v)", (method of complete induction). 
We consider that the old frames are first bolted together 
into v+1 (or fewer) groups, knowing how many ways that 
can be done, and then these groups are bolted to the new 
frame. The resulting complicated sum can be shown to be 
equal to (/+yv)”". Use must be made of the fact that 


v=m 


p> [*) (—1)"(a+»)?=0 if P<m. 


v=0 
The final answer then becomes /!/"/IIn,!(v!)"¥l?, which 


is the same as Eq. (19) if it is remembered that v! is in- 
cluded in the definition of 8, as a normalization factor. 
this proof is primarily due to Dr. M. Goeppert Mayer to 
whom the authors are greatly indebted. The original 
equation was deduced by Dr. P. G. Ackermann from a 
numerical calculation of the first several members. With- 
out this hint as to the correct solution, the proof would 
have been impossible. 
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It is shown that the theoretical equations of the preceding paper lead to the prediction that, 
in a finite density and temperature range above the temperature of disappearance of the 
meniscus, a range of microscopic densities are equally stable, and in a sealed bomb different 
densities may be observed at a given temperature and average density. The densities observed 
in various parts of the bomb will depend in a predictable way on its previous history. This is in 
agreement with the published results of O. Maass on systems in the neighborhood of the critical 


point. 





N the previous paper the authors have deduced 

certain equations for the thermodynamic 
properties of a system of chemically saturated 
molecules. Their deductions were based on 
general equations and did not involve any special 
assumptions as to the behavior of the system. 

In the neighborhood of the critical point the 
equations lead to a P—V diagram for the system 
qualitatively similar to that of Fig. 1. In this 
figure pressure is plotted as a function of volume, 
at various temperatures: the heavy solid lines. 
The curves at two characteristic temperatures 
T,, and 7, are especially marked. The shaded 
area enveloped by the dashed line is the region in 
which the lines are horizontal, that is where 
(@P/0V)7=0. The right-hand boundary of this 
region is called v., the volume of the saturated 
vapor, the left-hand boundary 2; is the volume of 
the condensed phase. The dotted lines inside the 
shaded area represent the extrapolated pressure 
of the supersaturated vapor. The dot-dash line on 
the left-hand side of the shaded area at JT» 
represents a discontinuous transition in the 
properties of the condensed phase.* 

Below T,, the isothermals show discontinuous 
changes in slope at v, and v;. At 7 the surface 
tension of the condensed phase disappears, and 
above T,, there will be no visible sharp meniscus 
between condensed and dispersed phase. The 
discontinuous transition in the condensed phase 
shows no heat of transition, at least at the vapor 
pressure P,,. In general the transition tempera- 
ture will be expected to vary with volume, and it 
is shown here increasing with decreasing volume 
(the dot-dashed line lies above the T,,, isothermal), 





* This transition is probably present. It is not a necessary 
consequence of the equations of the previous article. 


but of this there is no certainty. The isothermals 
go continuously with continuous change of slope 
through the shaded region between T,, and T,. 
At T, the isothermal has a horizontal tangent at 
only the one point v,, P.. Above T, the isothermal 
is never horizontal. 

The part of the diagram to the right of V,, 
and all of it above T,, is deduced strictly from the 
equations of the previous article. The portion to 
the left of V. and below 7, represents a logical 
guess from the nature of the kinetic picture 
deduced from these equations. The authors, 
however, are less certain of the details of this 
portion of the diagram. 

Below T,, and above 7, the diagram corre- 
sponds to the usual conception of the diagram 
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below and above the critical point respectively. 
Between 7, and 7, a system in the shaded area 
will be expected to have certain unusual 
properties. 

Below T,, compression of the system through 
the volume region between v, and vy; will be 
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accompanied by the separation of two distinct 
phases of density 1/v,; and 1/v,. The microscopic 
density corresponding to a volume between v, 
and v; is unstable. This is indicated on the 
diagram by the dotted lines showing the higher 
pressure of the supersaturated vapor. As is well 
known this phenomenon is associated with the 
fact that the surface tension increases the vapor 
pressure of very small droplets of liquid. 

Between 7, and T.., however, the denser phase 
has no surface tension, and corresponding to this, 
the isothermals pass smoothly through the 
condensing region, allowing no extrapolation to a 
higher pressure for the supersaturated vapor. 
Indeed all microscopic densities in the shaded 
region between T,, and 7, will correspond to the 
same pressure (at the same temperature), and to 
the same Gibb’s free energy. Compression of 
the system at constant temperature 7 with 
Tm<T<T, will be accompanied by a uniform 
increase in density throughout the whole system, 
the pressure remaining unchanged between 2, 
and vy. 

It is true that after infinite time, owing to the 
gravitational gradient in a real system, and the 
consequent higher pressure at the bottom than at 
the top, the system should attain an equilibrium 
state with the density higher at the bottom than 
at the top of the tube. However, the gravitational 
pressure gradient amounts to about one part in 
10° per cm of height, and it would therefore not 
be surprising if uniform density could be main- 
tained for weeks or months. 

Moreover, should a density difference once be 
obtained by any means, with higher density at 
the bottom of the system than at the top, while 
the temperature was between T,, and 7., and 
both densities were well within the shaded area 
between v; and v,, then uniform compression of 
the system should result in a uniform logarithmic 
change in the densities of both portions without 
a change in their relative abundance. Similarly a 
change in temperature should result in no change 
of the two densities, but be accompanied by an 
increase in the pressure of the system. 

If, however, by means of a volume or tempera- 
ture change, the density in one portion of the 
tube reached a value corresponding to v, or vs on 
the envelope of the shaded region, then that 
portion should resist, by increase in relative 
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abundance, any effort to change its density to 
one outside the shaded region of the diagram. 
Cooling of such a system should be expected to 
result in the formation of the meniscus at some 
temperature lower than 7;,, a typical super- 
cooling effect. The extent of the supercooling 
would be expected to be greater the further the 
densities in various parts of the system were 
from the values corresponding to the volume v,,, 
and vsm at 7. A system with uniform density 
1/v., approximately midway between 1/vsm and 
1/vsm, would be expected to show the greatest 
supercooling. 

Experimental evidence for the correctness of 
these concepts is good. Various experimenters 
have observed discontinuities in different prop- 
erties of the condensed phase in the neighborhood 
of 7;,, the temperature of disappearance of the 
meniscus, which is that listed in all existing tables 
as the critical temperature. These results are 
summarized in an article by O. Maass.! Solubility, 
dielectric constant and adsorption values have 
all been found to change discontinuously or 
rapidly at this temperature. Perhaps the most 
convincing change is that in the reactivity of HCI 
with propylene, which shows the usual increase in 
reaction rate with temperature as long as the 
HCI is liquid, but the rate drops suddenly to 
zero as the critical temperature (T7;,) of HCI is 
exceeded even though the density of the gaseous 
HCI is made comparable to that obtaining 
previously in the liquid. This has been investi- 
gated extensively by Maass.’ 

The evidence that the densities of the “‘liquid”’ 
and “gas” are not identical at temperatures 
considerably above the temperature of the 
disappearance of the meniscus is extensive and 
convincing. Again the paper of Maass previously 
referred to! contains a good historical summary. 
Apparently all investigators who have made 
measurements just above the critical tempera- 
ture have found differences of density in the 
material remaining above the meniscus level and 
that below. Any contrary evidence would, in 
addition, have to be regarded extremely critically 


in view of the experimental results of Maass and 


( oy Steacie and Maass, Can. J. Research 9, 217-239 
1 ’ 

2 Sutherland and Maass, Can. J. Research 5, 48-63 
(1931); Holder and Maass, ibid., 15B, 345-351 (1937). 
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STATISTICAL MECHANICS 





the conclusions of the authors that a uniform 
density, if once attained, would persist almost 
indefinitely. 

Maass! has drawn a simple conclusion from a 
common observation. It is found that the 
disappearance of the meniscus in a closed tube 
can be observed with varying total quantities of 
material in the tube. This can only be true if the 
density above and below the meniscus is not 
equal at the time of its disappearance. 

Callendar*® has stated that the dividing line 
between the material in bottom and top of a 
quartz tube containing pure water can be seen 
above 374°C, the temperature at which the 
meniscus disappears. He reports a difference in 
density in the ratio 5 to 3 at this temperature for 
water, and indicates that the true critical point is 
somewhat over 6°C higher. 

The best and most complete experimental 
evidence is contained in the excellent work of 
Maass!: *: § and several co-workers. Densities were 
measured to about 0.2 percent by means of float 
suspended on a quartz spiral which could be 
moved to different parts of the system. Tempera- 
ture control and the purification of materials 
were attended to with the greatest care. The 
density difference observed above the tempera- 
ture 7, was as much as 100 times the experi- 
mental error. The work was done on methyl 
ether, propylene, and ethylene. 

The experimental results which appear to us 
to be most significant, perhaps under the influence 
of the fact that they are in perfect agreement 
with our theory, were as follows. Upon heating a 
closed bomb containing definite quantities of 
material through 7,,, the temperature at which 
the meniscus was seen to disappear, the densities 
in top and bottom of the tube followed curves 
corresponding roughly to the envelope v, and v; 
of the shaded area in our figure. After reaching a 
temperature between 7, and 7, and subsequent 
cooling, the density in top and bottom remained 
constant until a temperature somewhat below 
T,, was reached. At this temperature the 
meniscus reappeared and the density difference 





* Proc. Roy. Soc. A120, 460-472 (1928); see also Engi- 
neering 126, 594-5, 625-7, 671-3 (1928). 


— and Maass, Can. J. Research 9, 613-629 
. Maass and Geddes, Phil. Trans. Roy. Soc., London 
A236, 303-332 (1937). 





OF CONDENSING SYSTEMS 








103 


increased rapidly with falling temperature. 
Density differences persisted for as long as 6 
hours without any tendency to decrease. The 
density difference was unaffected by mechanical 
stirring. This appears to be reasonable, since 
ordinary stirring, although effective in raising 
macroscopic portions of the high density material 
into the top of the tube, would not be expected to 
break up the material into the submicroscopic 
size droplets which would be prevented from 
settling by Brownian motion. If the density was 
once equalized by raising the temperature above 
T., or by cooling the top of the tube 1/10°C, a 
uniform density persisted unchanged. With 
uniform densities the supercooling phenomenon 
was most marked. The pressure exerted by the 
system at any temperature above 7, was the 
same whether the system was cooling or the 
temperature was rising. The behavior of any one 
bomb, with a given filling, was reproducible day 
after day. 

The following observations appear to be at 
variance with the authors’ theory: 

The behavior of the different bombs was not in 
complete numerical agreement. This may possi- 
bly be due to a conceivable failure of the float 
spiral balance to record absolute densities due to 
an unpredictable reversible variation in con- 
stants with pressure and temperature. (The 
instruments were apparently checked after the 
operations.) More likely, this lack of agreement 
between different bombs may be due to the effect 
of varying amounts of some impurity present in 
small amounts. For various reasons the influence 
of a trace of impurity might be expected to be 
considerable. However, since the total amount of 
impurity can scarcely exceed 0.1 percent, we 
believe that the effect is not caused by impurities 
alone. 

The most extreme case of varying behavior of 
different bombs was that of bomb No. 13'. The 
over-all volume of this bomb corresponded to one 
exactly midway between vs» and Vsm so that the 
meniscus remained stationary before disappear- 
ance. In this case no density differences were 
observed, but the supercooling was marked. The 
authors failed to reproduce this later with a 
different set-up.‘ This is certainly unaccountable, 
although one may assume that the stationary 
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meniscus was peculiarly favorable to super- 
heating. 

The curves of density at the top and bottom of 
the vessel plotted against temperature, on 
heating, do not show the sharp break in slope 
predicted by us at T7,,. Since we cannot, as yet, 
estimate the magnitude of this change it may 
not exceed the experimental error. The effect 
may also be masked in some way by the influence 
of impurity traces or by other minor effects not 
considered here. The experimental curves do 
show a marked curvature at 7;,. 

Certain observations on the effect of varying 
the total volume of the system‘ are not in 
numerical agreement with the v,, v;, T, curves 
deduced from the behavior of other bombs. This 
may, however, be linked with the fact that 
individual bombs are not in exact numerical 
agreement with each other. 

The P-—V isothermals of Maass do not show the 
horizontal portion demanded by us. The experi- 
mental points are scarcely sufficiently close 
together to decide this, and here, certainly, the 
effect of impurities can be predicted to be con- 
siderable, since total and not partial pressures 
are measured. It is inconceivable that pressure 
differences can persist in different parts of the 
same bomb, and this, coupled with the observa- 
tion that varying density differences (obtained 
on heating and cooling) lead to the same pressure, 
appears to us to require the assumption of a 
horizontal region in the P-V diagram. 

In short, the authors believe that the obser- 
vations of Maass and co-workers lead directly to 
conclusions essentially the same as those deduced 
from the equations of the preceding paper. The 
observations, however, show slight deviations 
from the predictions. Whether these can be 
accounted for by the influence of impurity traces, 
or by some unexpected failure of the experi- 
mental technique, or whether they indicate 
additional complications in the theory cannot be 
easily decided at the present time. 


S. F. HARRISON AND J. E. 
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The paper of Verschaffelt® who shows that 
variations in the density in top and bottom of 
the tube may be due to impurities, and who is 
quoted as stating that differences as large as 
40 percent may be caused by 0.1 mole percent of 
impurity, must be considered. Verschaffelt shows 
that, from the known influence of Os on the 
critical constants of COe, we can conclude that if 
0.1 percent Oz is in one “‘phase”’ and none in the 
other, then just above the critical point the 
density difference of the two regions would be 
36 percent. However, if we assume 0.1 mole 
percent Os as impurity in the total COs, and 10 
percent difference in concentration above and 
below the meniscus, we arrive (using Verschaffelt’s 
figures) at about 6 percent density difference just 
above the critical point, and only about 1 
percent density difference at 1°C above the 
critical point. Maass’ material, especially in the 
case of dimethyl ether, can scarcely have con- 
tained as much as 0.1 percent impurity, and 
since the impurity must have about the same 
vapor pressure as the substrate it seems reason- 
able to assume a concentration difference in 
liquid and gas of approximately 10 percent. We 
must compare Verschaffelt’s 1 percent predicted 
density difference with a difference of 15 percent 
observed by Maass 1° above 7;,, and of 10 
percent some 6°C above T,,, where Verschaffelt’s 
predicted difference would be negligible. Al- 
though it seems that Verschaffelt showed ade- 
quately that work before 1904 could be explained 
by the influence of impurities, the work of 
Callendar and especially the careful experiments 
of Maass indicate that a more fundamental 
phenomenon is at hand than can be due to 
impurities alone. Nevertheless, the large effect of 
impurity traces is regarded by us as lending 
weight to our supposition that some of the 
numerical discrepancies of Maass’ different 
bombs are to be explained in this way. 


6 J. E. Verschaffelt, Comm. Phys. Lab., Leiden, Supp. 


Nr. 10 (1904). 
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Note on the Second or Gaussian Approximation in the Heisenberg Theory of Ferro- 
magnetism When S>3 


J. H. Van VLEcK 
Harvard University, Cambridge, Massachusetts 


(Received January 12, 1937) 


N the Heisenberg theory, it is customary, ex- 

cept at very low temperatures, to write the 
expression for the intensity of magnetization in 
the form 


M=NgbBi(g8pH+apJ+bp°s?+---), (1) 


where J is the exchange integral, 8 is the Bohr 
magneton, /7 is the applied field, and 


+8 
p=1/kT, B,(0)= > mre-™/de-™*. (2) 


m=—S 


Eq. (1) shows that (aJ+bpJ*+ ---)/g6 plays the 
role of an internal or ‘“‘molecular’’ field. The coeffi- 
cients a, b, --- are functions of MM. It is quite 
easy to prove, as has often been done, that 


a=22M/N¢gp, (3) 


where z is the number of neighbors. The evalua- 
tion of 6 involves the second or so-called Gaussian 
approximation, since it can be shown that (1) 
would be exact, without the unwritten higher 
order terms, were the energy values for a given 
spin distributed in accordance with a Gaussian 
error law. If S= 34, then a calculation made either 
by means of group theory! or the vector model? 
shows that 


b= —(2zM/NgB)+(42M*/N%g%6%), (4) 


The proper formula for 6 has‘not previously been 
given when S>}3, as neither method can be used 
for this purpose,’ at least in a simple way, except 
when S=}. 

A quite different approach has, however, re- 
cently been developed by Opechowski.‘ He de- 
velops the partition function as a power series in 
J. The constant b in (1) is uniquely fixed by the 
requirement that the expansion of (1) agree with 





'W. Heisenberg, Zeits. f. Physik 49, 619 (1928). 
_° J. H. Van Vleck, The Theory of Electric and Magnetic 
Susceptibilities, Chap. XII. 

* See especially footnote 16, p. 329 of reference 2. 

*\W. Opechowski, Physica 4, 181 (1937). 
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this series development to J? inclusive. Ope- 
chowski shows that in this way formula (4) can 
be derived quite easily when S=}. His work re- 
veals perhaps more clearly than does the ordinary 
calculation by means of an assumed Gaussian 
spread,> that the practical convergence of the 
form (1) is very poor at ordinary temperatures. 
Indeed, better agreement with experiment is 
usually obtained if one stops with apJ than if 
one includes bp?J*, or, as he shows, even dp*J*. 

It is the purpose of the present note to show 
that the calculation of the partition function 
which I gave in a previous paper on ferromag- 
netic anisotropy® permits extension of the results 
to arbitrary S. Also it is well to mention that the 
general mathematical procedure which I used is 
in many respects similar to that employed by 
Opechowski, although obtained independently. 
The essence of both schemes was the development 
of the partition function in a system of represen- 
tation in which each atom is separately space 
quantized in a field H. Because of the invariance 
of the diagonal sum, use of such a system ‘is al- 
lowable even with interatomic perturbations. 
Saturation effects are not overloooked, and so 7 
can be arbitrarily large. In one notable respect, 
Opechowski's calculations are more general than 
ours, as he carries the development far enough to 
include third and fourth powers of the exchange 
integral, corresponding to unwritten terms 
cp®J*+dp*J* in (1). On the other hand, my com- 
putation, though carried less far, was made for 
arbitrary S, and allowed the interaction between 
atoms to be of the dipole-dipole rather than ex- 
change type if desired. The introduction of dipole- 
dipole coupling is obviously not needed in con- 
nection with the present note, but was wanted 
for other purposes, and somewhat complicates 
the mathematics, as it makes the perturbing 
potential cease to commute with the main or 
Zeeman energy. 


5 Cf., however, D. Inglis, Phys. Rev. 42, 442 (1932). 
6 J. H. Van Vleck, Phys. Rev. 52, 1190 (1937). 
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For arbitrary S, to terms in J, the formula for 
the magnetic moment is 


M/NgB=B,+22pJB,(B2,—B,’) 
+2p°J?(B,—B,’)(—82B; 
+42B,B.+4B,*—4B,B2—B)) 
+2p°J?(B;—B,B:)(22B,?—2B?+B.—S?—S). 
: (5) 


Eq. (5) is derived from my previous expression 
for the partition function in almost exactly the 
same way as is Eq. (3) of my article on adiabatic 
demagnetization appearing in the present issue,’ 
to which the reader is referred for further detail. 
The only differences are that in (5) all dipole- 
dipole terms are omitted, and S is not specialized 
to the value S=}3. In (5) the argument of the B 
functions is merely® g8pH. On the other hand, in 
(1) it is g8pH+apJ+bp’J*---, or to our degree 
of approximation 


gbpH + 22pJB,+422p°J°B (B,—By)+b7?, 


as is seen from (3) and (5). By expanding (1), 
in essentially the same way as explained more 
fully on page 192 of Opechowski’s paper,‘ it is 
found that (1) and (5) are the same to J? provided 


b=2(4B,'—4B,B.—B,) 
+2(—2B,/°+B,—S?—S) 
X (B3;—B2B,)(B2—B,*)". (6) 


In Eq. (6), B: may be replaced by M/N¢g@, as in 
the second-order coefficient 0} it is unnecessary to 
distinguish between arguments exclusive and in- 
clusive of the inner field. 

In the special case S=}, one has B3;=}B,, 


7 J. H. Van Vleck, J. Chem. Phys. 6, 81 (1938). 

8Jn the application to anisotropy in reference 6, an 
inner field was included in the argument of the B functions, 
but is to be omitted in deriving Eq. (5) from our previous 
calculations, as (5) is obtained by treating the exchange 
energy as part of the perturbing potential. 


VAN VLECK 


B,=} and (6) reduces to the simple expres- 
sion (4). 
If S=1, the following formulas apply 


B;=B,, B,=2 sinh 0/(1+2 cosh @), 
Bz=2 cosh 6/(1+2 cosh @). (7) 


By elimination of @ between the last two relations 
of (7), a quadratic equation is obtained giving 
Bz in terms of B;= M/ Ng. It is thus possible to 
reduce (6) to an explicit function of M alone, 
but the resulting formula for 6 involves I/ 
through radicals and is too cumbersome to be of 
much value. 

When S>1, the interrelations between B,, Bz, 
B; yielded by (2) are more complicated than (7). 
It is still possible, in principle, to express Bo, B; 
in terms of B or M, but only by solving algebraic 
equations of high degree, viz. 2S. The practical 
use of (6) thus appears to be limited to S=3 and 
possibly S=1, unless the field is weak, as is the 
case above the Curie point. The limiting general 
form of (6) for weak fields is tractible, and has 
been given in a previous paper, as well as dis- 
cussed in relation to the critical conditions for 
ferromagnetism.’ The rather unpleasant result 
(6) of acalculation of 6 for S> 3 in arbitrary fields 
is of interest mainly in connection with the light 
which it throws on the failure of group theory or 
the simple vector model to furnish a method of 
computing 6 when S>}3. Either of these proce- 
dures appear capable of yielding only an answer 
such as (4) which is a polynomial in M (or at 
least this is what analogy to the way they work for 
S=4 would suggest). Hence they must fail when 
S>4, since then the true formula for 6 is far 
from a polynomial. 


9 J. H. Van Vleck, J. Chem. Phys. 5, 320 (1937), especially 
footnote 20. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


Bond Force Constants and Vibrations of Tetramethyl- 
methane and Tetramethylsilicon 


In order to get information regarding the bond force 
constants ‘of the molecules tetramethylmethane and 
tetramethylsilicon, we have treated some of the vibrations 
of these molecules by the normal coordinate method. These 
molecules possess a large number of frequencies, but only 
six for each are considered here. First the totally symmetric 
nondegenerate frequencies are rigorously treated and 
second, certain vibrations are calculated assuming that 
each methyl group behaves like a single atom. 

The totally symmetric nondegenerate frequencies of 
C(CH3)4 and Si(CHs3)4 can be conveniently treated by 
considering the symmetric vibrations of a single methyl 
group against a carbon or a silicon atom of infinite mass. 
This simplifying device can be used because the central 
atoms are not involved in the motion of those particular 
vibrations, except insofar as they move with the centers 
of gravity of the molecules. Assuming valence forces, the 
secular equations for the totally symmetric nondegenerate 
frequencies of C(CH3)4 and Si(CHs3),4 are the same as the 
equation for the symmetric (Raman active) nondegenerate 


TABLE J. Tabulation of results. 


























CALCULATED FREQUENCIES 
OBSERVED 
FREQUEN- 
CIES* 4 A B 
Tetramethylmethane 
331 cm7 (331) em (2) 
415 , 383 (3) 
732 (732) cm=! (732) (1) 
1248 1264 (3) 
1452 (1452) 
2892 2894 
Tetramethylsilicon 
202 (202) (2) 
239 220 (3) 
598 (598) (598) (1) 
863 793 (3) 
1427 1427 
2905 2893 
kco_y =4.76X 10° ko_-rp=4.71X 10° 
dynes/cm dynes/cm 
ke =0.468 koccra) =0.321 
ko_-c=5.14 ksi_-r=3.14 
ksi-c=3.31 ko si—rs) = 9.120 
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not later than the 15th of the month preceding that of the 
issue tn which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


frequencies of ethane,': ? except that the Hooke’s law force 
constant for the carbon-carbon single bond of ethane is 
replaced by one-half the Hooke’s law constant for the 
carbon-central atom bond of the tetramethyl molecule. 
The methods employed in arriving at the equation, except 
as noted, are the same as those described in an earlier paper." 

Assuming the methyl groups behave like single atoms 
of atomic weight 15, four vibrations for each of the mole- 
cules CR, and SiR, can be calculated. The equations used 
are the valence force type equations applicable to molecules 
of the methane type, XY,4. Of the four ‘‘methane-like”’ 
vibrations, one is totally symmetric and nondegenerate 
and accordingly was also considered by the first method of 
calculation described above. 

The results of the numerical computations are given in 
Table I. Under column A are given frequencies calculated 
by the first described (rigorous) method and under column 
B, those computed by the second (approximate) method. 
The figures in parentheses following the numbers in 
column B are the degeneracies of the vibrations. Beneath 
the columns are given the values of the force constants 
used. The constant kg refers to bending of angles in the 
methyl group, and the constants kg(CR,4) and kg(SiR,) 
refer to bending of the angles in the molecules indicated 
for which the central atom forms a vertex. The constant 
kc_u was taken from reference 1, while the other constants 
were determined by using those frequencies enclosed in 
parentheses. 

The general agreement between calculated and observed 
results is good enough to justify the assignment of fre- 
quencies. Rank and Bordner,? by means of analogies, made 
different assignments for the higher triply degenerate 
frequencies here considered. However, our calculations seem 
to point fairly conclusively to the assignments in the table. 

It will be observed that the values for kc_c and kgi_c 
are, respectively, greater than those for kc_p and kgj_p. 
This would be expected because of the failure to take into 
account interactions within the methyl groups in obtaining 
the values for the latter. This result is in accordance with 
recent observations in connection with ethane and related 


molecules.® 
Department of Chemistry, 
University of Illinois, 
Urbana, Illinois, 
January 7, 1938. 
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3 Rank and Bordner, J. Chem. Phys. 3, 248 tio3s). 

4 Kohlrausch and Képpl, Zeits. f. physik. Chemie B26, 209 (1934); 
D. H. Rank, J. Chem. Phys. 1, pt ve Kohlrausch and Barnés, 
Anales Dag m Y Quimica 30, 733 (1932). 

5 F. T. Wall, J. Am. Chem. Soc. 60, 71 (1938). 
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Alterations in the Infrared Absorption Spectrum of 
Water in Gelatin 


The absorption spectra in the region 1—2.5u4 have been 
recorded in the four following instances with a self- 
registering instrument of slit widths equivalent to 0.003: 
(1) a 0.20 mm specimen of gelatin dried to constant weight 
at 120°C; (2) a portion of the same specimen after it had 
absorbed an amount of water equal to 35.7 percent of the 
total weight and swollen to a thickness of 0.28 mm; 
(3) a 0.09 mm cell of water at room temperature; (4) a 
0.17 mm cell of water. The spectrum of (1), in spite of 
the complexity of the protein molecules involved, can 
readily be interpreted in terms of CH and NH vibrations; 
any absorption by OH groups or by possible residual 
water is not observable. Although there may be slight 
changes in the NH bands, the main differences in the per- 
centage transmission curves obtained from the reduction 
of (1) and (2) and of (3) and (4) can fairly certainly be 
ascribed to alterations in the water spectrum. Since equal 
effective thicknesses (0.08 mm) of water were involved, 
these alterations must arise from differences in the struc- 
ture or in the environment of the water. The following 
modifications are observed, and provisional causes are 
assigned. 

(a) The 1.444 band (v;+v,). The maximum absorption 
is slightly shifted to a longer wave-length in the gelatin 
specimen, apparently owing to the partial or complete 
removal of the short wave portion of the band usually 
associated with the water vapor molecule. This indicates 
that most or all of the water molecules are bound mole- 
cules. Absorption at the maximum is increased but total 
absorptions in the two instances are comparable owing to 
a diminution of continuous absorption in the 1.65 region. 

(b) The 1.964 band (v,z-+v53). The location of the 
maximum is practically unaltered. This need not be con- 
strued as indicating that unbound molecules are present. 
For, since this band is a combination of the fundamental 
deformation frequency (6.14) and the fundamental asym- 
metric valence frequency (2.74), it should show the 
characteristics of both fundamentals; and the temperature 
effects and the hydration studies of Ganz,! as well as vapor- 
liquid band shifts, clearly show opposite behaviors in 
these two bands. A statement analogous to the last sentence 
of (a) can also be made in connection with this band. 

(c) The 1.79 (vz+vs+vp). This band, whose frequency 
differs from that of 1.964 by an amount equal to the 
frequency of “hindered rotation,” ve,? is stronger and 
sharper than in pure water. This seems to indicate, first, 
that at least some of the water molecules are bound in 
such a way as to permit hindered rotation about the axis 
of least inertia and, second, that the average field of 
potential energy is equivalent to that in pure water but 
that it is more homogeneous. 

(d) A new 1.354 band (vz+ve+vpr). This band is 
analogous to the 1.79u band and is doubtless observable in 
the gelatin specimen but not in pure water because of a 
greater homogeneity in the field. 

(e) The 2.4u region. The absorption in this broad region, 
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more intense in the gelatin instance, is probably to be 
interpreted in terms of vz+vp. 

(f) The absorption bands at 4.7 (vgs+ve) and 20 (vp) 
should, if observed in gelatin, be sharper than in pure 
water. 

It is the authors’ opinion that the absorption at 3yu 
found by Buswell, Krebs and Rodebush® in a gelatin 
specimen dried for 14 hours at 122°C is caused by protein 
NH groups rather than by residual water. 

JosePH W. ELLIs 
JEAN BATH 
University of California, 
Los Angeles, California, 
January 13, 1938. 


1 Ganz, Ann. d. Physik. 28, 445 (1937). 

2 Ellis, Phys. Rev. 38, 693 (1931); Cartwright, Phys. Rev. 49, 470 
(1936). . 

3 Buswell, Krebs and Rodebush, J. Am. Chem. Soc. 59, 2603 (1937). 





Errata: Thermodynamic Functions of the Chloro- and 
Bromomethanes, Formaldehyde and Phosgene 


(J. Chem. Phys. 6, 25 (1938)) 


Because of a systematic error the values of —(F°— E,°)/T 


in Tables VI and VIII are too large by 1.987 cal./°K mole. 
The value of AE,’ for the reaction 


Cgg+H2+1/202= H.CO, 


as calculated from the work of Newton and Dodge should 
be —26.0 kcal./mole instead of —27.0 kcal./mole. This 
leads to a final value of AE, for this reaction of —26.5 
kcal./mole in place of —27.2 kcal./mole. 

For the free energy of formation of the liquid chloro- 
methanes 


Cao+H2+Cl2=CH:Cl(1), A Fans? = — 14.7, 
Cag +3H2+3Cla=CHCIs(1), AFae° = —15.9, 
Cggt+2Cl2=CCIi(2), A Fg.° = — 16.6. 


The free energies of formation of the gases are given in 
Table X. 


TABLE X. AF® of ideal gases. 



































= H2CO COCl CH:;Cl CHeCle CHCl; CCl 
298.1} —25.2 | —49.7 | —130 | —14.3 | —15.1 | —15.5 
350 —248 | —49.1 =§18 | —124 | —137 | —133 
406 =244 | —484 | —106 | —113 | ~—12.3 | —12.1 
500 —23.5 | —47.2 —8.0 —9.2 —9.5 —8.9 
600 —22.5 | —459 —5.4 -6.6 —6.7 $5 
700 =—215 | —64.7 —26 —4.0 -3.8 
800 =—5 | —&3A +0.1 —1,3 -0.7 
900 —194 | —420 +3.0 
1000 —18.3 +5.8 
1100 17.2 8.7 
1200 —16.1 11.5 
1300 —15.0 
1400 —13.9 
1500 —128 
D. P. STEVENSON 
J. Y. BEAcH 


Princeton University, 
Princeton, N. J., 
January 14, 1938. 








posi 
dete 
inte 
clus: 
mad 
are | 
whic 
Ir 
in a 
150- 
duce 
on t 
ther 
amo 
the 
cont 
there 


wate 
In F 
first 
is ca 
Oxyg 
Start 
addit 
perin 

Sir 
sodiu 
conti 
soluti 
synth 

Th 
fill th 


10 


iv 


4 





mm 


wns 





LETTERS 





Photosynthesis* 


Photosynthesis has long been known to be composed of 
two parts, one photochemical and temperature independent 
and the other thermal. The average life of the thermal 
reaction has been determined by measurements made in 
flashing light.} 2 

Several explanations of the nature of the thermal (or 
Blackman) reaction have been advanced, some postulating 
the return of the plant chlorophyll to its original state 
after having been involved in the reduction of carbon 
dioxide. Franck and Herzfeld* describe the mechanism of 
the reaction as the decomposition of per-acids and per- 
aldehydes to formic and carbonic acids in chain processes 
started by light. However, the results of Emerson and 
Arnold! show the Blackman reaction to be inhibited by 
hydrogen cyanide, a characteristic of certain enzyme 
processes. 

The mechanism of photosynthesis here presented is one 
in which the initial processes result in the formation of an 
organic compound of the general formula CnHoNOwy and a 
peroxide—probably hydrogen peroxide. The formation of 
hydrogen peroxide is not eliminated by Franck and Herz- 
feld’s* calculation of the energy involved in photosynthesis. 
The Blackman reaction is then the enzymatic decom- 
position of the peroxide. This latter would be the rate 
determining process in oxygen evolution at high light 
intensities and high carbon dioxide concentrations. Con- 
clusions about photosynthesis drawn from measurements 
made under such conditions are not reliable. Our opinions 
are the result of experiments, performed in this laboratory, 
which are described below. 

Into a suspension of cells (chlorella vulgaris), contained 
in a Warburg manometer apparatus illuminated by a 
150-watt tungsten lamp at a distance of 5 cm, was intro- 
duced a small amount of catalase solution and the effect 
on the production of oxygen was observed. In all cases 
there was an increase in oxygen pressure beyond the 
amount due to normal cell activity. Figs. 1 and 2 show 
the results of two typical experiments. In each case a 
control run was made parallel to the test. Into the control 
there was placed, at the same time as the catalase solution 
was added to the test, an equivalent amount of distilled 
water. The solid curve in Fig. 1 represents this control. 
In Fig. 2 we see the results of adding catalase solution 
first to one cell suspension and later to the other. Attention 
is called to the fact that, at the end, the difference in 
oxygen pressures was approximately the same as at the 
start although there was a large difference before the 
addition of catalase to the second cell suspension. The ex- 
periments were run at 25+0.02 degrees centigrade. 

Similar experiments we performed with solutions of 
sodium magnesium chlorophyllin gave like sharp dis- 
continuities in oxygen pressure on the addition of catalase 
solution and can be interpreted as a duplication of photo- 
synthesis in vitro without the intervention of living cells. 

Therefore, chlorophyll, in the plant, can be supposed to 
fill the dual role of light absorber and energy utilizer and, 
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Fic. 1. Effect of catalase on oxygen production during photo- 
synthesis. The broken lines represent oxygen production after the 
addition of catalase solution. 
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Fic. 2. Effect of catalase on oxygen production during photo- 
synthesis. The broken lines represent oxygen production after the 
addition of catalase solution. Both systems were opened to and equili- 
brated with air at the end of two and six minutes in order to add 
catalase and water. The curves are corrected for the lag. 


in conjunction with the fat-protein to which it is con- 
nected, as an enzyme which decomposes hydrogen peroxide. 
A complete account of the experimental methods and a 
fuller presentation of the theory will follow in a short time. 
ALEXANDER F, KNOLL 
FREDERIC L. MATTHEWS, JR. 
Ray H. Crist 


Columbia University, 
New York, New York, 
January 15, 1938. 


* Publication assisted by the Ernest Kempton Adams Fund for 
Physical Research of Columbia University. 

1 Emerson and Arnold, J. Gen. Phys. 15, 391 (1932). 

2 Emerson and Arnold, J. Gen. Phys. 16, 191 (1932). 

3 Franck and Herzfeld, J. Chem. Phys. 5, 237 (1937). 
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Infrared Absorption Spectrum of Heavy Acetic Acid 


The infrared absorption spectrum of acetic acid-d 
(CH;COOD) vapor has been investigated in the region 
1-15.54 by means of a rocksalt spectrometer, previously 
described in the literature.! This investigation was under- 
taken with the aim of bringing out the hydrogen bond 
frequency by a deuterium substitution. It is well known 
that acetic acid vapor at room temperature exists mostly 
as double molecules, the bonding occurring through the 
mechanism of ‘‘hydrogen bridges.” According to current 
ideas,? the hydrogen bonded O—H frequency is shifted 
from 2.75 into the region of the C—H valence vibrations 
at 3.34 and cannot be clearly discerned. In heavy acetic 
acid this frequency will be shifted approximately by a 
factor of v2 and will lie in a region usually unoccupied by 
absorption bands (4.34). That this is so, is clearly brought 
out by the curves in Fig. 1, and Table I which gives the 


TABLE I. 


1389+10 cm™ 
1317+10 
1250+10 
1056+ 5 


1000+ 5 cm 
840+ 5 
667 +20 


3079+40 cm7! 
2325+20 
2105+15 
1725+15 


positions of the band centers. It is of interest to note that 
the C—H valence frequency appears only weakly. This 
fact suggests a possible interpretation of the broad and 
deep absorption at 3.25 in light acetic acid as being chiefly 
due to the hydrogen bond vibration. 

Raman spectra of liquid CH;COOD#:* show no lines 
in the region 1655 to 2940 cm“, thus making impossible 
a complementary study in the Raman effect. 

Our sample of heavy acetic acid was prepared by reacting 
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Fic. 1. Absorption spectra of CHsCOOD. Curve A, pressure 15.5 
mm, temperature 24°C; curve B, pressure 4 mm, temperature 24°C. 
Absorption cells 15 cm in length. 


pure acetyl chloride with 99.6 percent heavy water by a 
method previously employed by Engler.‘ 

Further experiments are now in progress to determine 
how association through the deuterium bond is affected 
by temperature. 

We wish to thank Mr. R. R. Brattain for his encourage- 
ment and very kind assistance. 

R. C. HERMAN 


Princeton University, R.H 
. HOFSTADTER 


Princeton, New Jersey, 
December 17, 1937. 
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3W. R. Angus, A. H. Leckie and C. L. Wilson, Proc. Roy. Soc. 155, 
183 (1936). 

4W. Engler, Zeits. f. physik. Chemie B32, 471 (1936). 





